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ABSTRACT 


'A  mathematical  simulation  model  of  a  coastal  storm  has  been 
programmed  to  forecast  or  hindcast  wave  and  longshore  current  con¬ 
ditions  at  a  coastal  site.  Storm  parameters  for  the  model  are 
based  on  the  size,  shape  intensity  and  path  of  the  storm  as  de¬ 
rived  from  weather  maps.  An  elliptical  form  is  used  to  model  the 
size  and  shape  of  the  storm  which  are  controlled  by  varying  the 
length  and  orientation  of  the  major  and  minor  axes.  Storm  in¬ 
tensity  is  a  function  of  the  barometric  pressure  gradient  which 
is  modeled  by  an  inverted  normal  curve  through  the  storm  center. 
The  storm  path  is  based  on  actual  storm  positions  for  the  hind- 
cast  mode,  and  on  projected  positions  assuming  constant  speed 
and  direction  for  the  forecast  mode.  The  location,  shoreline 
orientation  and  nearshore  bottom  slope  provide  input  data  for 
each  coastal  site. 

For  each  storm  position,  the  geostrophic  wind  speed  and  di¬ 
rection  are  computed  at  the  shore  site  as  a  function  of  baro¬ 
metric  pressure  gradient  and  latitude.  The  geostrophic  wind  is 
converted  into  surface  wind  speed  and  direction  by  applying 
corrections  for  frictional  effects  over  land  and  sea.  fThe  sur¬ 
face  wind  speed,  fetch  and  duration  are  used  to  computd  the  wave 
period,  breaker  height  and  breaker  angle  at  the  shore  slite.  The 
longshore  current  velocity  is  computed  as  a  function  of  wave 
period,  breaker  height  and  angle  and  nearshore  bottom  slope. 


The  model  was  tested  by  comparing  hindcast  output  with  ob¬ 
served  data  for  several  coastal  locations.  Forecasts  were  made 
for  actual  storms  and  for  hypothetical  circular  and  elliptical 
shaped  storms. 
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COASTAL  STORM  MODEL 


INTRODUCTION 


Coastal  storms  which  provide  a  combination  of  high  winds, 
pounding  waves  and  rapid  longshore  currents  are  a  major  cause  of 
distructive  erosion  along  beaches  and  cliffs.  Beaches  which  are 
generally  composed  of  sand  or  cobbles  are  subject  to  sudden  changes 
during  storms.  During  one  storm  at  Chesil  Bank  near  Abbotsbury, 
England,  the  crest  of  a  shingle  beach  was  cut  back  1.53  meters  in 
3  hours  (Lev;is,  1931).  During  a  severe  storm  in  July  1969  at 
Stevensville,  Michigan,  the  beach  and  bluff  were  eroded  back  over 
5.5  meters  when  the  waves  reached  a  height  of  2  meters  (Fox  and 
Davis,  1970b).  On  the  Oregon  coast,  a  beach  was  stripped  of  a 
2  meter  thick  blanket  of  sand  and  the  wave  cut  terrace  was  exposed 
when  the  waves  reached  heights  of  8  meters  during  late  Novemeber 
storms  (Fox  and  Davis,  1974). 

For  any  operation  involving  the  coastal  zone,  it  is  essential 
to  make  predictions  of  wave  and  current  conditions  during  a  coastal 
storm.  General  wave  forecasts  on  a  worldwide  grid  are  available 
from  the  National  Weather  Service  and  Fleet  Numerical  Weather  Center. 
These  forecasts  provide  accurate  predictions  of  wave  conditions  on 
the  open  ocean,  but  do  not  provide  detailed  enough  predictions  for 
the  coastal  zone.  Therefore,  a  computer  simulation  model  was  de¬ 
veloped  to  fill  the  gap  between  wave  predictions  on  the  open  ocean 
and  surf  predictions  along  the  coast. 

The  coastal  storm  model  utilizes  an  ellipse  to  simulate  a  map 
of  barometric  pressure.  The  shape  of  the  storm  can  be  modified 
by  varying  the  length  and  orientation  of  the  major  and  minor  axes 
of  the  ellipse.  The  intensity  of  the  storm  is  controlled  by  in¬ 
creasing  or  decreasing  the  range  in  barometric  pressure.  The  actual 
storm  track  and  dimensions  are  read  in  as  data  for  hindcasting  wave 
and  current  conditions.  For  making  forecasts,  the  size,  shape  and 
intensity  of  the  storm  are  provided  as  input  data  for  the  model.  In 
forecasting,  the  storm  path  is  plotted  by  assuming  a  constant  azimuth 
and  velocity. 

One  of  the  initial  steps  in  developing  a  coastal  storm  model  is 
determination  of  the  barometric  pressure  gradient  at  any  point  on  the 
ground  surface  under  the  storm.  The  pressure  gradient  is  then  used 
in  conjunction  with  the  latitude  to  calculate  the  geostrophic  wind 
speed,  which  in  turn  is  used  to  compute  the  surface  wind  speed,  wave 
height  and  longshore  current  velocity.  In  the  model,  it  is  assumed 
that  a  profile  of  barometric  pressure  along  the  major  or  minor  axis 
of  the  storm  ellipse  can  be  represented  by  an  inverted  normal  curve. 
By  rotating  the  normal  curve  around  the  ellipse  and  taking  the  deriv¬ 
ative,  it  is  possible  to  calculate  the  pressure  gradient  at  any  point 
on  the  ground.  From  that  point  on,  conventional  methods  are  employed 
for  determining  the  geostrophic  wind  speed,  surface  wind  speed,  wave 
height  and  longshore  current  velocity. 
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PREVIOUS  WORK 


The  coastal  storm  model  is  based  on  a  series  of  detailed  field 
studies  which  extended  from  1969  through  1975.  The  studies  inclu¬ 
ded  the  analysis  and  synthesis  time-series  data  on  barometric  pres¬ 
sure,  wind  speed  and  direction,  wave  period  and  height  and  longshore 
current  velocity  for  15  days  to  1  year.  Topographic  profiles  across 
the  beach  and  nearshore  area  were  used  to  construct  topographic  maps 
and  maps  of  erosion  and  deposition.  The  sites  in  the  study  are 
plotted  by  date  in  Figure  1  and  included  in  the  following  list  along 
with  references  for  each  study. 


1969 

Stevensville,  Michigan 

Fox  and  Davis,  1970a,  1970b 

1970 

Holland,  Michigan 

Fox  and  Davis,  1971a 

Davis  and  Fox,  1971 

1971-2 

Mustang  Island,  Texas 

Davis  and  Fox,  1972c 

Davis  and  Fox,  1975 

1972 

Sheboygan,  Wisconsin 

Fox  and  Davis,  1972 

1973 

Cedar  Island,  Virginia 

Davis  and  Fox,  1974a 

1973-4 

South  Beach,  Oregon 

Fox  and  Davis,  1974 

1974 

Zion,  Illinois 

Davis  and  Fox,  1974b 

1974 

South  Haven,  Michigan 

Davis  and  Fox,  1974b 

1975 

Plum  Island,  Massachusetts 

The  models  have  evolved  from  a  geometric  model  called  the 
area-time  prism  (Davis  and  Fox,  1972a)  through  a  conceptual  model 
(Davis  and  Fox, 1972b)  to  a  process-response  model  for  Lake  Michigan 
(Fox  and  Davis,  1971b  and  1973  ).  The  simulation  model  developed 
for  Lake  Michigan  was  limited  to  the  local  geographic  area  where 
the  storms  moved  directly  onshore  to  the  north  of  the  study  area. 

The  proposed  coastal  storm  model  is  an  outgrowth  of  the  earlier 
model  but  is  more  generalized  with  broader  application  under  a 
wider  range  of  storm  conditions  and  shoreline  orientations. 

Several  different  types  of  computer  models  have  been  proposed 
for  the  coastal  environment.  Probablistic  models  were  developed  to 
reproduce  gross  coastal  features  such  as  a  recurved  spit  on  the 
south  coast  of  England  (McCallagh  and  King,  1970)  and  the  Mississippi 
River  Delta  (McCammon,  1971).  A  markcv  process  was  used  to  simu¬ 
late  the  sequence  of  bar  formation  and  migration  across  a  beach 


(Sonu  and  van  Beek.  1971).  A  deterministic  model  resembling  a 
wave  tank  experiment  was  proposed  to  simulate  the  interaction 
between  a  prograding  delta  and  waves  (Komar,  1973).  A  statis¬ 
tical  model  with  a  relatively  simple  beach  topography  to  compute 
breaker  height,  longshore  current  velocity  and  wave  setup  (Collins, 
1971)  was  followed  by  a  more  deterministic  approach  to  model  the 
nearshore  circulation  patterns  employing  monochromatic  waves  and 
more  complex  beach  topographies  (Noda,  et  al_,  1974).  An  explicit 
finite  difference  model  for  predicting  time-dependent,  wave  in¬ 
duced  nearshore  circulation  was  developed  by  Birkemeier  and 
Dalrymple  (1976). 

On  a  larger  scale,  Resio  and  Hayden  (1973)  proposed  an  inte¬ 
grated  storm  model  which  combines  three  scales  of  atmospheric  mo¬ 
tion,  large  scale,  synoptic  scale  and  small  scale  into  an  estima¬ 
tion  of  a  winter  wave-surge  climate  for  the  mid-Atlantic  coast. 

At  a  similar  scale.  Goldsmith,  et_  a_l_  (1974)  developed  a  wave  cli¬ 
mate  model  for  the  mid-Atlantic  coa:t  by  using  Dobson's  (1967) 
wave  refraction  program  to  project  offshore  waves  into  the  coastal 
zone. 


The  coastal  storm  model  proposed  in  this  report  provides  a 
link  between  the  large-scale,  seasonal  wave-climate  models  and 
the  dynamic  surf  zone  models.  By  tracking  a  storm  across  a  shore¬ 
line,  the  wave  parameters  which  are  output  from  the  storm  model 
furnish  input  for  the  surf  zone  models.  Therefore,  the  proposed 
storm  model  could  be  combined  with  other  computer  models  to  pro¬ 
vide  an  integrated  process  model  for  the  coastal  zone. 


4 


COMPUTER  PROGRAMS  -  COASTAL  STORM  MODEL 


Program  STORM 

Program  STORM  is  a  mathematical  simulation  model  which  has 
been  programmed  for  the  computer  to  forecast  or  hinacast  wave 
conditions  at  a  coastal  site  during  a  storm.  The  actual  storm 
as  represented  by  the  isobars  on  a  weather  map  is  modeled  by 
an  elliptical  storm  with  major  and  minor  axes  at  right  angles 
and  passing  through  the  center  of  the  low.  The  size,  shape, 
intensity  and  path  of  the  storm  as  determined  from  weather  maps 
are  used  to  generate  the  surface  wind  pattern,  wave  height  and 
period,  and  longshore  current  velocity  as  the  storm  moves  across 
the  coast. 

The  computer  program  is  divided  into  a  main  program,  STORM, 
and  a  series  of  11  subroutines.  The  main  program  is  used  to  read 
in  the  data,  call  the  various  subroutines  for  computing  the  wind, 
wave  and  current  conditions,  print  out  the  predictions  at  one  hour 
intervals.  All  the  input  and  output  is  handled  by  the  main  pro¬ 
gram  while  the  calculations  are  carried  out  by  the  various  sub¬ 
routines.  In  this  way,  any  portion  of  the  model  can  be  indepen¬ 
dently  tested  by  using  a  small  main  program  to  call  each  subrou¬ 
tine  individually.  Therefore,  if  any  problem  arises  in  the  main 
program,  it  can  be  narrowed  down  to  a  particular  subroutine,  and 
that  subroutine  can  be  tested  under  a  variety  of  conditions  with¬ 
out  using  the  main  program.  Also,  if  a  portion  of  the  simulation 
model  is  to  be  incorporated  into  another  program,  any  of  the  indi¬ 
vidual  subroutines  can  be  removed  and  used  separately  with  the 
appropriate  calling  arguments. 

The  theory  and  mechanics  of  the  program  will  be  explained  in 
detail  starting  with  the  MAIN  program  and  proceeding  through  each 
of  the  subroutines  as  they  are  called  by  the  MAIN  program.  The 
program  was  written  in  FORTRAN  IV  for  an  8k  IBM  1130  at  Williams 
College.  A  full  listing  of  the  programs  with  appropriate  comment 
cards  is  included  in  Appendix  A.  A  second  version  of  the  MAIN  pro¬ 
gram  was  written  for  the  Xerox  530  which  includes  a  graphics  pack¬ 
age  for  a  29  inch  plotter.  The  graphics  package  is  used  to  plot 
barometric  pressure,  surface  wind,  onshore  wind,  alongshore  wind, 
breaker  height,  wave  period  and  longshore  current. 
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The  main  program  is  used  to  read  the  input  data  for  the  storm 
and  shoreline  conditions,  call  the  various  subroutines  and  print 
out  the  results  at  one  hour  intervals.  A  listing  of  the  input  cards 
is  included  in  Appendix  I  with  a  description  of  each  of  the  input 
variables.  The  program  is  dimensioned  to  make  predictions  up  to 
130  hours  or  5  days  and  10  hours.  If  a  longer  prediction  is  desired, 
it  is  necessary  to  increase  the  dimension  of  U (1 30 )  and  V ( 1 30 )  to 
the  required  number  of  hours.  In  the  model,  130  hours  was  selected 
because  of  core  limitations  on  the  8k  IBM  1130.  For  most  of  the 
storms,  the  130  hour  limitation  was  not  a  serious  constriction, 
however,  a  larger  dimension  statement  would  be  helpful  in  some  cases. 

The  first  two  data  cards  are  used  to  read  in  the  title,  starting 
time,  date  and  input/output  options.  The  title  used  for  the  location 
of  the  shore  site  can  be  up  to  80  spaces  long  filling  card  1.  On 
the  second  card,  the  starting  hour,  ISTRT,  is  included  in  columns  1 
and  2  followed  by  the  date,  DAY,  in  columns  3  to  22.  The  starting 
time  is  read  in  as  an  interger  and  must  be  right  justified.  If  ISTRT 
is  read  in  as  0,  the  program  will  terminate.  For  the  input  option, 
INAUT,  in  column  23,  a  0  is  used  for  metric  units  and  a  1  is  used 
for  nautical  units  including  nautical  miles,  knots,  and  feet  per  sec¬ 
ond.  The  output  option,  NAUT,  in  column  24  is  separate  from  the  in¬ 
put  option  but  uses  the  same  code,  0  for  metric  units  and  1  for  nau¬ 
tical  units.  Although  metric  units  are  becoming  the  standard  and 
are  now  required  for  scientific  reports,  it  may  be  desirable  at  times 
to  have  the  input  or  output  in  nautical  units.  With  separate  input 
and  output  options,  it  is  possible  to  have  the  input  in  metric  or 
nautical  units  and  convert  to  the  other  units  with  the  output. 

The  first  3  columns  of  card  3  are  used  to  select  the  major 
options  for  the  program.  For  the  first  option,  IN0PT,  a  1  in  column 
1  will  call  the  hindcasting  mode,  while  a  2  in  column  1  will  call 
the  forecasting  mode.  The  hindcasting  mode  is  used  when  storm  posi¬ 
tions  are  available  at  six  hour  intervals.  The  hourly  positions  of 
the  storm  are  determined  by  a  linear  interpolation  between  the  6-hour 
positions.  For  the  forecasting  mode,  the  initial  position  of  the  storm 
along  with  a  constant  velocity  and  azimuth  are  used  to  calculate 
successive  positions  at  1-hour  intervals.  The  variables  for  the 
hindcasting  and  forecasting  modes  are  read  in  on  card  6.  The  second 
option  on  card  3  is  the  tide  prediction  option,  IFTID.  If  a  0  is 
punched  in  column  2,  the  tide  prediction  option  is  suppressed  and 
card  4  is  not  included  in  the  data  set.  The  tide  option  is  omitted 
for  a  non-tidal  body  of  water,  such  as  the  Great  Lakes.  Where  tide 
data  are  available  from  the  tide  tables  or  from  observations,  a  1 
is  punched  in  column  2,  and  the  tide  data  are  included  on  card  4. 

The  longshore  current  equation  for  the  simulation  run  is  selected 
in  option  3,  LSC0P.  Four  different  longshore  current  equations  are 
included,  (1)  Fox  and  Davis  (1972),  (2)  Longuet-Higgins  (1970),  (3) 


Coastal  Engineering  Research  Center  (1973),  and  (4)  Komar  and 
Inman  (1970).  The  longshore  current  equations  are  called  in 
subroutines  SURF  and  their  differences  will  be  discussed  under 
that  subroutine. 

The  number  of  storm  positions,  NX,  are  punched  in  columns 
4  to  6  of  card  3,  for  6-hour  intervals  in  the  hindcasting  mode, 
and  for  1-hour  intervals  in  the  forecasting  mode.  For  example, 
if  the  hindcasting  mode  is  used,  and  3  days  of  data  are  included, 
the  initial  position  and  4  positions  for  each  day  would  give  a 
total  of  13  for  NX.  For  the  forecasting  mode,  a  3  day  forecast 
would  use  a  73  for  NX,  1  for  the  initial  position,  and  72  for  the 
72  hour  forecast.  The  maximum  value  for  NX  is  22  for  the  hind- 
cast  mode  and  130  for  the  forecast  mode. 

The  average  basin  fetch  in  kilometers,  BNFCH,  is  punched  in 
columns  7  to  12  on  card  3.  The  average  basin  fetch  is  used  as 
the  limiting  fetch  in  determining  the  wave  height  and  period  from 
the  wind  speed.  Where  the  basin  fetch  is  smaller  than  the  maximum 
storm  axis,  the  waves  are  fetch  limited.  However,  where  the  fetch 
is  significantly  larger  than  the  storm  size,  the  average  basin  fetch 
will  not  be  a  limiting  factor  in  determining  the  wave  parameters. 

The  basin  fetch  is  considered  in  an  offshore  direction  from  the  shore 
site.  In  the  case  of  a  large  ocean,  the  approximate  width  of  the 
ocean  can  be  used  as  the  basin  fetch. 

In  columns  13  to  17  of  card  3,  the  time  interval  between  storm 
positions  TINT  is  normally  set  at  1.0.  The  time  interval  refers  to 
the  printout  spacing  for  the  forecast  modes.  For  the  hindcast  mode, 
the  values  are  read  in  at  6-hour  intervals,  and  the  results  are 
printed  out  at  1-hour  intervals. 

The  minimum  barometric  pressure  in  millibars,  PMIN,  taken  at 
the  center  of  the  low  pressure  cell  is  punched  in  columns  18  to  24 
on  card  3.  Usually,  the  minimum  barometric  pressure  is  interpolated 
within  the  smallest  isobar.  Thus,  if  the  smallest  isobar  is  1004, 
and  the  isobar  spacing  is  4  millibars,  the  minimum  pressure  would 
be  estimated  at  1002  millibars.  The  pressure  at  the  largest  en¬ 
circling  isobar,  PMAXR,  is  used  to  determine  the  intensity  of  the 
storm.  If  the  storm  is  circular  or  oval  shaped,  the  largest  isobar 
which  encloses  the  storm  center  is  used  for  PMAXR  and  punched  in  co¬ 
lumns  25  to  31  of  card  3.  If,  however,  the  storm  has  a  wave  form 
extending  down  from  the  north, a  line  is  drawn  along  the  storm  path 
through  the  storm  center  to  the  margins  of  the  storm.  The  largest 
isobar  which  the  line  crosses  on  both  sides  of  the  storm  is  then 
considered  the  largest  encircling  isobar,  PMAXR.  In  the  program, 
the  largest  encircling  isobar  is  defined  as  2  standard  deviations 
away  from  the  center  of  the  storm.  Therefore,  the  total  storm  radius 
would  be  1.5  times  the  radius  of  the  largest  encircling  isobar.  The 
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pressure  range  would  be  1.145  times  the  range  within  the  largest 
encircling  isobar.  The  latitude  at  the  shore  site,  SLAT,  punched 
in  columns  32  to  36  is  used  in  subroutine  WIND  to  compute  the 
geostrophic  wind  speed. 

The  geographic  size  of  the  storm  is  defined  in  terms  of  an 
ellipse  with  a  major  half-axis  and  minor  half-axis  corresponding 
to  radius  of  a  circle.  The  major  half-axis,  AR,  of  the  storm 
ellipse  measured  when  the  storm  is  closest  to  the  study  site  is 
punched  in  columns  37  to  42.  If  the  storm  ellipse  is  assymetrical , 
the  longest  half-axis  on  the  side  toward  the  shore  location  is  used 
as  the  major  half-axis.  The  minor  half-axis,  BR,  is  measured  at 
right  angles  to  the  major  half-axis  through  the  center  of  the  low. 

The  major  and  minor  half-axes  are  measured  from  the  center  of  the 
low  to  the  largest  encircling  isobar  PMAXR.  The  minor  half-axis, 

BR,  is  punched  in  columns  43  to  48  on  card  3. 

The  orientation  of  the  major  half  axis  EAZ  is  punched  in  col¬ 
umns  49  to  54  of  card  3.  The  orientation  of  the  major  axis  is 
measured  in  degrees  from  true  north  to  the  northern  end  of  the 
major  axis  ranging  from  -90°  on  the  west  to  90°  on  the  east.  For 
a  front  or  trough  related  to  a  low  pressure  system,  the  major  axis 
is  usually  several  times  longer  than  the  minor  axis  and  the  orien¬ 
tation  of  the  major  axis  lies  along  the  line  of  the  front.  In  a 
circular  or  oval  storm,  the  major  axis  is  usually  1  to  1.5  times 
as  long  as  the  minor  axis.  The  major  and  minor  half  axis  are  mea¬ 
sured  when  the  center  of  the  storm  is  at  its  nearest  position  to 
the  shore  study  site. 

Variables  for  hourly  tide  prediction  are  contained  on  card  4. 

The  spring  tide  range  in  meters,  ST,  is  punched  in  columns  0  to  5 
and  the  neap  tide  range,  TN,  in  columns  6  to  10.  The  spring  and 
neap  tide  ranges  and  the  time  of  the  last  spring  high  tide  are  a- 
vailable  in  the  tide  tables  which  are  published  annually  by  N.O.A.A. 
In  making  the  hourly  predictions  for  the  model,  it  is  necessary  to 
punch  the  number  of  days  since  the  last  spring  tide,  TDAY,  in  col¬ 
umns  11  to  15.  The  hour  of  the  last  spring  high  tide,  THR,  preceding 
the  start  of  the  run  is  punched  in  columns  16  to  20.  The  tidal  form 
number,  FN,  punched  in  columns  12  to  25  is  used  to  reproduce  a  semi¬ 
diurnal,  mixed-semidiurnal,  mixed-diurnal  or  diurnal  tide  with  the 
right  spacing  and  tidal  beat.  The  nearshore  bottom  slope  at  low 
tide,  SLPLO,  is  punched  in  columns  26  to  32  and  the  slope  at  high 
tide,  SLPHI,  is  punched  in  columns  33  to  39.  The  nearshore  bottom 
slope  which  varies  with  tidal  elevation  is  used  for  computing  long¬ 
shore  current  velocity  in  subroutine  SURF.  The  preferred  method  of 
determining  bottom  slope  is  to  fit  a  linear  surface  to  the  nearshore 
map  at  low  tide,  and  repeat  the  process  at  high  tide.  The  linear 
slope  should  extend  to  a  depth  of  at  least  twice  the  breaker  height. 
For  the  high  tide  range,  the  foreshore  slope  and  low  tide  terrace 
should  be  included  in  the  slope  calculation.  Where  it  is  not  pos¬ 
sible  to  fit  a  linear  surface  because  of  lack  of  data,  it  is  possible 


to  get  an  approximation  of  the  nearshore  slope  by  measuring  the 
depth  at  some  predetermined  distance  from  the  shore  at  low  tide 
and  at  high  tide.  By  dividing  the  depth  by  the  distance,  a  good 
approximation  of  nearhsore  slope  can  be  estimated  for  low  and  high 
tides.  The  nearshore  slope  is  an  initial  factor  for  the  deter¬ 
mination  of  the  longshore  current  velocity,  so  care  should  be 
taken  in  estimating  nearshore  slope  at  low  and  high  tide.  It  also 
should  be  pointed  out  that  nearshore  bars  have  a  significant  in¬ 
fluence  on  nearshore  currents  and  must  be  considered  in  making  an 
estimate  of  the  nearshore  slope.  The  final  variable  on  the  tide 
prediction  card  is  the  mean  tide  level,  TMEAN,  which  is  the  dif¬ 
ference  between  mean  sea  level  and  mean  low  tide  as  reported  in 
the  annual  tide  tables. 

Data  for  the  shore  site  location  including  geographic  coor¬ 
dinates,  onshore  direction,  average  bottom  slope  and  offshore  is¬ 
land  option  are  punched  on  card  5.  The  shore  site  location  is 
given  in  a  X  and  V  coordinate  system  where  the  X-axis  runs  east- 
west  with  positive  X  in  the  east  direction,  and  the  Y-axis  north- 
south  with  positive  Y  in  the  north  direction.  The  X-Y  coordinate 
system  is  measured  in  kilometers  with  the  origin  located  at  the 
southwest  of  the  study  site.  In  practice,  a  piece  of  10  to  the 
inch  rectangular  grid  graph  paper  is  laid  over  the  weather  map 
with  the  Y  axis  parallel  to  the  longitude  line  nearest  the  study 
site.  The  origin  of  the  graph  paper  is  placed  several  inches  to 
the  southwest  of  the  shore  location  so  that  the  X-axis  runs  east- 
west  and  the  Y-axis  runs  north-south  parallel  to  the  latitude  and 
longitude  lines  through  the  study  site.  The  X  and  Y  coordinates 
are  read  off  the  map  in  inches  and  converted  to  kilometers  before 
they  are  used  in  the  program.  The  X  coordinate,  ULOC,  is  punched 
in  columns  1  to  7  and  the  Y  coordinate,  VLOC,  is  punched  in  col¬ 
umns  8  to  14  on  card  5. 

The  orientation  of  the  shoreline  given  by  the  onshore 
azimuth,  SHAZ,  and  the  average  nearshore  bottom  slope,  SLOPE,  are 
punched  on  columns  15  to  21  and  22  to  38  respectively.  The  onshore 
direction  measured  in  degrees  in  a  clockwise  direction  from  north 
is  used  to  give  the  orientation  of  the  shoreline.  Since  the  storm 
is  considered  a  regional  feature,  it  is  necessary  to  give  the  re¬ 
gional  orientation  for  the  shoreline.  An  east-west  shoreline  with 
land  to  the  north  would  give  a  0  azimuth.  If  a  shoreline  is  running 
north-south  with  the  land  to  the  east  and  water  to  the  west,  the 
shoreline  azimuth  would  be  90  degrees.  Similarlly,  if  the  shoreline 
is  running  north-south  with  the  land  on  the  west  and  the  water  on 
the  east,  the  onshore  azimuth  would  be  270  degrees. 

An  option  is  available  with  the  program  for  an  offshore  island 
which  is  not  influenced  by  a  large  continental  land  mass.  When  a  0 
is  punched  in  column  30  of  card  5,  the  offshore  island  option  is  sup¬ 
pressed  and  a  normal  continental  coast  or  barrier  island  is  assumed. 
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For  a  coast  backed  by  land,  land  corrections  are  used  in  computing 
the  surface  wind  speed  when  the  wind  blows  offshore.  Therefore, 
when  the  island  option  is  used  and  a  1  is  punched  in  column  30  of 
card  5,  the  wind  is  assumed  to  be  blowing  from  the  sea  in  all 
directions  and  the  land  correction  is  not  used.  For  a  barrier 
island  which  lies  roughly  parallel  to  the  coast,  the  island  option 
is  not  used  because  an  offshore  wind  blows  over  land  for  a  long 
distance  before  it  hits  the  lagoon  and  barrier  island. 

The  storm  positions  for  the  hindcasting  and  forecasting  modes 
are  punched  on  card  6.  For  option  1,  the  hindcasting  mode,  the  X 
and  Y  coordinates  are  punched  in  columns  1  to  7  and  8  to  14  respec¬ 
tively.  The  X  and  Y  coordinates  are  from  the  rectangular  grid  dis¬ 
cussed  for  the  shore  site  location  on  card  5.  The  coordinates  for 
the  storm  are  given  for  the  initial  storm  position  and  at  successive 
6  hour  intervals,  with  one  pair  of  coordinates  per  card.  The  num¬ 
ber  of  pairs  of  coordinates  is  specified  by  NX,  the  number  of  storm 
positions  on  card  3.  For  the  forecasting  mode,  option  2,  the  storm 
positions  are  determined  at  1  hour  intervals  from  the  storm  velocity, 
storm  azimuth  and  initial  X  and  Y  coordinates.  The  storm  velocity, 
SVEL,  in  columns  1  to  7,  is  given  in  kilometers  per  hour.  Reason¬ 
able  storm  velocities  would  vary  from  about  25  to  75  kilometers  per 

hour  for  a  slow  to  fast  moving  storm.  In  the  forecasting  mode,  it 

is  not  possible  to  vary  the  storm  velocity,  so  the  initial  storm 
velocity  must  be  maintained  for  the  entire  forecast  run.  The  storm 
azimuth  or  path  is  measured  in  degrees  clockwise  from  north.  As 
with  the  storm  velocity,  it  is  not  possible  to  vary  the  storm  azi¬ 
muth  in  the  middle  of  a  forecast  run.  An  azimuth  of  0  degrees  would 

have  the  storm  moving  due  north,  and  a  90  degree  azimuth  would  have 

the  storm  heading  east.  In  the  forecast  mode,  the  initial  X  and  Y 
coordinates  for  the  storm  are  punched  in  columns  15  to  21  and  22  to 
28  respectively.  It  is  possible  to  make  a  map  for  each  predicted 
variable  by  making  a  series  of  runs  with  different  initial  coordin¬ 
ates. 


It  is  possible  to  run  a  series  of  models  for  different  coastal 
situations  by  including  a  new  data  set  for  each  model  starting  with 
the  title  card.  To  terminate  the  run,  two  blank  cards  are  included 
at  the  back  of  the  data  deck.  Since  the  second  card  of  the  new  data 
set  is  blank,  ISTRT  is  read  in  as  0  and  the  program  will  finish. 

Different  versions  of  the  main  program  are  used  for  making  fore¬ 
casts  directly  from  the  console,  and  for  printing  a  map  using  the 
forecast  mode.  Listing  of  the  programs  are  Included  in  Appendix  A 
along  with  explanations  of  the  input  options. 
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Subroutine  LOCAT 


As  a  storm  moves  across  a  coastline,  subroutine  LOCAT  is 
used  to  determine  the  position  of  the  coastal  site  relative  to 
the  storm  center  for  each  increment  of  time.  In  preparing  the 
input  data  for  the  program,  the  location  of  the  coastal  site  and 
a  sequence  of  storm  positions  are  plotted  on  a  rectangular  grid 
referred  to  as  the  map  coordinate  system.  For  the  map  coordinate 
system,  the  X  axis  points  east,  the  Y  points  heading  north  and 
the  origin  is  located  to  the  southwest  of  the  Initial  storm  po¬ 
sition.  When  the  shore  location  and  storm  positions  are  plotted 
on  a  weather  map,  the  coordinates  are  measured  in  kilometers  or 
nautical  miles,  whichever  are  the  most  convenient  units  for  the 
project. 

In  computing  the  geostrophic  wind  speed,  it  is  necessary  to 
determine  the  gradient  in  barometric  pressure  at  the  coastal  site. 
Therefore,  a  storm  coordinate  system  is  established  with  the  ori¬ 
gin  at  the  center  of  the  storm,  the  XI  axis  parallel  to  the  shore, 
the  Y1  axis  perpendicular  to  the  shore,  and  the  positive  Y1  di¬ 
rection  heading  onshore.  When  facing  the  land  from  the  sea,  the 
positive  XI  direction  is  along  the  shore  to  the  right,  and  the 
negative  XI  direction  is  to  the  left  (Figure  2).  Each  time  the 
storm  moves,  the  origin  of  the  storm  coordinate  system  Is  also 
moved  to  the  new  location  for  the  center  of  the  storm.  However, 
the  orientation  of  the  XI  and  Y1  axes  remains  the  same  with  the 
XI  axis  parallel  to  the  shore  and  the  Y1  axis  at  right  angles  to 
the  coast. 

In  the  storm  model,  the  units  for  the  XI,  Y1  coordinate  sys¬ 
tem  are  converted  from  kilometers  or  nautical  miles  to  storm  radii 
by  dividing  by  the  radius  of  the  storm.  In  an  elliptical  storm,  the 
length  of  the  major  half  axis  Is  used  in  place  of  the  storm  radius. 

In  subroutine  LOCAT,  ULOC  and  VLOC  are  the  map  coordinates  for 
the  coastal  site,  and  UST  and  VST  are  the  map  coordinates  for  the 
storm  center  (Figure  2).  Vectors  are  computed  parallel  to  the  X 
axis  (U  =  UST-ULOC)  and  parallel  to  the  Y  axis  (V  =  VLOC-VST).  The 
resultant  vector  U?  +  V  )  gives  the  map  distance  from  the  storm 
center  to  the  coastal  site.  The  counterclockwise  angle  (ANG)  be¬ 
tween  the  positive  X  axis  and  the  Z  vector  is  computed  by  the  arc¬ 
tangent  subroutine  ARCTA.  The  onshore  azimuth  (SHAZ)  is  the  onshore 
direction  normal  to  the  shoreline  measured  in  a  clockwise  direction 
from  north.  The  angle  A  (A  =  ANG  -  SHAZ)  is  used  for  converting 
coordinates  from  the  map  system  to  the  storm  system.  The  shore 
position  is  then  determined  In  the  storm  coordinate  system  for  dis¬ 
tances  along  the  XI  axis  (XI  =  -Z  *  cos  (A))  and  along  the  Y1  axis 
(Y1  -  Z  *  sin  (A)). 

A  third  coordinate  system  Is  set  up  for  dealing  with  an  ellip¬ 
tical  storm.  In  the  elliptical  coordinate  system,  the  P  axis  lies 
along  the  major  half-axis  and  the  Q  axis  lies  along  the  minor  half- 


11 


axis.  The  distances,  P  and  Q,  within  the  storm  ellipse  are  used 
to  locate  the  shore  site  relative  to  the  center  of  the  ellipse. 
The  orientation  of  the  storm  ellipse  is  given  by  the  ellipse  azi¬ 
muth  (EAZ)  which  Is  the  azimuth  of  the  major  half-axis  plus  or 
minus  90  degrees  from  true  north. 


A 


Figure  2.  A.  Map  coordinate  system  (X-Y)  for  locating  storm 

center  and  shore  site,  and  B.  Storm  coordiante  system 
(Xl-Yl)  with  origin  at  storm  center  and  X  axis  parallel 
to  the  shore. 
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Subroutine  EL  I PS 

The  wind  angle  and  barometric  pressure  gradient  at  any 
position  within  an  elliptical  storm  are  determined  by  subrou¬ 
tine  ELIPS.  A  is  the  length  of  the  major  half  axis,  B  is  the 
length  of  the  minor  half  axis  and  EAZ  is  the  angle  from  true 
north  to  the  path  of  the  major  half  axis.  The  shore  site  is 
located  at  point  XI,  Y1  within  the  storm  ellipse.  To  determine 
the  pressure  gradient,  a  second  ellipse  is  plotted  with  axes 
A1  and  B1  which  passes  through  point  XI,  Yl.  The  second  el¬ 
lipse  has  the  same  axial  ratio  B/A  and  the  same  origin  as  the 
storm  ellipse  (Figure  3).  The  tangent  to  an  ellipse  is  defined 
by  equation  1. 

Xj  X  Yx  Y 

- +  -  =  i 

a2  b2  (1) 

therefore,  the  intersection  of  the  tangent  to  the  ellipse  with 
the  X  axis  can  be  found  by  equation  2. 

a2  a2  Y,  Y 

X  =  —  -  -  (2) 

Xj  b2  Xi 

Equation  3  is  used  to  generate  the  line  normal  to  the  tangent 

b2Xj 

X  =  X0  +  -  Y.  (3) 

a2Yj 

The  intersection  of  the  line  normal  to  the  tangent  and  the  A  axis 
is  found  by  equation  4. 

b2  Xj 

x0  =  x - Y  (4) 

a2  Y  i 

which  is  terms  of  the  point  XI,  Yl  would  be 

b- 

Xg  *  X, - Y,  (5) 


The  distance  X,  is  measured  from  the  center  of  the  ellipse  to  the 
point  where  the  line  normal  to  the  tangent  through  XI,  Yl  intersects 
the  A  axis.  The  point  X2,  Y2  is  the  intersection  of  the  line  normal 
to  the  tangent  at  XI  ,Y1  and  the  path  of  the  ellipse. 

To  determine  the  gradient  of  the  barometric  pressure,  it  is 
assumed  that  the  pressure  gradient  follows  a  normal  curve  along 
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Figure  3.  Location  of  the  shore  site  (Xl.Yl)  within  a  storm  ellipse 
(AB)  and  on  a  minor  ellipse  (A1,B1). 
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the  major  axis  of  the  ellipse.  PI  is  the  barometric  pressure 
at  point  X0  along  the  A  axis.  The  pressure  at  X0  is  used  in 
determining  the  pressure  gradient  normal  to  the  isobar  at  point 
XI, Yl.  The  final  pressure  gradient  calculation  Is  made  in  sub¬ 
routine  WIND.  XA  and  YA  are  used  to  plot  the  tangent  to  the 
ellipse  at  XI  ,Y1  for  determining  the  wind  direction.  The  wind 
direction  is  assumed  to  be  parallel  to  the  tangent  to  the  el¬ 
lipse  at  X1,Y1  and  in  a  counterclockwise  direction  around  the 
center  of  the  ellipse. 
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Subroutine  WIND 


The  geostrophic  wind  speed  and  direction  for  each  storm 
position  are  computed  in  subroutine  WIND.  The  equation  for  geo¬ 
strophic  wind  speed  V  is  based  on  latitude  and  barometric  pres¬ 
sure  gradient.  9 


Vg  =  FlTsinT  rfr  <6> 

where  S  is  the  specific  volume,  {779  cm3/gm),  s';  is  the  angular 
velocity,  (7.29  x  10~5  rad. /sec.),  <{>  is  the  latitude  in  degrees, 
and  aP/aN  is  the  barometric  pressure  gradient  normal  to  the  iso¬ 
bar  at  the  shore  location  (Godske  et  al_,  1957,  p.  370).  The  baro¬ 
metric  pressure  gradient  is  computed  at  right  angles  to  the  tan¬ 
gent  of  the  ellipse  through  the  shore  site,  (point  XI,  Yl).  To 
compute  the  gradient,  a  normal  curve  is  constructed  perpendicular 
to  the  tangent  through  point  XI,  Yl .  The  derivative  of  the  normal 
curve  is  taken  at  point  XI,  Yl  to  compute  the  barometric  pressure 
gradient.  The  geostrophic  wind  direction  is  assumed  to  be  parallel 
to  the  tangent  of  the  ellipse  at  point  XI,  Yl  and  heading  in  a  coun¬ 
ter-clockwise  direction  around  the  center  of  the  ellipse.  It  is 
assumed  that  the  small  ellipses  within  the  storm  ellipse  are  par¬ 
allel  to  isobars.  Therefore,  wind  direction  can  be  determined  if 
the  geostrophic  wind  is  directed  along  the  isobars  with  the  high 
pressure  to  the  right  and  low  pressure  to  the  left  of  the  motion 
in  the  northern  hemisphere.  By  means  of  geostrophic  wind  equations, 
the  wind  direction  can  be  estimated  with  error  of  less  than  10°, 
and  speed  with  an  error  of  less  than  20%  (Cole,  1970,  p.  185). 

An  approxiamte  relationship  exists  between  the  speed  and  direction 
of  the  surface  wind  measured  at  anemometer  level  and  the  upper  quasi- 
geostrophic  wind.  Owing  to  differences  in  surface  roughness,  this 
relationship  varies  from  one  station  to  another,  and  also  varies  at 
a  single  station  with  stability.  Thus  it  is  rather  difficult  to  de¬ 
termine  the  surface  wind  accurately  from  the  upper  quasi-geostrophic 
wind.  Under  average  conditions,  a  rough  method  may  be  applied  which 
makes  use  of  the  horizontal  friction  force,  SR,  near  the  ground 
(Figure  4).  The  equation  for  horizontal  motion  can  be  used  for  com¬ 
puting  the  friction  force  at  a  given  station  (Godske,  1957,  p.  453). 


sR  =  v^  +  svhP+  2.<z  k  x  v^ 


(7) 
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The  horizontal  component  of  the  wind  speed,  v^  and  the  barometric 
pressure  gradient,  sv^p,  can  be  measured  directly  and  inserted  in¬ 
to  the  equation.  Although  it  would  be  possible  to  determine  v^,  the 

vertical  variation  in  wind  speed,  it  would  involve  measurements  in 
time  and  space  which  would  be  vary  laborious  and  not  to  reliable. 

In  practice,  v^  includes  both  a  convective  term  and  a  local  term. 

However,  the  influence  of  vh  is  small  when  mean  values  over  rela¬ 
tively  large  areas  are  included  (Hasselberg  and  Sverdrup,  1915). 
Computations  of  the  average  of  R  based  on  time-series  averages  of 
a  series  of  synoptic  series  of  maps  was  carried  out  over  land  (Baur 
and  Phillips,  1938)  and  over  sea  (Westwater,  1943).  Based  on  their 
computations,  the  frictional  force  SR  is  directed  backward  to  the 
right  of  the  wind  vh,  and  is  proportional  to  the  wind  velocity, 

sR  =  bvh  (8) 

and  forms  an  angle  B  with  -vh  as  shown  in  Figure  4.  To  make  cor¬ 
rections  for  wind  speed  and  direction  over  land  and  sea,  the  mean 
values  of  b  and  e  are:  b  =  1.9  x  lO’^sec"1,  3  =  29°  over  land  and 
b  =  0.65  x  10’4sec’’1,  e  =  50°  over  sea  (Baur  and  Phillips,  1938). 

Using  b  and  b,  a  simple  diagram  has  been  constructed  to  show 
the  configuration  of  forces  when  vh  =  0.  The  balance  of  forces  in 

the  diagram  along  v^  and  normal  to  it  gives  the  following  equalities 

with  the  angle  between  the  geostrophic  wind  and  the  surface  wind 
denoted  by  a. 


sin  a  =  sR  cos  e 


bv^  cos 


cos  a  =  sR  sin  i  +  2s;zvh 


therefore,  by  division 

cot  a  =  tans  + 


bvh  sin  b  +  2“zvh 


2  n  sin  $ 
b  cos  B 


where  i  is  the  latitude.  Therefore,  the  angle  u  between  geostrophic 
wind  and  the  surface  wind  can  be  determined  directly  from  the  latitude 
't  when  b  and  A  are  known  for  a  given  station. 

The  ratio  between  the  surface  wind  v^  and  the  geostrophic  wind 
Vg  can  be  derived  from  equations  9  and  10  and  are  given  below  according 


to  Godske  et  al ,  1957. 


1 


therefore. 


_ 

T 

2p  sin  4> 


s  1  vp  |  sin  a  2 u  sin  ^ 
b  cos  6  s  |  vp  j 


vh  _  2  P  sin  $  sin  a 

b  cos  3  (11) 


Once  the  angle  a  between  the  surface  wind  and  geostrophic  wind 
has  been  computed,  it  can  be  inserted  into  equation  11  to  compute 
the  ratio  between  the  surface  wind  and  the  geostrophic  wind  over 
land  or  sea.  Table  1  gives  values  at  different  latitudes  for  u, 
the  angle  between  surface  wind  and  geostrophic  wind,  and  vh/vg, 

which  are  plotted  in  Figure  5.  In  subroutine  WIND,  the  correction 
factors  for  computing  surface  wind  speed  and  direction  are  com¬ 
puted  following  statement  50. 

Since  the  values  for  b  and  3  are  given  for  wind  blowing  over 
land  or  over  sea,  intermediate  values  must  be  computed  for  winds 
blowing  along  the  shore.  Winds  blowing  directly  onshore  with  a 
wind  angle  of  zero  would  have  values  of  b  =  .000065  and  6  =  50. 

In  this  case,  the  wind  is  blowing  from  the  sea  and  the  land  does 
not  have  any  frictional  effect  on  the  wind.  In  like  manner,  if 
the  wind  is  blowing  over  the  land  in  an  offshore  direction,  values 
for  the  land,  b  =.000190  and  p  =  20  are  applied.  Fo,r  the  transi¬ 
tion  zones,  a  cosine  transformation  is  used  to  compute  the  inter¬ 
mediate  values.  Within  the  transition  zone,  angle  A  is  computed 
from  0  to  90°  with  0  being  land  and  90°  being  sea.  The  new  angle 
A  is  used  in  equations  12  and  13  to  compute  the  transition  values 
for  b  and  beta. 


b  =.0001  •  (1.275  +  .625  (sin  A))  (12) 

3  =  39.5  -  10.5  •  sin  A  (13) 

The  computed  values  for  b  and  3  are  substituted  in  equations  10  and 
11  to  compute  the  surface  wind  speed  and  direction  from  the  geo¬ 
strophic  wind. 

The  final  step  in  subroutine  WIND  is  to  compute  the  effective 
wind  speed  which  is  carried  over  into  subroutines  FETCH  and  WAVES 
for  determining  effective  fetch  length  and  wave  height.  The  effec- 
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tive  wind  speed  is  that  which  will  generate  waves  which  will  in  turn 
have  an  effect  on  the  beach.  If  the  wind  is  blowing  directly  on¬ 
shore,  the  full  force  of  the  wind  is  used  in  generating  waves  which 

will  hit  the  coast.  However,  if  the  wind  is  blowing  directly  off¬ 

shore,  small  waves  will  be  generated  in  the  nearshore  area  (Resio 
and  Hayden,  1973).  Based  on  empirical  observations,  onshore  winds 
are  about  three  times  as  effective  in  generating  local  waves  as 
offshore  wind  (Davis  and  Fox,  1974  and  Owens,  1975).  Therefore, 
a  cosine  transformation  is  used  in  subroutine  WIND  to  compute  ef¬ 
fective  wind  speed  from  the  surface  wind  speed  and  direction. 

When  the  wind  is  blowing  directly  onshore,  the  effective  wind  is 
equal  to  the  surface  wind.  On  the  other  hand,  when  the  wind  is 
blowing  along  the  shore,  the  effective  wind  is  equal  to  2/3  of 
the  surface  wind,  and  when  the  surface  wind  is  blowing  directly 
offshore,  the  effective  wind  is  1/3  the  surface  wind.  Although 
the  values  for  effective  wind  may  be  rough  in  some  cases,  they 
seem  to  give  good  estimates  where  comparative  wind  and  wave  data 
are  available. 
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Latitude 


\  _^h 

vg  L  vg  s 


0 

61 

40 

— 

— 

10 

55 

29 

— 

— - 

20 

49.5 

23 

— 

— 

30 

45 

19 

0.31 

0.56 

40 

42 

16 

0.38 

0.63 

50 

39 

14.5 

0.42 

0.67 

60 

37 

13 

0.46 

0.70 

70 

36 

12.5 

0.485 

0.715 

80 

35 

12 

0.495 

0.723 

90 

35 

12 

____ 

— 

Table  1.  The  angle  between  surface  wind  and  geostrophic 
over  land  ct^  and  over  sea  ,  and  the  ratio 

between  surface  wind  speed  and  geostrophic 
wind  speed  over  land  and  sea  for  different 
latitudes  <(>,  according  to  Baur  and  Phillip 
(1938,  p.  292). 
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/GEOSTROPHIC  WIND) 


0  10  20  30  40  50  60  70  80  90 

LATITUDE  (°) 


LATITUDE  (•) 

Figure  5.  Angles  between  surface  and  geostrophic  winds  and  ratio  of 
surface  to  geostrophic  wind  speed. 


Subroutine  DECAY 


Subroutine  DECAY  is  used  to  determine  the  height  of  a  given 
wave  after  it  has  decayed  for  a  specified  length  of  time.  Snod¬ 
grass,  et  al  (1966)  presented  empirical  data  on  the  attenuation 
of  selected  frequencies  which  they  observed  in  their  study  of 
propagation  of  ocean  swell  in  the  Pacific.  In  general,  they 
found  the  attenuation  to  be  large  within  the  limits  of  the  wind 
area  of  the  storm,  and  small  outside  the  storm  area.  The  em¬ 
pirical  attenuation  data  were  logarithmic  coefficients  reported 
in  units  of  decibels  per  latitude  degree  of  propagation  distance 
For  the  range  of  frequencies  0.06  to  0.08  Hertz,  these  data  fit 
an  attenuation  function  of  the  form 


e 


-2ax 


(14) 


where  a  is  the  modulus  of  amplitude  decay  in  degree'  =0.1151  e, 
where  a  is  the  logarithmic  attenuation  coefficient  in  decibels/ 
degree,  and  where  x  is  the  propagation  distance  in  degrees. 

The  logarithmic  attenuation  coefficient  versus  frequency  was 
plotted  on  a  graph  (Kaufman,  1973),  and  equation  (15)  was  de¬ 
rived  from  the  line  on  the  graph. 


a  =  10 


(15) 


where  F  is  the  frequency  of  the  wave  being  decayed. 

The  propagation  distance  x  is  found  by  multiplying  the  wave 
velocity  1.5606*T  (where  T  is  the  wave  period)  times  the  time 
interval  TINT.  This  distance  is  then  reduced  to  degrees  by 
multiplying  it  by  the  constant  value 


360°/circle 

40074km/circle 

(using  the  circumference  of  the  earth  at  the  equator).  To  find 
the  decayed  wave  height  then,  the  original  wave  height  is  multi¬ 
plied  by  the  attenuation  function.  This  decay  factor  was  tested 
with  several  different  wave  heights,  periods,  and  time  intervals, 
giving  very  reasonable  decay  results,  but  there  were  no  empirical 
data  against  which  to  check  the  results. 


Subroutine  ETIME 


Subroutine  ETIME  is  used  to  determine  the  amount  of  time 
(referred  to  as  effective  duration)  which  would  be  required  to 
produce  waves  of  a  certain  height  by  wind  blowing  at  a  given  wind 
speed.  A  wave  forecasting  procedure  developed  by  Sverdrup  and 
Munk  (1947),  and  revised  by  Bretschneider  (1952,  1958)  with  ad¬ 
ditional  empirical  data  is  called  the  sverdrup-Munk-Bretschenider 
(SMB)  method  (C.E.R.C.,  1973).  The  SMB  curves  for  forecasting 
wave  height  are  based  on  equation  16  from  Bretschneider  (1958). 


jjjr  =  0.283  Tanh  [_  0.0125  (  ]  (16) 

where  g  is  the  acceleration  due  to  gravity,  H  is  the  wave  height, 

F  is  the  effective  fetch  length,  and  U  is  the  wind  speed.  Solving 
equation  16  for  F  gives 

E  rARCTANH(  0.a3ag  u.  n 
9  L  0.0125  J  (17) 

Therefore,  F  is  the  effective  fetch  that  it  would  take  to  generate 
waves  of  height  H  with  a  wind  speed  of  U. 

In  terms  of  storm  duration,  the  effective  fetch  equation  is 


F 


S 

10 


0.72 


10°'3x  D1*25 


(18) 


Where  F  is  effective  fetch,  W  is  wind  speed,  and  D  is  storm  duration. 
Solving  this  for  0,  we  have: 


D  = 


(™-  100.3 


(19) 


Therefore,  D  is  the  effective  duration  that  it  would  take  to  build 
waves  of  a  certain  height  (used  to  find  the  effective  fetch)  with 
winds  of  a  given  speed.  Then,  this  effective  duration  is 
added  to  the  current  time  increment  of  the  storm  to  give  a  duration 
which  takes  into  account  the  wave  built  in  previous  time  increments. 

Subroutine  ETIME  was  tested  by  running  the  data  from  subroutine 
WAVES  back  through  it  to  arrive  at  the  original  data.  For  example, 


the  following  wave  heights  and  wind  speeds  were  tested  using 
the  C.E.R.C.  charts  and  the  subroutine  (Table  2). 


Wave  Height 

Wind  Speed 

Effective  Duration 

From  SMB  Curves  Derived 

14  ft 

80  kts 

1.1  hrs 

1.05  hrs 

1  ft 

12  kts 

1 .4  hrs 

1.38  hrs 

45  ft 

80  kts 

9.9  hrs 

10.23  hrs 

3  ft 

12  kts 

21.0  hrs 

21.40  hrs 

14  ft 

30  kts 

21.0  hrs 

21.76  hrs 

Table  2.  Test  for  subroutine  ETIME  comparing  effective 
duration  available  from  SMB  curves  with  derived 
duration  for  selected  wave  heights  and  wind 
speeds. 
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Subroutine  FETCH 


Subroutine  FETCH  is  used  to  determine  the  maximum  wind  gen¬ 
erating  area  and  the  average  wind  speed  within  that  area  for  a 
circular  storm  on  the  open  ocean.  Three  cases  are  considered  for 
determining  the  actual  storm  fetch  as  the  storm  passes  across  a 
shoreline. 

For  a  storm  on  the  open  ocean,  the  wind  speed  varies  with 
the  slope  (first  derivative)  of  a  normal  curve  as  one  crosses 
the  storm.  The  area  of  the  storm  with  wind  speed  greater  than 
one  half  of  the  maximum  lies  within  a  ring  with  an  outside  radius 
of  1.92  standard  deviations  and  an  inner  radius  of  0.32  standard 
deviations.  This  ring  between  1.90  and  0.30  was  divided  into  13 
smaller  rings.  The  area  of  each  ring  was  found  and  multiplied 
by  the  wind  speed  at  the  midpoint  of  the  ring.  These  products  were 
summed  and  then  divided  by  the  total  area  to  give  an  average  wind 
speed  of  0.8847  of  the  maximum  wind  speeu.  The  total  area  was 
found  to  be  7.7598  square  standard  deviations.  Since  the  wind  is 
being  generated  in  a  circular  pattern,  approximately  one  quarter 
of  the  wind  is  blowing  along  each  axis  of  a  grid  with  its  center 
at  the  origin.  Since  the  storm  fetch  area  has  a  shape  resembling 
an  ellipse,  the  fetch  for  winds  blowing  in  any  one  direction  can 
be  considered  an  ellipse  with  an  area  equal  to  one  quarter  of  the 
maximum  wind  generating  area  (1-9400  square  standard  deviations) 
and  centered  on  the  maximum  wind  speed  circle  (1  standard  devi¬ 
ations).  The  short  half  diameter  was  taken  to  be  0.7  standard 
deviations  (1.0-0. 3). 

Area  =  i  x  a  x  b 
1.9400  =  3.1415  x  a  x  .7 
a  =  0.8822  standard  deviations  (20) 

Therefore,  the  maximum  storm  fetch  length  Is  twice  that,  or  1.7644 
standard  deviations.  Converting  to  storm  radii,  the  fetch  is  0.5881 
radi i . 

Where  the  storm  crosses  the  shoreline  three  cases  must  be  con¬ 
sidered;  Case  1,  where  1/3  R  -*  x  >_  0;  Case  2,  where  0.4444  R  >  x  >_  1/3 
R  and  Case  3,  where  x  >  0.4444  R,  where  R  is  the  storm  radius  and  x 
is  the  distance  from  the  center  of  the  storm. 

For  Case  1,  (1/3  R  >  x  ^  0)  (Figure  6),  it  is  assumed  that  the 
fetch  area,  centered  on  the  1/3  R  circle,  swings  around  so  that  the 
long  axis  of  the  ellipse  is  always  pointing  at  the  beach.  We  have 
maximum  storm  fetch  in  this  case  until  Y  is  small  enough  that  the 
fetch  begins  to  cross  the  beach  as  it  continues  to  swing.  This 
point  is  Y1 .  From  the  Case  One  diagram,  D?  =  X;  +  Yl:  and  D; = 
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0/3  R)2  +  (1/2  F)2  so  we  have  X2  +  Yl2  =  {1/3  R)2  +  (1/2  F)2  or 


Y1  =  /(1/3  R)?  +  (1/2  F)2  -  X2  (21) 

As  the  storm  continues  to  cross  the  beach  and  rotate,  the  fetch 
decreases.  When  the  1/3  R  circle  of  the  storm  crosses  the  beach, 
the  storm  has  swung  so  that  the  beach  is  in  the  center  of  the  fetch 
area.  Therefore,  the  storm  fetch  is  one  half  of  maximum  at  this 
point  (Y2).  As  the  storm  continues  inland,  the  fetch  length  reaches 
a  minimum  as  Y  goes  to  zero  (Y3).  This  minimum  is  determined  by 
the  equation  22. 


1X1 

FMIN  =  •  FMAX.  (22) 

This  gives  a  diminishing  minimum  as  f x !  goes  to  zero.  The  fetch 
then  increases  to  a  high  equal  to  1 / 2  *  FMAX  as  the  1/3  R  circle  a- 
gain  crosses  the  beach  (Y4).  The  fetch  then  declines  again,  fin¬ 
ally  dropping  to  zero  at  Y4,  at  which  point  the  storm  is  completely 
onshore.  This  gives  the  "eye  of  the  storm"  effect  where  the  winds 
drop  to  a  minimum  when  the  storm  is  at  its  closet  point  of  approach 
and  then  increase  again  as  the  storm  moves  on,  and  finally  drop  as 
the  storm  moves  away.  Cosine  smoothing  is  used  to  smooth  out  the 
increases  and  decreases  in  fetch. 

For  Case  2  (0.4444  R  ■  x  1/3  R)(Figure  7),  again  we  have  full 

fetch  until  the  fetch  area  comes  into  contact  with  the  beach  at  Y1 . 

The  fetch  then  decreases  as  the  storm  moves  onshore,  hitting  zero 
at  Y2.  The  reason  that  the  fetch  doesn't  drop  and  then  rise  again, 
as  in  Case  One,  is  that  the  eye  of  the  storm  (the  part  inside  the 
1/3  R  circle)  never  passes  over  the  beach. 

This  case  ends  at  X  =  0.4444  R  because  D  in  the  case  two  dia¬ 
gram  =  ✓'(1/2  R):’+  (1/2  F)2  =  0.4444  R.  When  X  _  D,  then  the  geo¬ 

metric  relationships  that  allow  us  to  find  Y1  no  longer  hold,  since 
0,  the  distance  from  the  origin  to  the  beach,  can  no  longer  equal 
both  /X2  +  Yl2  and  »'(l/3  R)2  +  (1/2  F)2.  Again  consine  curves 
are  used  to  smooth  out  the  changes  in  storm  fetch. 

For  Case  3,  (  x  ■  0.4444  R)(Figure  8)  the  storm  fetch  area  never 
actually  crosses  the  beach,  but  the  fetch  decreases  as  the  storm  pas¬ 
ses  over  the  shore  farther  along  the  coast.  Thus  the  storm  fetch  is 
at  its  maximum  until  the  fetch  area  starts  going  onshore  at  Yl .  Here 
the  fetch  is  perpendicular  to  the  storm  path.  Yl  is  at  zero  since 
at  that  point,  the  whole  half  of  the  storm  that  is  generating  along¬ 
shore  waves  is  still  over  the  water,  so  we  still  have  full  fetch. 

As  the  entire  storm  passes  onshore  (Y2),  the  fetch  drops  to  zero. 

Y2  is  -vW  -  (0.4444  ITp  .  Y2  is  frozen,  using  an  X  value  of  0.4444R. 
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Case  2  -  storm  fetch  when  the  distance  from  center  of  storm 
X  is  less  than  0.4444  and  qreater  than  1/3  storm  radius,  R. 


Figure  8.  Case  3  -  storm  fetch  when  the  distance  from  center  of  storm 
X  is  greater  than  0.4444  storm  radius,  R. 

■n 


The  reason  that  this  is  valid,  is  that  from  X  =  0.4444  R  on  out, 
the  fetch  area  never  actually  passes  over  the  beach,  so  that  all 
these  cases  are  essentially  the  same,  as  far  as  the  size  of  their 
storm  fetches  goes.  Only  the  distance  from  the  storm  fetch  to 
the  beach  changes. 

Subroutine  FETCH  was  tested  in  two  ways.  First  X  and  Y  were 
varied  independently,  running  Y  from  +R  to  -R  for  each  value  of  X. 
This  simulated  storms  with  paths  perpendicular  to  the  shoreline. 
Second,  X  and  Y  were  varied  simultaneously  by  a  constant  amount, 
simulating  storms  crossing  the  shoreline  at  an  angle.  In  both 
cases,  the  program  produced  continuous  curves  of  fetch  varying  with 
storm  location,  with  smooth  transitions  between  all  cases. 


Subroutine  WAVES 


The  equations  for  predicting  significant  wave  height  and 
period  in  subroutine  WAVES  are  based  on  the  Sverdrup-Munk- 
Bretschneider  (SMB)  method  revised  by  Bretschneider  (1958) 
and  plotted  as  a  series  of  curves  by  C.E.R.C.  (1973).  The  SMB 
wave  forecasting  curves  for  fetches  of  1  to  1000  miles  are  given 
in  Figure  9.  The  wave  prediction  curves  use  the  wind  speed  in 
knots,  storm  duration  in  hours,  and  storm  fetch  to  calculate  the 
significant  wave  height  and  significant  wave  period. 

In  order  to  use  the  SMB  method  in  the  model,  the  first  task 
was  to  find  equations  to  approximate  the  effective  fetch  from  the 
wave  prediction  curves.  The  effective  fetch  is  the  limiting  fetch 
which  corresponds  to  a  given  wind  speed  and  duration.  The  effect¬ 
ive  fetch  is  determined  by  moving  to  the  left  across  the  chart  at 
the  level  of  the  wind  speed  until  you  hit  the  appropriated  storm 
duration  line.  Then  drop  straight  down  to  the  fetch  length  axis 
from  the  intersection  of  the  wind  speed  line  with  the  duration 
line.  This  value  on  the  fetch  length  axis  is  the  effective  fetch, 
if  it  is  less  than  the  actual  fetch. 

To  develop  an  equation  for  effective  fetch,  the  first  problem 
is  to  determine  the  intercept  of  the  proper  duration  line.  To  do 
this,  the  intercept  of  each  line  with  the  fetch  length  axis  was 
plotted  against  the  storm  duration.  Log  scales  were  used  on  both 
axes  since  the  original  fetch  length  axis  had  a  log  scale  and  the 
duration  lines  themselves  were  spaced  logarithmically.  These 
formed  a  nearly  linear  trend  and  the  equation  for  the  line  was 
found,  in  terms  of  the  log  axes,  to  be: 

I  =  100,3  X  D1 '25  (23) 

where  I  is  the  intercept  and  D  is  the  storm  duration. 

The  next  step  was  to  determine  the  slope  of  the  storm  duration 
lines  given  in  the  following  equation. 

/  s  .0.72 

f  =(-fo)  x  I  (24) 

where  F  is  the  effective  fetch,  S  is  the  wind  speed,  and  I  is  the 
intercept  of  the  duration  line  with  the  fetch  length  line.  Com¬ 
bining  equations  23  and  24  we  get  equation  25  for  the  effective 
fetch  in  terms  of  wind  speed  and  storm  duration: 

F-(V'72  *  O0°'3  x  o’ ,25)  (25) 
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If  this  value  is  less  than  the  actual  storm  fetch  then  it  decreases 
the  relevant  fetch  length  which  is  used  in  the  rest  of  the  equations. 


The  SMB  forecasting  curves  were  constructed  from  equations  26 
and  27,  which  were  empirically  derived  by  Bretschneider  (1958). 


jjr  =  0.283  TANH  (  0.0125(^)°*4" )  (26) 


2§jJ=  T-20  ™NH  (  0.077  (  )°"  )  (27) 

where  g  is  the  gravitational  constant,  p  is  PI  (3.1459),  H  is 
the  significant  wave  height,  U  is  the  wind  speed,  F  is  the  effective 
fetch,  and  T  is  the  significant  wave  period.  Solving  these  equations 
for  H  and  T,  we  get  equations  28  and  29. 


H  = 


U2  x  0.283  x  TANH \  0.125 
9 


(28) 


T  =  2pU  x  1.20  x  TANH  (  0.077  ifo  (29) 

g 


The  values  for  wind  speed,  storm  duration,  and  effective  fetch  are 
then  inserted  into  equations  28  and  29  to  yield  the  wave  height  and 
period.  For  example  in  the  case  of  a  storm  with  a  duration  of  10 
hours,  a  wind  speed  of  35  knots,  and  an  actual  fetch  of  200  nautical 
miles,  this  gives  us  an  effective  fetch  of  87.44  nautical  miles,  a 
wave  height  of  12.78  feet,  and  a  wave  period  of  7.85  seconds.  But 
suppose  that  we  have  a  storm  the  same  as  the  last  one,  but  with  a 
fetch  of  only  80  nautical  miles.  In  this  case  the  actual  fetch  is 
smaller  than  the  computed  fetch,  so  it  remains  as  the  effective 
fetch.  This  gives  us  a  wave  height  of  12.36  feet  and  a  wave  period 
of  7.72  seconds. 

The  program  has  an  option  so  that  the  results  can  either  be 
metric  or,  to  facilitate  checking  the  results  against  the  SMB  curves, 
the  results  can  be  in  nautical  miles  and  feet.  The  subprogram  was 
tested  with  numerous  combinations  of  wind  speeds,  durations,  and 
fetch  lengths,  with  the  results  agreeing  very  well  with  the  SMB 
forecasting  curves. 
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F«'C*  (  Stetut*  Mm**) 


Subroutine  TIDE 


Subroutine  TIDE  is  used  to  determine  the  tide  level  at  each 
hour,  the  spring  tide  range  (ST),  the  neap  tide  range  (NT),  the 
number  of  days  since  the  last  spring  high  tide  (TDAY),  the  hour 
of  the  last  high  spring  tide  (THR),  and  the  tidal  form  number 
(FN).  Four  principal  tidal  components,  -  Principal  lunar, 

S2  -  Principal  solar,  K-j  -  Lunar-solar  diurnal,  and  0-j  -  Prin¬ 
cipal  lunar  diurnal  are  used  for  making  a  prediction  of  the  hourly 
tide  level.  The  periods  of  the  semi-diurnal  components  (TM2  = 
12.42  hours  and  TS2  =  12.00  hours)  and  the  diurnal  components 
(TK1  =  23.93  hours  and  T01  =  25.82  hours)  are  constants  in  the 
subroutine.  The  tidal  form  number  FN  is  used  to  classify  the 
tides  of  a  locality  according  to  equation  30  (Defant,  1960,  p. 
306). 


FN  = 


K1  +01 

m2  +  s2 


(30) 


The  following  classification  based  on  Dietrich  (1944,  p.  69)  is 
used  to  classify  tides  according  to  their  form  number. 


FN  =  0  -  0.25  Semi-diurnal  tide 

FN  =  0.25  -  1.50  Mixed-  mainly  semi-diurnal  tide 

FN  =  1.50  -  3.00  Mixed-  mainly  diurnal  tide 

FN  =  greater  than  3.0  Diurnal  tide 

As  examples  of  the  different  types  of  tides,  Immingham,  England 
has  a  semi-diurnal  tide  with  a  form  number  of  0.11.  San  Francisco, 
California  has  a  mixed,  dominately  semi-diurnal  tide  with  a  form 
number  of  0.90.  Manila  has  a  mixed,  dominately  diurnal  tide  with  a 
form  number  of  2.15,  and  Do  San, Viet  Nam  on  the  Gulf  of  Tonkin  has 
a  pronounced  diurnal  tide  with  a  form  number  of  19.2.  The  four  ma¬ 
jor  components  are  responsible  for  the  general  form  of  the  tides 
and  generallyaccount  for  about  70  percent  of  the  total  variance. 

If  the  next  three  most  important  tidal  components,  N^,  K2  and  P-j 

are  included,  the  percentage  of  the  total  variance  increases  to 
about  83%  (Defant,  1960). 


If  some  simplifying  assumptions  are  made  concerning  the  major 
tidal  components,  it  is  possible  to  make  a  good  approximation  of 
the  hourly  tide  level  from  the  maximum  spring  tide  range,  minimum 
neap  tide  range,  and  the  tidal  form  number.  First,  it  is  necessary 
to  assume  that  the  diurnal  components,  K-|  and  0-j  are  approximately 

equal.  Second,  assume  that  the  maximum  spring  tide  range  is  equal 
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to  the  sum  of  the  four  major  components  according  to  equation  31. 


ST  =  M2  +  S2  -  (^  +  Oj)  (31) 

Next  it  is  assumed  that  the  neap  tide  range  is  approximately  equal 
to  the  lunar  components  minus  the  solar  components. 

TN  *  Mg  +  K-|  -  (S2  +  0-,)  (32) 


If  K-j  and  0-j  are  approximately  equal,  it  follows  that 

TN  =  M2  -  S2  (33) 


The  form  number  is  the  diurnal  components  over  the  semidiurnal 
components  in  equation  30,  however,  if  K-|  =  0-| ,  then 


FN  = 


2  Kj 

M2+  S2 


(34) 


By  combining  the  equations  for  the. form  number  (Equation  30), 
the  spring  tide  range, and  the  neap  tide  range.it  is  possible  to 
solve  for  M-j  and  S2< 


M 


ST  +  TN 
1  '  4  (1  +  FN) 


(35) 


and 


_  ST  -  TN 
2  4  (1 -FN) 


(36) 


If  it  is  now  assumed  the  lunar  components  are  proportional, 
an  approximation  of  the  K-|  component  can  be  derived  from  the  M2 

component 


K1  =  FN  •  M2  (37) 

and  the  solar  components  are  related  in  like  manner,  therefore, 

01  =  FN  •  S2  (38) 

The  amplitude  of  the  maximum  spring  tide  was  taken  at  the  last 
previous  spring  tide  for  each  run,  therefore,  the  phases  for  the  four 
major  components  are  considered  0  at  that  time.  By  computing  the 
time  differences  from  the  last  spring  high  tide  to  the  hour  for  the 
prediction,  the  contribution  for  each  tidal  component  can  be  calcu¬ 
lated.  The  argument  (ARG)  is  equal  to  2  pi  times  the  number  of  hours 
since  the  last  spring  tide.  The  tide  is  computed  by  adding  together 
the  contribution  for  each  of  the  tidal  components. 
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Although  some  rough  approximations  were  made  in  deriving  the  ji 

major  tidal  components  from  the  spring  tide  range,  neap  tide  range  1 

and  form  number,  the  resulting  tide  predictions  work  out  quite  i 

closely  with  the  tide  tables.  The  tide  tables  give  the  time  of  l 

high  and  low  tides  for  each  day,  and  the  predicted  times  of  high  j 

and  low  tides  fall  within  1  hour  using  subroutine  TIDE.  The  sub-  ! 

routine  was  tested  for  Plum  Island,  Massachusetts,  Cedar  Island,  i 

Virginia,  Sapelo  Island,  Georgia  and  Mustang  Island,  Texas,  and  H, 

gave  satisfactory  predictions  for  each  of  the  areas.  , 

One  of  the  major  reasons  for  making  tidal  predictions  in  the 
model  is  to  determine  the  effect  that  tides  have  on  the  nearshore 

bottom  slope,  The  slope  at  low  tide  SLBOT,  the  slope  at  high  tide  H 

SLTOP  and  the  tide  level  are  used  to  determine  the  intermediate 
slope  between  high  and  low  tide. 

SLOPE  =  SLBOT  +  TLOC  *(SL0TP  -  SLBOT)  (40)  j 

The  final  tide  level  which  is  included  with  the  output  is 
computed  by  adding  the  relative  tide  level  TIDX  to  the  mean  low  j‘ 

tide  level  TMEAN.  j 


Subroutine  SURF 


Breaker  height,  angle  and  longshore  current  velocity  are 
computed  in  Subroutine  SURF.  The  critical  value,  H^/h^  =  .78 

where  Hb  and  hb  are  breaker  height  and  depth,  are  used  for  a 

breaking  criterion  (Munk,  1949).  Applying  any  wave  theory  and 
assuming  conservation  of  energy  flux,  Komar  and  Gaughan  (1972) 
derived  the  relationship 


Hb  =  0.73  cm  +  .383  g1/5  (T  H2)2/5  (41) 

where  Hb  is  the  breaker  height,  g  is  gravity,  T  is  wave  period 
and  is  the  deep  water  wave  height. 

The  breaker  angle  ab  is  computed  by  first  finding  the  shallow 

water  wave  length  and  then  taking  the  ratio  of  shallow  water  to 
deep  water  wave  length  using  Snell's  law  to  determine  the  breaker 
angle. 


sin  «b 


sin 


(42) 


where  ac  is  the  deep  water  wave  angle  which  is  assumed  to  be  the 
same  as  the  wind  angle,  Hb  is  the  breaker  height  and  l_b  is  the 
breaker  depth. 

The  refracted  breaker  height,  HR,  is  obtained  from  the  re¬ 
fraction  coefficient,  KR, 

cos  (a0) 

KR  =  - 7 — r — 

C0S  (43) 

which  is  multiplied  times  the  breaker  height,  Hb 


Hr  =  KR  •  Hb 


(44) 


Four  different  options  are  available  for  computing  the  long¬ 
shore  current  velocity.  The  longshore  current  equations  by  Longuet- 
Higgins  (1970),  Komar  and  Inman  (1970),  Fox  and  Davis  (1972)  and 
Coastal  Engineering  Research  Center  (C.E.R.C.,  1973)  used  basic 
different  assumptions  with  the  same  set  of  variables.  The  variation 


in  longshore  current  velocity  across  the  surf  zone  and  along  the 
shore,  as  well  as  differences  in  nearshore  topography  brought  about 
by  bars  and  rip  channels  make  any  prediction  of  average  longshore 
current  velocity  very  difficult.  However,  in  making  predictions 
about  the  surf  zone,  it  is  essential  to  at  least  have  a  good  esti¬ 
mate  of  the  maximum  longshore  current  velocity. 

The  radiation  stress  theory  of  Longuet-Higgins  (1970)  has 
been  tested  with  laboratory  data  from  Galvin  and  Eagleson  (1965), 
and  field  data  from  Putman,  Munk  and  Traylor  (1949).  The  long¬ 
shore  current  velocity  in  the  surf  zone,  Vb>  is  a  function  of  the 

bottom  slope,  m,  the  breaker  height,  Hb>  and  breaker  angle,  ab> 

between  the  wave  crest  and  the  shoreline  (Longuet-Higgins,  1970). 

Vb  =  Mlm  <9Hb)V2  sin  2ctb  (45) 

where  Mp  the  friction  factor  is: 


The  longshore  current,  Vb  is  measured  at  the  breaking  position  and 

r  is  a  mixing  coefficient  with  a  range  of  0.17  (little  mixing)  to 
0.5  (complete  mixing)  with  a  mode  at  about  0.2.  The  depth  to 
height  ratio  in  shallow  waver,  3,  is  taken  to  be  1.2  and  ff  the 

friction  coefficient  is  set  at  0.01.  By  inserting  the  above  values 
in  equation  46,  the  value  for  becomes  9.0.  Therefore,  the  long¬ 
shore  current  equation  according  to  Longuet-Higgins  (1970)  can  be 
reduced  to: 

Vb  =  9.0  m  (gHb)1/2  sin  2%  (47) 

When  equation  47  was  applied  to  test  sets  of  field  and  lab¬ 
oratory  data  by  C.E.R.C.  (1973),  the  data  yields  predictions  that 
average  about  0.43  of  the  measured  values.  The  measured  values 
were  taken  in  the  fastest  field  of  flow  shoreward  of  the  breaker 
zone,  whereas  the  predictions  were  made  for  longshore  current  at 
the  line  of  breakers.  Therefore,  it  has  been  proposed  by  C.E.R.C. 
(1973)  that  the  Longuet-Higgins  equation  be  multiplied  by  2.3  to 
yeild  the  C.E.R.C.  equation: 

V  =  20.7  m  (gHfa)1/2  sin  2*b  (48) 
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Komar  and  Inman  (1970)  derived  a  longshore  current  equation 
based  on  radiation  stress.  Where  the  radiation  stress  components 
defined  by  Longuet-Higgins  and  Stewart  (1964)  is  the  excess  flow 
of  momentum  due  to  the  presence  of  waves.  The  Komar  and  Inman 
(1970)  equation  is: 

v  3  clum  si"  -b  cos  «b  <49> 

where  V  is  the  longshore  current  velocity.  Tan  $  is  the  beach 
slope,  is  the  bottom  frictional  drag  coefficient.  Um  is  the 

maximum  horizontal  component  of  the  orbital  velocity  of  the  waves 
and  C-|  is  a  dimensional  coefficient  of  proportionality.  However, 

Komar  (1969)  suggested  that: 


(Tan  6  cos  ab)/Cf  =  constant  (50) 

indicating  that  the  variation  in  beach  slope  does  not  produce  a 
change  in  longshore  current  velocity.  Therefore,  the  Komar  and 
Inman  (1970)  longshore  current  equation  becomes: 


V  =  C, U  sin  a. 

I  m  b 


(51) 


A  fourth  equation  developed  by  Fox  and  Davis  (1972)  uses  em¬ 
pirical  data  subjected  to  linear  regression  analysis  to  predict 
longshore  current  velocity.  The  linear  regression  analysis  is 
based  on  3  sets  of  data  collected  at  Stevensvi 1 le ,  Michigan  (Fox 
and  Davis,  1970),  Holland,  Michigan  (Fox  and  Davis,  1971a)  and 
Sheboygan,  Wisconsin  (Fox  and  Davis,  1972).  Each  set  of  data 
consists  of  360  observations  taken  at  2  hour  intervals  for  30  days 
of  longshore  current  speed  and  direction,  breaker  height,  period 
and  breaker  angle.  Using  a  stepwise  regression  analysis,  the  con¬ 
tribution  of  each  variable  was  tested  separately,  and  then  in  var¬ 
ious  combinations.  The  ratio,  H^/T  is  related  to  the  mass  flux  on 

volume  of  water  which  enters  the  surf  zone  and  must  be  removed  by 
the  longshore  current.  The  breaker  angle,  a^,  defines  the  angle 

between  the  breaker  crest  and  the  shoreline  and  is  therefore  re¬ 
lated  to  the  momentum  transfer  in  the  longshore  direction.  Using 
the  regression  program,  a  series  of  combinations  was  tested  for 
the  sin  of  the  angle  including  sin  A,  sin  2A,  sin  3A,  sin  4A... 

sin  8A.  The  closest  fit  was  obtained  when  sin  4A  was  used  for  the 

angles.  For  the  1969  set  of  data  from  Stevensville,  Michigan,  the 
following  equation, 


V  =  5.42^-—^  sin  4A 


(52) 


gave  the  best  fit  and  accounted  for  83.5  percent  of  the  total  sum 
of  squares.  For  the  1970  data  from  Holland,  Michigan,  the  coef¬ 
ficient  of  proportionality  was  3.47  and  the  equation  accounted  for 
78.8  percent  of  the  total  sum  of  squares.  For  the  1972  data  from 
Sheboygan,  Wisconsin,  the  coefficient  was  2.98  and  the  equation 
accounted  for  77.8  percent  of  the  total  sum  of  squares. 

The  three  areas  differed  in  the  nearshore  bottom  slope  and 
the  occurrence  of  sand  bars  which  influenced  the  coefficient  of 
proportionality.  The  coefficient  for  each  case  was  approximately 
equal  to  100  times  the  bottom  slope.  Therefore,  the  longshore 
current  velocity  according  to  Fox  and  Davis  (1973)  is 

V  =  100  mf-y  ^  sin  4  y  (53) 


When  the  four  longshore  current  equations  were  tested  in  the 
model,  the  equation  by  Longuet-Higgins  (1970)  and  Fox  and  Davis 
(1973)  gave  very  similar  results  for  breaker  angles  up  to  about 
20  degrees.  For  higher  breaker  angles,  the  predicted  results  from 
the  Fox  and  Davis  (1973a)  equation  were  too  low.  The  values  for 
longshore  current  predicted  by  Komar  and  Inman  (1970)  and  C.E.R.C. 
(1973)  were  consistently  too  high.  Although  the  four  equations  are 
available  as  options,  it  is  recommended  that  the  Longuet-Higgins 
(1970)  equation  be  used  for  making  predictions.  If  possible,  it 
is  best  to  test  predictions  with  hindcast  data  from  the  same  area. 
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Subroutine  ENRGY 


Subroutine  ENRGY  is  used  to  determine  the  wave  energy  during 
each  hour  of  the  storm  which  is  summed  to  give  the  total  wave 
energy  for  the  storm.  The  deep  water  wave  energy  E  (C.E.R.C., 

1 973 )  is  given  by 


Eo  = 


p  g  H  L0 


-  5.12  ,  g  (HT)- 

8 


(54) 


where  is  the  mass  density  of  the  water  which  is  1.94  slugs/cubic 
foot  for  fresh  water  and  2.0  slugs/cubic  foot  for  salt  water,  H  is 
the  deep  water  wave  height  and  T  is  the  wave  period.  Conversion 
factors  are  included  to  change  from  foot  pounds/foot  to  Joules/ 
meter.  The  subroutine  was  tested  using  wave  energy  calculation 
from  previous  studies  (Fox  and  Davis,  1971b). 


Subroutine  ARCTA 


Subroutine  ARCTA  is  a  customized  arctangent  subroutine  for 
finding  the  angle  in  radians  from  the  arctangent  of  a  function 
(Louden,  1967,  p.  119).  The  library  arctangent  function  ATAN 
accepts  as  an  argument  the  tangent  of  an  angle  (sin/cos)  and  pro¬ 
duces  as  output  the  angle  in  radians.  Since  the  tangent  of  an 
angle  repeats  itself  every  180  degrees,  it  is  not  possible  to 
use  the  library  function  ATAN  to  determine  a  full  range  of  angles 
from  0  to  360  degrees.  To  compute  the  correct  angle  for  all  pos¬ 
sible  combinations  of  X  and  Y,  it  is  necessary  to  test  for  positive 
near  zero  and  negative  X,  and  positive  near  zero  and  negative  Y. 

The  IF  statements  accomplish  these  test  and  produce  an  angle  in 
radians  ranging  from  0  to  2  pi. 


HINDCAST  ANALYSIS  WITH  COASTAL  STORM  MODEL 


Hindcast  Tests  of  Model 

The  coastal  storm  simulation  model  can  be  used  to  hindcast 
wind,  wave  and  current  conditions  at  a  shore  site  durinq  the 
passage  of  a  coastal  storm.  Hindcast  analysis  differs  from 
forecast  analysis  discussed  in  the  previous  section  because 
exact  storm  positions  are  known  in  hindcasting,  whereas  a  con¬ 
stant  azimuth  and  storm  velocity  are  used  in  forecasting.  The 
results  of  hindcast  analysis  at  several  sites  are  included  in 
Appendix  C.  On  the  Great  Lakes,  the  sites  include  Holland  and 
Stevensvi 1 le,  Michigan,  and  Sheboygan,  Wisconsin.  On  the  east 
coast  of  the  United  States  and  Canada,  sites  include  the  Mag¬ 
dalen  Islands  on  the  Gulf  of  Saint  Lawrence;  Plum  Island,  Mass¬ 
achusetts;  Cedar  Island,  Virginia  and  Sapelo  Island,  Georgia. 

Mustang  Island,  Texas  was  studied  on  the  Gulf  Coast.  On  the 
west  coast  of  the  United  States,  hindcasts  were  made  for  Mon¬ 
terey,  California  and  South  Beach,  Oregon. 

Sites  were  selected  for  hindcast  analysis  which  had  weather 
and  wave  data  available  for  several  storms.  Several  of  the  sites 
were  studied  by  Davis  and  Fox  using  time  series  analysis  from  1969 
through  1975.  Other  sites  were  chosen  in  which  there  was  good  beach 
profile  data  which  could  be  correlated  with  wave  and  current  condi¬ 
tions  during  a  storm. 


Stevensville,  Michigan,  July  1969 

A  storm  which  passed  over  Lake  Michigan  in  late  July  1969 
has  been  choosen  as  an  example  of  hindcast  analysis.  When  the 
storm  passed  over,  a  30  day  time-series  study  was  being  conducted 
at  Stevensville,  Michigan  by  Fox  and  Davis  (1970a  and  b).  Stevens¬ 
ville  is  located  on  the  southeastern  shore  of  Lake  Michigan  about 
11  kilometers  south  of  Benton  Harbor,  Michigan.  The  shoreline  is 
oriented  roughly  north-south  with  an  average  nearshore  slope  of 
about  0.033. 

The  storm  which  affected  the  Stevensville  area  was  tracked 
from  2000  on  July  26,  1969  through  0800  on  July  30  (Table  3). 

The  size,  shape,  intensity  and  path  of  the  storm  were  interpreted 
from  weather  maps  for  July  26  through  30,  1969.  When  the  storm 
was  closest  to  the  coastal  site  at  Stevensville,  the  barometric 
pressure  at  the  center  of  the  low  was  estimated  as  994  millibars. 
The  pressure  at  the  largest  encircling  isobar  was  1012  millibars, 
and  therefore, the  maximum  pressure  included  in  the  storm  was  1014.6 
millibars.  The  storm  had  an  elliptical  shape  with  the  length  of 
the  major  half  axis  equal  to  960  kilometers  and  the  minor  half  axis 
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equal  to  700  kilometers.  The  orientation  of  the  major  half  axis 
was  30  degrees  west  of  north.  For  this  particular  example,  the 
equation  by  Fox  and  Davis  (1972)  was  used  to  compute  the  long¬ 
shore  current  velocity. 

The  position  of  the  shoreline  at  Stevensville  is  given  in 
the  X-Y  coordinate  system  with  X  equal  to  1043  and  Y  equal  to  486 
kilometers.  The  X  axis  runs  east-west  and  the  Y  axis  north-south 
with  the  origin  located  to  the  southeast  of  Stevensville.  The 
latitude  at  the  shore  location  is  42°  north  and  the  onshore  azi¬ 
muth  is  90  degrees.  The  nearshore  slope  from  the  shoreline  out 
across  the  nearshore  bars  is  0.033.  The  average  fetch  distance 
for  the  southeastern  shore  of  Lake  Michigan  is  about  200  kilo¬ 
meters  . 

The  storm  positions  were  plotted  at  6  hour  intervals  in  kilo¬ 
meters  on  the  X-Y  coordinate  system  (Table  3).  The  initial  posi¬ 
tion  of  the  storm  at  2000  on  July  26  was  to  the  northwest  of  the 
shore  site  with  X  equal  to  333  and  Y  equal  to  1229  kilometers. 

The  storm  passed  over  the  shoreline  about  0230  on  July  28.  At 
that  time,  the  storm  center  was  located  178  kilometers  north  of 
the  shore  site  (664-486  =  178  kilometers).  When  the  storm  track¬ 
ing  was  complete,  the  final  storm  position  at  0800  on  July  30  was 
X  equal  to  1907  and  Y  equal  to  1830  kilometers.  In  general,  the 
storm  made  a  loop  swinging  down  from  the  northwest,  passing  east¬ 
ward  across  the  shore,  then  moving  off  to  the  northeast. 

The  XI,  Y1  coordinate  is  oriented  with  the  XI  axis  parallel  to 
the  shore  and  the  Y1  axis  normal  to  the  coast  (Figure  2).  The 
origin  of  the  XI,  Y1  coordinate  system  is  at  the  center  of  the  storm 
with  the  positive  XI  direction  to  the  right  and  the  positive  Yl 
direction  toward  the  coast.  The  XI,  Yl  coordinate  system  is  used 
to  locate  the  shore  position  with  reference  to  the  center  of  the 
storm.  The  units  of  the  XI,  Yl  coordinate  system  are  in  terms  of 
storm  radii.  The  storm  radius  is  1.5  times  the  length  of  the 
major  half  axis  which  is  measured  from  the  center  of  the  storm 
to  the  largest  encircling  isobar.  Using  an  inverted  normal  curve 
to  simulate  the  storm  cross-section  the  largest  encircling  isobar 
is  defined  as  2  standard  deviations  from  the  center  of  the  storm. 

The  full  radius  would  extend  out  3  standard  deviations  from  the 
center  of  the  low.  As  the  storm  approaches  shoreline,  the  Yl  value 
decreases,  and  it  becomes  negative  after  the  storm  has  passed  over 
the  coast.  When  the  storm  is  to  the  north  of  shore  site  the  XI 
values  are  positive.  Therefore,  the  storm  at  Stevensville  remained 
to  the  north  of  the  shore  site  for  the  entire  run. 

In  the  hindcast  analysis,  the  barometric  pressure  decreases 
from  1012.4  millibars  on  July  26  to  a  minimum  of  995.2  at  2200  on 
July  27,  then  increased  to  1014.6  on  July  30,  1969  (Table  3  and 
Figure  10).  In  the  actual  barometric  pressure  record  at  Stevens¬ 
ville,  the  pressure  dropped  to  1000.2  millibars  (29.54  inches)  at 
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2000  on  July  27.  Since  small  scale  weather  maps  were  used  for 
determining  the  storm  position  and  pressure  values,  more  accurate 
results  could  be  obtained  by  using  large  scale  maps  available  at 
3  hour  intervals.  The  plot  for  barometric  pressure  is  a  function 
of  accurate  plotting  of  the  storm  positions  and  a  careful  estimate 
of  the  size  and  intensity  of  the  storm.  Therefore,  a  correspondence 
of  the  observed  and  hindcast  curves  for  barometric  pressure  is  not 
a  test  of  the  predictive  capabilities  of  the  model,  but  is  more  a 
test  of  the  accuracy  of  the  weather  maps  and  of  the  plotting  abil¬ 
ity  of  the  operator. 

In  the  wind  observations  at  Stevensville,  the  maximum  wind 
speed  during  the  storm  was  12.4  meters/second  (28.4  knots)  at  1600 
on  July  28  (Fox  and  Davis,  1970a).  For  the  hindcast,  the  maximum 
wind  speed  was  12.0  meters/second  at  1400  on  July  28.  The  overall 
pattern  for  the  observed  and  hindcast  winds  are  also  quite  similar 
(Figure  10 ) . 

The  observed  longshore  current  velocity  at  Stevensville  reached 
-116  centimeters/second  (northerly)  at  2200  on  July  27,  and  215  centi¬ 
meters/second  (southerly)  at  1800  on  July  28  (Fox  and  Davis,  1970a). 
The  maximum  hindcast  values  were  -62.9  centimeters/second  at  2100  on 
July  27  and  101.9  centimeters/second  at  1700  on  July  28  (Figure  10 
and  Table  3).  The  instantaneous  longshore  current  values  were  much 
higher  in  the  observed  than  the  hindcast  values,  however,  the  minimum 
and  maximum  values  occurred  within  one  hour  of  the  hindcast  values. 

The  observed  longshore  current  during  the  storm  was  exceptionally 
high  and  may  be  accounted  for  by  the  well  developed  trough  between 
the  nearshore  bar  and  the  shore  which  channeled  the  current  along 
the  coast.  The  value  used  for  nearshore  slope  could  also  be  too 
low  in  the  hindcast.  When  the  longshore  current  curve  was  smoothed 
using  the  15  term  Fourier  plot,  the  hindcast  more  closely  resembles 
the  observed  curve  for  longshore  current  (Fox  and  Davis,  1970a  and 
b). 


The  observed  and  hindcast  curves  for  breaker  height  are  very 
close  (Figure  10).  For  the  observed  curve,  the  maximum  breaker 
hieght  of  1.82  meters  (6  feet)  occurred  at  1800  on  July  28.  For 
the  hindcast,  the  maximum  height  was  1.83  meters,  also  at  1800  on 
July  28.  The  overall  shapes  of  the  observed  and  hindcast  curves 
for  breaker  height  are  also  very  similar  but  the  hindcast  curve 
drops  off  more  rapidly  than  the  observed  curve  (Figure  10). 
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Additional  Hindcast  Examples 


Additional  examples  of  the  hindcast  tests  are  presented  in 
Appendix  B  and  Figures  26,  27,  28  and  29.  Since  the  input  para¬ 
meters  are  given  in  Appendix  B  a, id  comparative  plots  are  in  Figure 
26  through  29,  a  full  discussion  will  not  be  included  for  each  of 
the  examples. 

The  observed  data  for  the  Holland,  Michigan  examples  was  ex¬ 
tracted  from  a  30-day  time  series  study  during  July  1970  (Fox  and 
Davis,  1971a).  Holland  is  located  on  the  southeastern  shore  of 
Lake  Michigan  about  96  kilometers  north  of  Stevensville.  The  hind¬ 
cast  values  for  longshore  current  and  breaker  height  were  quite 
close  for  July  3,  1970  (Figure  11  ),  but  were  somewhat  low  for  July 
18,  1970  (Figure  12  ). 

The  observed  data  for  the  Sheboygan,  Wisconsin  examples  come 
from  a  30-day  time  series  study  conducted  at  Sheboygan  during  July 
1972  (Fox  and  Davis,  1973a).  Sheboygan  is  located  on  the  western 
shore  of  Lake  Michigan  about  72  kilometers  north  of  Milwaukee, 
Wisconsin.  At  this  location,  the  storms  moved  in  a  northeast  di¬ 
rection  and  generally  offshore.  Therefore,  the  characteristic  re¬ 
versal  of  longshore  current  direction  observed  on  the  eastern  shore 
of  Lake  Michigan  was  not  observed  at  Sheboygan.  The  curves  for  July 
16  show  generally  low  observed  and  hindcast  curves  for  longshore 
current  and  breaker  height  with  pronounced  peaks  on  July  17,  1972 
(Figure  13).  The  curves  for  July  22  show  a  substantial  drop  in 
barometric  pressure,  but  very  low  values  for  longshore  current  and 
breaker  height  (Figure  14).  There  is  a  reversal  in  longshore  cur¬ 
rent  direction  on  July  24  for  both  the  observed  and  hindcast  curves. 
Since  the  winds  were  blowing  predominately  offshore,  the  waves  and 
longshore  currents  are  quite  subdued. 

Several  additional  tables  of  hindcast  results  are  given  in 
Appendix  B.  The  output  for  the  Atlantic,  the  Gulf  and  the  Pacific 
coasts  of  the  United  States  includes  a  variety  of  conditions  for 
storms  of  varying  sizes,  shapes  and  intensities.  Tidal  predictions 
are  also  included  for  the  oceanic  sites  where  the  tide  tables  are 
available.  Some  of  the  hindcast  results  followed  quite  closely 
with  the  observed  data,  while  at  other  places,  the  fit  was  not  as 
good  as  expected. 
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Figure  11.  Observed  and  hindcast  curves  for  Holland,  Michigan,  July 
3-5,  1970. 
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Figure  13.  Observed  and  hindcast  curves  for  Sheboygan,  Wisconsin, 
July  16-18,  1972. 


Figure  14.  Observed  and  hindcast  curves  for  Sheboygan,  Wisconsin, 
July  22-26,  1972. 
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FORECAST  ANALYSIS  WITH  THE  COASTAL  STORM  MODEL 


Short-term  Forecasts 


Short-term  predictions  of  wave  height  and  longshore  current 
velocity  can  be  made  using  the  forecast  mode  of  the  coastal  storm 
model.  Given  the  initial  coordinates  of  the  storm,  its  size,  shape 
and  intensity,  along  with  the  position  and  orientation  of  the  shore¬ 
line,  it  is  possible  to  predict  the  wind,  waves  and  currents  as  the 
storm  passes  over  the  coast.  In  the  forecast  mode,  the  storm  azi¬ 
muth  and  velocity  are  used  to  plot  a  straight  storm  path  as  the 
storm  proceeds  toward  the  shore.  In  general,  the  forecast  is  run 
for  72  hours  which  is  the  usual  limit  for  short-term  weather  pre¬ 
diction.  Examples  of  output  using  the  forecast  mode  are  given  in 
Appendix  C. 

In  the  forecast  mode,  it  is  assumed  that  the  size,  shape  and 
intensity  of  the  storm  remain  constant,  as  well  as  the  direction  of 
the  storm  path  and  the  speed  of  the  storm  along  the  path.  The 
operators  experience  with  weather  prediction  plays  an  important 
role  in  estimating  the  speed  and  path  of  the  storm.  For  12  to  24 
hours,  the  speed  and  path  may  remain  fairly  constant,  but  for 
longer  periods  of  time,  the  storm  may  veer  off  on  another  path 
or  change  its  speed  along  the  path.  Since  it  is  almost  impossible 
to  predict  the  path  of  a  storm  for  several  days,  a  series  of  dia¬ 
grams  have  been  devised  for  predicting  wave  and  current  conditions 
for  a  storm  with  a  given  size,  shape  and  intensity,  but  without  a 
fixed  storm  path. 

Circular  Storm  Test 


A  circular  storm  test  is  used  as  an  example  to  explain  how 
the  forecast  mode  works  with  the  coastal  storm  model.  The  cir¬ 
cular  storm  test  is  based  on  a  series  of  intense  storms  which 
crossed  over  the  Oregon  coast  during  late  fall  of  1973  (Fox  and 
Davis,  1974).  The  size,  shape  and  intensity  of  the  storm,  and 
the  orientation  of  the  shoreline  are  similar  to  those  encountered 
at  South  Beach,  Oregon  in  November  1973  (Fox  and  Davis,  1974). 

The  computer  listings  for  the  circular  storm  are  given  in  Ap¬ 
pendix  C. 

For  the  circular  storm,  the  barometric  pressure  at  the  center 
of  the  low  was  set  at  1000  millibars.  The  pressure  at  the  largest 
encircling  isobar  was  placed  at  1020  millibars  to  give  a  ranqe  of 
20  millibars  within  the  central  portion  of  the  storm.  The  plot  of 
barometric  pressure  for  the  circular  storm  model  is  generated  by 
rotating  an  inverted  normal  curve  around  its  center.  Therefore, 
the  barometric  pressure  surface  has  a  basin  shape  with  the  low 
pressure  at  the  center,  the  steepest  pressure  gradient  at  one 
standard  deviation  out  from  the  center,  and  gradually  reaches  a 


maximum  pressure  at  3  standard  deviations  away  from  the  center. 

Since  the  outer  margin  of  a  storm  often  interfers  with  other  high 
or  low  pressure  systems,  it  is  assumed  that  the  largest  encircling 
isobar  occurs  at  2  standard  deviations  away  from  the  center,  or 
2/3  of  the  total  storm  radius.  Therefore,  the  total  storm  radius 
would  be  1.5  times  the  radius  measured  at  the  largest  encircling 
isobar,  and  the  maximum  pressure  at  the  margin  of  the  storm  would 
be  1.145  times  the  pressure  range  from  the  center  to  the  largest 
isobar.  In  the  circular  storm  test,  the  pressure  at  the  2  stan¬ 
dard  deviations  is  1020  and  the  maximum  pressure  included  in  the 
storm  is  1022.9  millibars  at  3  standard  deviations. 

The  size  of  the  storm  is  determined  by  measuring  the  lengths 
of  the  major  and  minor  half  axes  of  the  storm  For  an  elliptical 
or  wave-shaped  storm,  the  length  of  the  major  half  axis  is  mea¬ 
sured  from  the  center  of  the  low  to  the  largest  isobar,  where  it 
is  farthest  from  the  storm  center.  The  length  of  the  minor  half 
axis  is  measured  at  right  angles  to  the  major  half  axis  from  the 
storm  center  to  the  largest  isobar.  A  circular  storm  exists  when 
the  major  and  minor  half  axes  are  equal.  For  the  circular  storm 
test,  lengths  of  the  major  and  minor  half  axes  are  300  kilometers. 
When  the  major  and  minor  half  axes  have  different  lengths,  the 
orientation  of  the  major  half  axis  is  plotted  in  degrees  from 
north.  Therefore,  for  the  circular  storm  test,  the  orientation  of 
the  major  half  axis  is  0  degrees. 

For  the  circular  storm  model,  the  storm  velocity  is  set  at  40 
kilometers/hour  with  a  storm  azimuth  of  90  degrees.  This  means  that 
the  storm  will  proceed  from  its  initial  position  along  a  path  90 
degrees  east  of  north  at  40  kilometers/hour.  Therefore,  if  the 
storm  is  tracked  for  30  hours,  it  will  move  a  distance  of  1200  kilo¬ 
meters  . 

For  the  shore  position  coordinates  in  the  circular  storm  test, 

X  is  set  at  1000  and  Y  is  set  at  0  kilometers.  In  the  X-Y  coor¬ 
dinate  system  used  for  plotting  storm  position  and  shore  location 
on  a  weather  map,  the  shore  site  is  located  1000  kilometers  east 
along  the  X  axis  and  0  kilometers  north  along  the  Y  axis.  The 
shore  latitude  for  the  test  case  is  43  degrees  north,  the  approxi¬ 
mate  latitude  of  the  Oregon  coast.  The  onshore  azimuth  is  90 
degrees  which  indicates  that  the  coastline  runs  north-south  with 
land  to  the  east  and  sea  to  the  west,  the  same  orientation  as  the 
Oregon  coast.  The  nearshore  slope  is  at  0.033  which  is  close  to 
the  average  nearshore  slope  in  the  different  areas  studied.  A 
value  of  1000  kilometers  was  used  for  the  average  fetch  which  would 
indicate  open  ocean.  If  the  average  fetch  is  greater  than  the 
storm  radius,  the  fetch  distance  will  not  have  an  effect  on  the 
wave  and  current  calculations.  However,  if  the  fetch  distance  is 
smaller  than  the  storm  radius,  the  waves  will  be  fetch  limited  when 
the  average  fetch  is  less  than  the  effective  fetch. 
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In  compiling  the  output  for  the  circular  storm  test,  the  storm 
outlined  above  was  tracked  along  12  paths  normal  to  the  shoreline 
and  parallel  to  each  other  (Figure  15).  The  location  of  the  shore 
site,  orientation  of  the  shoreline,  and  three  storm  tracks  are 
given  in  Figure  15A.  The  diagram  extends  1200  kilometers  in  a 
north-south  direction,  600  kilometers  to  the  north  (positive)  and 
600  kilometers  to  the  south  (negative)  of  the  study  site  which 
is  located  in  the  center  of  the  diagram.  The  diagram  also  extends 
1200  kilometers  in  an  east-west  direction,  600  kilometers  offshore 
(negative)  and  600  kilometers  onshore  (positive)  of  the  shore  site. 

A  time  scale  is  included  along  the  bottom  of  the  diagram  to  indicate 
the  length  of  time  in  hours  from  the  initial  tracking  of  the  storm 
to  any  position  along  the  storm  path.  The  speed  of  the  storm  was 
set  at  40  kilometers/hour,  therefore,  if  the  storm  started  at  the 
left  edge  of  the  diagram,  its  center  would  pass  over  the  shoreline 
after  15  hours  and  would  move  off  the  right  side  of  the  diagram 
after  30  hours. 
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BAROMETRIC  PRESSURE  -  MILLIBARS 


Figure  15.  Map  of  storm  tracks  and  time-distance  plot  of  barometric 
pressure  for  a  circular  storm. 


Barometric  Pressure 


The  barometric  pressure  diagram  (Figure  15B)  is  set  up  with 
the  same  coordinate  system  as  the  storm  track  diagram  (Figure  15A). 
The  barometric  pressure  which  would  be  recorded  at  the  shore  site 
is  plotted  at  the  storm  location  as  the  storm  moves  across  the  dia¬ 
gram.  After  12  storms  were  tracked  normal  to  the  shore  at  100 
kilometer  spacings  along  the  shore,  the  barometric  pressure  values 
plotted  along  the  storm  tracks  were  contoured  to  produce  the  baro¬ 
metric  pressure  diagram  (Figure  15B). 

The  use  of  the  barometric  pressure  diagram  can  be  explained 
by  examining  a  series  of  profiles  through  the  diagram  (Figure  16). 
Three  storm  tracks  are  plotted  on  the  pressure  diagram  indicating 
storms  which  moved  from  west  to  east  across  the  shoreline.  Baro¬ 
metric  pressure  profiles  are  shown  when  the  storm  track  is  200 
kilometers  north,  directly  over  the  shore  site,  and  200  kilometers 
south  of  the  site  (Figure  16).  When  the  storm  passes  200  kilo¬ 
meters  north  of  the  shore  site,  the  pressure  at  the  shore  location 
drops  from  1022.8  millibars  to  1013.4,  then  increases  again  to 
1022.8.  When  the  storm  passes  directly  over  the  shore  location, 
the  pressure  at  the  shore  site  drops  from  1022.8  millibars  to 
1000.0,  then  increases  to  1022.8.  When  the  storm  passes  200  kilo¬ 
meters  south  of  the  shore  site,  the  pressure  profile  is  identical 
to  the  profile  which  was  made  200  kilometers  north  of  the  shore 
site.  For  a  circular  storm,  the  barometric  pressure  diagram  is 
symmetrical,  so  that  profiles  cut  through  the  storm  a  given  dis¬ 
tance  north  or  south  of  the  shore  site  will  result  in  identical 
patterns. 

It  should  be  emphasized  that  the  values  plotted  at  the  storm 
locations  are  for  observations  recorded  at  the  shore  site  when 
the  storm  follows  along  a  given  path.  Therefore,  the  pressure 
profile  in  Figure  16  are  profiles  of  the  pressure  at  the  shore 
site  when  the  storm  passes  to  the  north,  over  the  site,  or  to  the 
south.  Although  the  diagram  for  barometric  pressure  is  identical 
to  a  weather  map  with  storm  center  located  directly  over  the  shore, 
it  should  not  be  interpreted  in  that  way.  The  X  axis  represents 
time,  and  distance  along  the  X  axis  is  used  to  show  the  storm  po¬ 
sition  at  a  given  time.  With  the  other  diagrams,  such  as  wind 
speed  and  breaker  height,  it  is  impossible  to  use  the  time  distance 
diagram  as  a  map  with  the  storm  center  located  at  the  shore  site. 
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Wind  Speed  and  Direction 

The  time  distance  diagram  for  surface  wind  speed  is  not 
symmetrical  because  the  maximum  wind  speeds  occur  when  the  storm 
track  is  between  100  and  200  kilometers  north  of  the  study  site 
(Figure  17).  When  the  storm  passes  200  kilometers  north  of  the 
study  site,  the  surface  wind  speed  reaches  28  meters/second  at 
the  shore  site.  However,  when  the  storm  passes  directl>  ever 
the  shore  site,  the  surface  wind  speed  reaches  18.8  meters/sec¬ 
ond  as  the  storm  approaches,  drops  down  to  zero  as  the  center  of 
the  storm  passes  over  the  coast,  then  increases  to  3C.P  meters/ 
second.  When  the  storm  track  is  200  kilometers  to  the  south  of 
the  shore  site,  the  surface  wind  speed  reaches  17.0  meters/second. 

The  highest  wind  velocities  are  recorded  when  the  storm  passes  to 
the  north,  and  after  the  storm  has  passed  over  the  coastline. 

To  understand  the  wind  pattern  during  a  coastal  storm,  it  is 
necessary  to  consider  both  wind  speed  and  direction.  Time-distance 
diagrams  are  plotted  for  surface  wind  speed,  wind  direction,  on¬ 
shore  component,  and  alongshore  component  of  the  wind  in  a  circular 
storm  (Figures  18A,  B,  C  and  D).  The  surface  wind  speeds  in  the 
storm  were  computed  by  the  geostrophic  wind  equation  (Equation  6) 
with  corrections  applied  for  speed  and  direction  (Equations  10  and 
11)  to  account  for  the  frictional  effects  of  land  or  sea  (Figure  5). 
For  example,  at  40°  north  latitude,  the  angle  between  the  surface 
wind  and  the  geostrophic  wind  would  be  42c  if  the  wind  is  from  the 
land,  and  16°  if  the  wind  is  from  the  sea  (Figure  5  and  Table  1). 

When  the  storm  passes  to  the  north  of  the  shore  site,  the  wind 
blowing  around  the  storm  center  in  a  counterclockwise  direction  is 
generally  onshore  at  the  shore  site  (Figure  18B).  The  contour  lines 
indicating  wind  direction  radiate  out  from  the  center  of  the  dia¬ 
gram  and  the  arrows  along  each  contour  line  point  in  the  direction 
the  wind  is  blowing  along  that  line  (Figure  18B).  The  dark  lines 
on  the  diagram  indicate  the  major  wind  directions  with  onshore  winds 
=  0°,  north  winds  =  90°,  offshore  winds  =  180°  and  south  winds  - 
270°.  When  a  storm  follows  a  path  200  kilometers  north  of  the  study 
site,  the  wind  direction  at  the  study  site  starts  out  from  the  south 
(274°),  slowly  shifts  over  to  onshore  (0°),  and  ends  up  out  of  the 
northwest  (52°).  When  a  sto"m  passes  200  kilometers  to  the  south  of 
the  shore  site,  the  wind  starts  off  blowing  offshore  (209  )  and  then 
shifts  over  to  the  north  (71°).  Different  patterns  are  used  to  show 
area  where  the  winds  are  blowing  generally  onshore,  from  the  north, 
offshore  are  from  the  south.  These  patterns  related  to  similar  areas 
for  the  onshore  and  alongshore  components  of  the  wind  (Figures  18C 
and  0). 

Diagrams  for  the  onshore  and  alongshore  components  of  the  sur¬ 
face  wind  speed  are  given  in  Figures  18C  and  D.  The  onshore  com¬ 
ponent  of  the  surface  wind  is  obtained  by  taking  the  ccsine  of  the 
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wind  direction  times  the  wind  speed,  and  the  alongshore  component 
is  produced  by  taking  the  sine  of  the  wind  direction  times  the 
wind  speed.  Storms  which  follow  a  path  to  the  north  of  the  shore 
site  generally  have  a  strong  onshore  wind  component,  while  those 
following  a  path  to  the  south  of  the  shore  site  are  predominately 
offshore  (Figure  18C).  The  dividing  line  between  the  onshore  and 
offshore  components  of  the  wind  follows  north-south  wind  direction 
lines  (0J  to  180“ )  (Figure  18B).  The  onshore  wind  reaches  a  max¬ 
imum  of  28.2  meters/second  along  a  storm  path  100  kilometers  north 
of  the  shore  >ite.  The  offshore  wind  reaches  18.6  meters/second 
when  the  storm  passes  100  kilometers  south  of  the  shore  site  and 
before  its  center  moves  across  the  coast.  The  onshore  wind  speed 
is  greater  than  the  offshore  wind  speed  because  the  friction  is 
less  when  the  wind  is  blowing  over  the  water. 

Tne  time  distance  diagram  for  alongshore  wind  component  in¬ 
dicates  that  the  winds  are  from  the  south  as  the  storm  approaches 
the  coast  and  shift  over  to  the  north  after  the  storm  passes  (Fig¬ 
ure  182).  To  the  north  of  the  shore  site,  the  shift  in  the  along¬ 
shore  component  from  south  to  north  takes  place  after  the  storm 
has  passed  over  the  coast,  but  to  the  south,  the  shift  takes  place 
before  the  storm  reaches  the  coast.  The  boundary  line  between  the 
north  and  south  components  of  the  longshore  wind  (Figure  18D)  fol¬ 
lows  the  onshore-offshore  line  (90°-270°)  in  the  wind  direction 
diagram  (Figure  18B).  The  reversal  from  south  to  north  wind  is 
very  abrupt  near  the  center  of  the  storm,  and  more  gradual  near 
its  margin.  When  the  storm  passes  to  the  north,  the  low  in  baro¬ 
metric  pressure  reaches  a  minimum  when  the  storm  passes  over  the 
coast  (Figure  17),  but  the  shift  in  wind  direction  from  south  to 
north  does  not  take  place  until  several  hours  after  the  low  has 
passed.  This  lag  in  wind  direction  reversal  behind  the  low  in 
pressure  was  observed  during  storms  at  Holland  and  Stevensville, 
which  passed  to  the  north  of  the  study  areas  (Fox  and  Davis,  1970 
and  1971a).  Also,  the  maximum  wind  speeds  were  observed  after 
the  lows  had  passed  and  the  wind  shifted  over  to  the  north. 
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gure  18.  Time-distance  plots  of  A  -  surface  wind  speed,  B  -  wind 
direction,  C  -  onshore  wind,  and  D  -  alongshore  wind  in 
a  circular  storm. 
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Wave  Period  and  Breaker  Height 

The  diagrams  for  wave  period  and  breaker  height  are  quite 
similar  (Figures  21 A  and  B) ,  and  both  resemble  the  plots  for 
surface  wind  and  onshore  wind  (Figures  18A  and  C).  Three  pro¬ 
files  were  plotted  across  the  breaker  height  diagram  to  show 
what  the  height  would  be  at  the  shore  site  as  the  storm  moves 
across  (Figure  19).  If  the  storm  moves  200  kilometers  north 
of  the  shore  site,  the  breaker  height  will  reach  5.30  meters 
3  hours  after  the  storm  crosses  the  coast.  If  the  storm  passes 
directly  over  the  shore  site,  the  breaker  height  will  reach 
.86  meters  as  the  storm  approaches,  drop  down  as  the  center  of 
the  storm  passes,  then  reaches  3.63  meters  5  hours  after  the 
center  passes  over  the  coast.  The  decrease  in  wave  height  as 
the  storm  center  passes  directly  over  the  shore  site  corresponds 
to  the  zero  wind  velocity  at  the  center  of  the  storm.  While 
wind  velocity  may  drop  to  zero  at  the  storm  center,  the  zero 
wave  height  is  probably  an  artifact  of  the  computer  model  and 
does  not  occur  in  nature.  Residual  waves  would  most  likely  re¬ 
main  in  the  area  and  could  be  built  model  if  so  desired. 

The  surface  wind  speed  is  not  used  directly  for  determining 
breaker  height  and  wave  period,  because  strong  onshore  wind  is 
effective  in  generating  waves  which  will  reach  the  coast,  and  a 
strong  offshore  wind  tends  to  subdue  existing  waves.  On  the  Texas 
coast  during  studies  made  at  Mustang  Island,  offshore  wind  was 
accompanied  by  a  sharp  drop  in  breaker  height  (Davis  and  Fox,  1972c). 
On  Lake  Michigan,  where  a  single  storm  system  was  studied  as  it 
moved  offshore  at  Zion,  Illinois  and  onshore  at  South  Haven,  Mich¬ 
igan,  breaker  height  was  over  2  times  as  great  where  wind  was 
blowing  onshore  than  where  it  was  blowing  offshore  (Davis  and  Fox, 
1974b).  Therefore,  the  effective  wind  speed  was  used  in  deter¬ 
mining  wave  height  and  period  in  place  of  the  surface  wind.  For 
an  onshore  wind,  the  effective  wind  is  equal  to  the  onshore  wind. 
However,  for  an  offshore  wind,  the  effective  wind  is  about  one 
third  of  the  surface  wind  speed.  For  a  wind  blowing  along  the 
shore,  the  effective  wind  speed  is  two  thirds  of  the  surface  wind 
speed.  A  cosine  transformation  was  used  to  produce  a  smooth  gra¬ 
dient  in  effective  wind  from  onshore  through  alongshore  to  off¬ 
shore. 

The  plot  for  wave  period  closely  resembles  the  plot  for  breaker 
height  with  the  maxima  to  the  north  of  the  shore  site  and  displaced 
landward  of  the  shoreline  (Figure  21 A ) .  The  maximum  wave  period  of 
9.5  seconds  occurs  at  the  same  time  as  the  maximum  breaker  height, 

3  hours  after  the  storm  has  passed  over  the  coast.  The  plot  for 
wave  period  has  a  broad  relatively  flat  area  surrounding  the  maxi¬ 
mum,  while  the  plot  for  breaker  height  is  much  steeper,  reaching  a 
peak  and  rapidly  dropping  off  after  the  peak  has  been  passed.  The 
wave  height  and  periods  forecast  in  the  model  correspond  closely  to 
those  encountered  on  the  Oregon  coast  during  November  1973  (Fox  and 
Davis,  1974). 
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Breaker  Angle  and  Longshore  Current  Velocity 


The  plots  of  breaker  angle  and  longshore  current  velocity 
(Figures  21C  and  D)  are  similar  in  many  respects  to  the  plot  for 
alongshore  wind  (Figure  180).  The  boundary  line  which  separates 
the  north  and  south  components  of  the  wind  is  the  same  as  the 
boundary  which  separates  the  north  and  south  breaker  angles  and 
longshore  currents.  For  field  studies  conducted  at  Holland  and 
Stevensville,  Michigan,  there  was  also  a  close  correspondence 
between  the  longshore  component  of  the  wind  and  longshore  current 
velocity  (Fox  and  Davis,  1970a  and  b,  and  1971a). 

The  breaker  angle  is  defined  as  the  acute  angle  between  the 
wave  crest  and  the  shoreline  as  the  wave  passes  over  the  nearshore 
bar.  In  deep  water  the  wave  direction  is  roughly  parallel  to  the 
wind  direction,  and  the  wave  crests  are  about  normal  to  the  wind. 
The  dominate  wind  direction  is  often  used  in  wave  refraction  com¬ 
puter  programs  to  determine  the  deep  water  wave  angle  (Dobson, 
1967).  As  a  wave  enters  shallow  water,  the  celerity  decreases  and 
the  wave  crest  is  refracted  so  that  it  becomes  closer  to  parallel 
to  the  beach.  Snell's  Law  of  geometrical  optics  is  used  for  com¬ 
puting  the  refraction  coefficient  and  breaker  angle  in  the  surf 
zone. 


In  the  breaker  angle  diagram,  the  area  to  the  left  of  the  zero 
line  has  breaker  angles  open  to  tue  north,  and  to  the  right  of  the 
zero  line,  the  breaker  angles  are  open  to  the  south  (Figure  21C). 

The  largest  breaker  angles  are  about  32  degrees  when  the  wind  is 
blowing  directly  out  of  the  north  or  the  south.  As  the  wind  di¬ 
rection  swings  around  from  alongshore  to  onshore  or  offshore,  the 
breaker  angles  decreases  from  30  degrees  to  zero  (Figures  18B  and 
21C).  When  the  wind  is  blowing  directly  onshore  or  offshore,  the 
breakers  are  parallel  to  the  beach  and  breaker  angle  is  zero. 

The  plot  for  longshore  current  velocity  is  very  similar  to 
the  plot  for  the  alongshore  component  of  the  wind  (Figure  ISO  and 
21D).  Three  profiles  are  plotted  when  the  storm  passes  200  kilo¬ 
meters  north,  over  the  shore  site  and  200  kilometers  south  (Fig¬ 
ure  20).  When  the  storm  path  is  200  kilometers  north,  the  long¬ 
shore  current  velocity  reaches  97.4  centimeters/second  to  the  north, 
reverses  direction  after  the  low  has  passed,  and  increases  to  81.6 
centimeters/second  to  the  south.  When  the  storm  passes  directly 
over  the  study  site,  the  current  reverses  from  63.0  centimeters/ 
second  to  the  north  to  129.7  centimeters/second  to  the  south.  How¬ 
ever,  when  the  storm  path  is  200  kilometers  south  of  the  shore  site, 
the  current  to  the  north  is  only  4.8  centimeters/second  anu  t no- 
southward  current  is  58.1  centimeters/second.  When  the  stcni;  passes 
to  the  south,  the  reversal  in  current  direction  takes  place  before 
the  low  in  barometric  pressure  passes  the  shore.  The  maximum  long- 
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Figure  19.  Time-distance  plot  of  breaker  height  and  three  profiles  j" 

of  breaker  height  in  a  circular  storm. 
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Figure  20.  Time-distance  plot  of  longshore  current  and  three  profiles 
of  longshore  current  in  a  circular  storm. 
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Time-distance  plot  of  A  -  wave  period,  B  -  breaker  height, 
C  -  breaker  angle,  and  D  -  longshore  current  in  a  circular 
storm. 
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shore  current  was  144  centimeters/second  to  the  south  4  hours  after 
the  storm  passed  over  the  shoreline  on  a  path  100  kilometers  north 
of  the  shore  site  (Figure  20). 

The  longshore  current  velocity  is  a  function  of  nearshore  slope, 
breaker  height  and  breaker  angle  (Longuet-Higgins ,  1970).  The  in¬ 
fluence  of  both  breaker  angle  and  breaker  height  can  be  seen  in  the 
plot  for  longshore  current  velocity  (Figures  19B,  C  and  D).  Other 
longshore  current  equations  were  tested  which  gave  similar  patterns, 
but  different  absolute  velocities. 

The  circular  storm  test  illustrates  the  general  patterns  which 
emerge  in  barometric  pressure,  wind,  waves  and  longshore  currents 
as  a  storm  passes  over  a  coast.  If  the  shape  and  path  of  the  storm, 
and  the  orientation  of  the  shoreline  are  held  constant  while  the 
size  or  intensity  of  the  storm  are  varied,  the  same  patterns  will 
persist,  but  the  absolute  values  will  change  for  each  of  the  vari¬ 
ables.  However,  if  the  shape  or  path  of  the  storm  are  changed, 
the  patterns  as  well  as  the  absolute  values  will  change  for  each 
of  the  variables.  In  the  next  section  an  elliptical  storm  is  used 
to  demonstrate  the  effect  which  a  change  in  storm  shape  would  have 
on  the  wind,  waves  and  currents. 
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Elliptical  Storm  Test 

With  the  coastal  storm  model,  it  is  possible  to  vary  the 
size  and  shape  of  the  storm  while  holding  the  intensity  constant. 
The  elliptical  storm  test  provides  a  good  example  of  an  oval  shaped 
storm  which  has  its  long  axis  extending  to  the  north-northeast 
(Figure  22).  In  both  the  circular  and  elliptical  tests,  the  shore¬ 
line  orientation  and  nearshore  bottom  slope  are  the  same. 

In  the  circular  and  the  elliptical  storms,  the  minimum  baro¬ 
metric  pressure  at  the  center  is  1000  millibars,  and  the  pressure 
at  the  largest  encircling  isobar  is  1020  millibars  (Figures  16  and 
22).  In  the  circular  storm,  the  major  and  minor  axes  are  the  same, 
300  kilometers.  In  the  elliptical  storm,  however,  the  major  axis 
(500  kilometers)  is  twice  the  length  of  the  minor  axis  (250  kilo¬ 
meters).  The  major  axis  in  the  elliptical  storm  is  oriented  30c 
east  of  north.  Therefore,  a  low  pressure  trough  extends  in  a 
north-northeast  direction  with  the  lowest  value  at  the  storm 
center . 

The  time-distance  plot  of  barometric  pressure  is  identical 
to  a  weather  map  made  when  the  storm  center  is  over  the  shore 
site  (Figure  22).  When  the  storm  track  is  located  to  the  south 
of  the  shore  site,  the  low  pressure  trough  reaches  the  shore  site 
before  the  low  pressure  center  passes  over  the  coast.  However, 
when  the  storm  track  is  to  the  north  of  the  shore  site,  the  low 
pressure  center  reaches  the  coast  before  the  trough  passes  over 
the  shore  site.  Therefore,  the  long  axis  of  the  storm  marks  the 
time  when  the  low  pressure  trough  passes  over  the  shoreline. 

Although  the  range  in  barometric  pressure  is  the  same  in  the 
circular  storm  and  the  elliptical  storm,  the  pressure  gradient  is 
steeper  in  the  constricted  part  of  the  elliptical  storm.  The  pres¬ 
sure  gradient  is  a  function  of  the  size  of  the  storm  and  the  range 
in  barometric  pressure.  In  the  circular  storm,  both  the  major  and 
minor  axis  have  lengths  of  300  kilometers,  and  therefore  the  pres¬ 
sure  gradient  is  equal  on  all  side  of  the  storm.  In  the  elliptical 
storm,  the  major  axis  is  500  kilometers  and  the  minor  axis  is  250 
kilometers.  Therefore,  along  the  minor  axis  the  pressure  gradient 
is  steeper,  while  it  is  more  gentle  along  the  major  axis. 

The  time-distance  plot  of  surface  wind  speed  has  an  elliptical 
shape  with  the  high  winds  concentrated  on  the  right  side  of  the 
diagram  (Figure  23A).  The  high  wind  speeds  are  a  function  of  the 
steeper  pressure  gradient  along  the  minor  axis  and  differences  in 
surface  friction  over  land  and  sea.  With  the  higher  pressure  gra¬ 
dient,  the  wind  speed  reaches  36.8  meters/second  in  the  elliptical 
storm,  while  in  the  circular  storm,  it  only  reaches  30.9  meters/ 
second.  The  winds  greater  than  20  meters/second  are  split  into 
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two  areas  in  the  elliptical  storm,  a  major  area  down  the  right 
side  of  the  storm,  and  a  minor  area  in  the  northwest  quadrant. 

At  the  center  of  the  storm  along  the  major  axis,  the  surface  wind 
speed  drops  down  to  zero. 

The  pattern  for  wind  direction  in  the  elliptical  storm  is 
significantly  different  from  the  pattern  in  the  circular  storm 
(Figures  18B  and  23B).  In  the  circular  storm,  the  wind  direction 
contours  radiate  out  from  the  center  and  are  rotated  in  a  clock¬ 
wise  direction  from  14°  to  45°.  The  zero  wind  direction  indicates 
onshore  wind  and  180°  is  an  offshore  wind.  The  elliptical  storm 
is  constructed  along  the  minor  axis  and  extended  along  the  major 
axis.  Therefore,  the  wind  direction  contours  are  gathered  around 
the  major  axis  which  is  pointed  30°  east  of  north.  Similarly, 
the  contours  are  spread  out  from  the  minor  axis.  In  the  elliptical 
storm  the  zero  contour  extends  45°  east  of  north,  and  the  180°  con¬ 
tour  extends  67°  west  of  south.  The  270°  contour  indicating  south 
winds  points  17°  east  of  north  and  the  90°  contour  for  north  winds 
extends  25°  west  of  south. 

The  difference  in  wind  direction  pattern  results  in  major 
changes  in  the  onshore  and  alongshore  wind  patterns  in  the  elliptical 
storm  (Figures  23C  and  D).  In  the  circular  storm,  the  boundary  line 
between  onshore  and  offshore  winds  runs  in  a  generally  east-west 
direction  with  onshore  winds  when  the  storm  track  is  to  the  north 
and  offshore  winds  when  the  storm  track  is  to  the  south  of  the  study 
site.  In  the  elliptical  storm,  on  the  other  hand,  the  boundary  bp- 
tween  onshore  and  offshore  winds  has  shifted  so  that  it  runs  general¬ 
ly  north-south  with  onshore  winds  to  the  east  and  offshore  winds  to 
the  west  (Figure  23C).  As  the  storm  approaches  the  coast  the  winds 
are  offshore,  and  after  the  storm  has  passed  over  the  shoreline  the 
winds  shift  to  onshore.  When  the  storm  passes  to  the  south  of  the 
shore  site,  the  shift  from  offshore  to  onshore  winds  takes  place 
shortly  after  the  low  pressure  trough  passes  over,  but  when  the 
storm  track  is  to  the  north,  the  shift  in  wind  direction  shortly 
preceeds  the  low  pressure  trough.  The  maximum  offshore  wind  is 
21.4  meters/second  as  the  storm  approaches,  and  the  maximum  onshore 
wind  is  30.2  meters/second  after  the  storm  has  passed  over  the  coast. 

The  boundary  for  the  alongshore  wind  which  separates  the  north 
wind  from  the  south  wind  extends  generally  in  a  northwest-southeast 
direction  (Figure  23D).  It  resembles  the  alongshore  wind  diagram 
for  the  circular  storm,  but  the  axes  are  rotated  about  30°  in  a 
clockwise  direction  (Figures  180  and  23D).  The  maximum  south  wind 
of  21.1  meters/second  occurs  when  the  storm  track  is  to  the  north 
of  the  shore  site  and  after  the  storm  has  passed  over  the  coast. 

For  the  north  wind,  the  maximum  of  28.1  meters/second  occurs  along 
a  storm  track  200  kilometers  south  of  the  shore  site  after  the  storm 
has  passed  the  coast. 
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Figure  22.  Time-distance  plot  of  barometric  pressure  in  an  elliptical 
storm. 
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Figure  23.  Time-distance  plot  of  A  -  surface  wind  speed,  B  -  wind 
direction,  C  -  onshore  wind  and  D  -  alongshore  wind  in 
an  elliptical  storm. 
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C.  ONSHORE  WIND  Q  ALONGSHORE  WIND 


Figure  24.  Time-distance  plot  of  A  -  wave  period,  B  -  breaker  height, 
C  -  breaker  angle  and  D  -  longshore  current  in  an  ellip¬ 
tical  storm. 


f 
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The  plots  for  wave  period  and  wave  height  in  the  elliptical 
storm  are  very  similar  and  closely  resemble  the  plots  for  wind  speed 
and  onshore  wind  (Figures  23A  and  C,  and  24A  and  B).  Both  wave 
period  and  breaker  height  reach  their  maxima  after  the  storm  trough 
has  passed  over  the  coast  and  the  wind  has  shifted  from  offshore 
to  onshore.  In  the  elliptical  storm,  the  maximum  wave  period  is 
8.9  seconds  and  the  greatest  breaker  height  is  5.1  meters. 

The  greatest  breaker  angles  for  the  elliptical  storm  occur 
when  the  wind  is  blowing  out  of  the  south  (270°)  or  the  north  (90°) 
(Figures  22B  and  23C).  When  the  storm  track  is  to  the  north  of 
the  shore  site,  the  greatest  northerly  breaker  angles  (32. 3‘ )  are 
present  just  before  the  low  pressure  trough  passes  over  the  coast. 
However,  when  the  storm  track  is  to  the  south  of  the  study  site, 
the  largest  southerly  breaker  angles  (32.4°)  occur  just  after  the 
storm  trough  has  passed  the  shore. 

The  plot  for  longshore  current  in  an  elliptical  storm  is  very 
similar  to  the  plot  for  alongshore  wind  (Figures  23D  and  24D).  The 
boundaries  between  north  and  south  winds  and  north  and  south  currents 
follow  the  same  line  and  the  maxima  are  in  the  same  position.  The 
maximum  northward  flowing  longshore  current  (44  centimeters/second) 
occurs  when  the  storm  is  on  a  track  300  kilometers  north  of  the  shore 
site,  while  the  maximum  southward  flowing  current  (143  centimeters/ 
second)  occurs  on  a  storm  track  100  kilometers  south  of  the  shore 
site. 


In  summarizing  the  comparison  between  a  circular  storm  and  an 
elliptical  storm  of  the  same  intensity,  the  differences  in  baro¬ 
metric  pressure  and  wind  direction  influence  the  other  environmental 
parameters.  In  the  circular  storm,  the  highest  surface  winds  are 
found  when  the  storm  track  is  to  the  north  of  the  shore  site  while 
in  the  elliptical  storm,  the  maximum  winds  occur  in  a  north-northeast 
trending  zone  to  the  right  of  the  low  pressure  trough.  The  boundary 
between  offshore  and  onshore  winds  lies  generally  east-west  for  the 
circular  storm  and  north-south  for  the  elliptical  storm.  For  along¬ 
shore  winds,  the  boundary  between  north  and  south  winds  is  rotated 
about  30°  in  a  clockwise  direction  in  the  elliptical  storm.  In  both 
the  circular  and  elliptical  storm,  the  wind  direction  contours  radi¬ 
ate  out  from  the  center  and  are  rotated  clockwise  due  to  surface 
friction,  but  in  the  elliptical  storm,  the  contours  are  gathered 
around  the  major  axis.  Wave  period  and  breaker  height  in  the  circ¬ 
ular  storm  resemble  wind  speed  and  form  pods  to  the  north  of  the 
shore  site,  while  in  the  elliptical  storm,  period  and  breaker  height 
form  linear  trends  to  the  east  of  the  shore  site.  For  both  storms, 
the  longshore  current  follows  the  same  pattern  as  the  alongshore  wind. 
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CONCLUSIONS 


A  mathematical  model  has  been  developed  and  programmed  for  a 
computer  to  forecast  barometric  pressure,  wind,  waves  and  longshore 
currents  during  passage  of  a  storm  across  a  coastal  site.  The  fol¬ 
lowing  set  of  conclusions  can  be  drawn  from  the  coastal  storm  model. 

1.  The  shape  of  a  coastal  storm  can  be  approximated  with  an  ellip¬ 
tical  model  by  specifying  the  lengths  of  the  major  and  minor  half 
axes  and  the  orientation  of  the  major  axis. 

2.  The  barometric  pressure  profiles  along  the  major  and  minor  axes 
of  the  ellipse  are  represented  by  a  series  of  inverted  normal  curves. 

3.  Geostrophic  wind  speed  and  direction  at  any  point  on  the  earth's 
sufrace  under  a  storm  are  computed  from  the  latitude  and  barometric 
pressure  gradient. 

4.  Geostrophic  wind  speed  and  direction  are  used  to  compute  surface 
wind  speed  and  direction  over  land  or  sea. 

5.  Wave  period,  height  and  direction  are  calculated  from  the  wind 
speed,  fetch  and  duration  as  the  storm  passes  over  the  coast. 

6.  Longshore  current  speed  and  direction  are  computed  from  wave 
height,  period  and  direction,  and  nearshore  bottom  slope. 

7.  Wave  and  current  hindcast  data  can  be  used  to  test  the  model  when 
the  size,  shape, intensity  and  track  of  the  storm  are  known. 

8.  If  the  storm  azimuth  and  velocity  are  assumed  to  be  constant,  fore¬ 
casts  of  wind,  wave  and  longshore  currents  can  be  made  for  a  coastal 
storm. 
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lard  a  -  riot  irediltIwN  -  fr_m  carl  ) 


if R  I  No  ’  IDE  RANGE  IN  FEET 
•»t*F  TIDE  RANGE  :  •«  ^lET 

Nl/MBlR  .1-  .  A»i  F  RiJM  LAST  SPRING  TIDE 

Ml  •.  F  WAAT  high  SPRING  T  l  _  t 


-  F,im  Nl'MBE  R 

T,.  M  I  L’  I  uR  N  A  l  TlDt 

,.N  ’.  :.*■  -  “  I  *  E  L  it  M  i  D  I  UR  N  Al  TIOt 
.  .  '  ■  ‘  .  -  m;  *tD  w  :  uR N Al.  7  IDE 

jRE*’*r  '‘'A-,  * .  -  Dl-RNAk  TIDE 


•  si  *-  «  -4  hl. 
S.  •  *♦!  -  •»}  A  fc  S '  •  ■  -i  *« 

w.  ft-,  -  .  ,  (  _  t 


Slope  at  iv.  TIDE 
v*t  A '  h!l»h  TILE 
*.  FEE’ 


I 


*-  .  gm.INAM  ;  N  NAjTJCA.  •*  |  i  l 
•  -  .  iSA’t  IS  SAw’itAv  M»ilS 

*•  H_-i  »  RUM  NORTH 

Sf»«‘-..K|  SLL.Pl 

-  L  S'  IS!'*’  4,  .  A  *,  *  O-R  B  An  k  1  E  **  ISLAND 

.  ~  -A.  “  S">t  IS.  ANL 

■  A  *-  m  .  •  |  .«  v.  •  PR*  .  *  ,  N 


1  -  s’  “*m  ;.;S',  a-  %  **OwR 

.S  1  R  ,M  ;  j  '  ;  al  -1,1. 

. -'  l  .  *  /*  -.  ."«»'(  is  •» a l ’  : 1  * i 

.  ^  »  -  -  .  .  ,  -  •  -  |  IN  « A  .  •  :  .  A  .  k;t(  i 

‘  -  R  .  -  1  k  *  A  i  •  I  n  J  -  .El  :t  i  •  •  .  A  i  ,  »v  ’  M  AS .  •  . s 

.  .  ,  -  -  i1  ,)  ,  i  *  *  ;  N  »  N  _  *  s 


L 


c 


CwLS.  0-U  -  STJ«m  A2!MyTrt  -  ClOOcISE  FROM  Nuk|m 

C  COLS.  ii>-2i  w(l)  -  IMTIAl  X-CClRDINATE  IN  NAUTICAL  M^ES 

C  COLS.  22-4*1  v(I)  -  INITIAL  v-COufcJiNATfc  IN  NAw.fJC.AL  «Uti 

CNA-1. 85319 
C«\«.?>396l 

C*r  •  J.2s08 
C fC  f  i  ■  1  •9‘*2t> 
tC'<r-  .  .032908 
CjDol  *  *737561 
CPU* 33. 86 39 
C 

RAD-  5  7.29  58 

8C90  READ  (2*915)  TUlL 

915  FORMAT ( 20A* ) 

«f*iTfcl3»91b  t  TITLE 
91®  FCRMat  (  INI // t  1*  »20A*/ ) 

H£  AC.  (  2  »  91  7  »  ISTRT  »DAr . 1 NAy T , nAoT 
917  FORMAT  l  12  »5An  »2  !  1  » 

IFilNAgT)  1.1.2 

1  CNN  ■  1 . 0 

F  M  *  1  •  0 

2  IFINAlT)  11.1  1. 12 
U  CRN-1. C 

CMF-i.o 
CRTS-1.0 
CCMF-l.O 
CjOLl-1. 0 
CPj-i.c 

12  IF  •  I STRT  .tiJ.0.0  J  siO  TO  1000 
-»l TE  »  3,  via  )  i  STM  1  .  0 A  I 

91fc  F  oS"A  T  |  IX  .  •  NWN  PtOlNB  AT  HOUR  1  .  J  ,  1  ON  *.5A*»I 

HEAD  i  2 .901  I  JNOPT  »  I  F  T  I U  »l  SCOP  •  N  A  *  'jnF  C  i  *  T  J  NT  .PM  IN  .PM  A.  < , Sl AT  .AR1.sk! 
1  .EA2 

90  1  FORMAT  I  31  1  *  13  .F6.0.F5.  .  ,2F  ,F  5.  JF6.0I 
•'NFCn.bNFC  1  «CMi 
A9«AR*B( S« 

BR  *  b  W 1 • »  Nk 

c 

SAHUA’t  TR  I  C  F  N  L*ioB  c.  AT  LARGEST  ENCIRClING  ISOBAR  lb  AiiuMtD  lw  0L 
C  *  T  f  »'J  i  T  A  N  0  A  k  J  L-tvlAT  f  V.-NS  FkC/M  IhE  CENTER  of  ThE  S  T  onM.  To 

C  FIND  ThL  A.rgAw  SI-.RA-  SUl.  ThE  majoK  AND  MINOR  AaES  ARE 

C  VyL  f  1**1.  Ir :  !3*  1.V  AND  Tm,  '*A«IH\,y  PRt.SSuRt  IS  MULTIPLIED  bir 

-  i  .t'rKvj.t,  TH  ••R\.SSgRt  A  -  ThE  MARGIN  of  ThE  STURM. 


•.  «*lAL  .  *  ,9?C.  Si  .  T  :  ;A  '  .  ThL  tf  \, -LPLC.S.  t-Hl  .  TMEAl 

>2  0  ►0R9»*  :  ’  •  2  •• 

TN- '  l*c‘ M 
TME  A  *«  *  T  vt A  .  *L  *  m 
»  RA'.si‘'<r  *‘£  <R-Pf*  J  N 

;,v  Aa-^.iW*!  ,  i«b«Pm|n 

«K  1  •  t.  I  3 . 905  I  PMlN.PMAxR  |M*» 

F'-R-at  (  ;mO  .  *ST  .>M  -  cARvMETHIC  PRESSURE  AT  CENTEh  OF  lOk  -'.F7.1. 
I'  M  ILL  ;tjA*$  ■  ,  9>  .  ’P^ESSLME  AT  iAKQEST  ENCIRCLING  1  SOBAh  ■••♦*•1. 

;•  “ ;  Ll  I  “A-i  ’  /  .  ' -AA  I  Mw'M  .'Rt'.SOPE  IN  L  oDED  IN  b  I  oRM  ’.A. 

3  *  m  l l  l  B  A  •  1  * 

A‘«  1  »  AR  .  J’> 

9  -  1  "  l  R  •  v  *  ’*♦•05 
i  ...  A  L  *  • 

^  p;  ' t  J.9.K.  A  R  I  .A  J  INN)  #A2  (  NN  i  .  A3  I  N  .)  »A*»  t  NN  1  *BR 1 • A  1  INN)  .  A  2  <  NN  I  » 

*  Aii  j  »  A  ••  l  N  .  E  **’ 

t  .  •  ’  I  »»  »  \  E  S  G  *  ”  QF  VA.J<  ►'A.  •  AXIS  •*  .  K7.1.4.AL/ 

.  •  •  .  '  » k  N  •»’  OF  v  I  N  w  n  hAl  F  A»!o  •'  .F  7.1.HAA, 

29*  .  •  ’  •  I  J.*.T  *  T  1  WN  .  >  “Aj.K  A»'i  7.  1  ♦  •  OEC.MEES  FROM  sjMTM'n 

'  •  (  81  Oi  .8  1 .8h  :  .LCCoP 

■ :  - - : * e  )i“iv 

5  1  .  .(.‘./.jy-  RL  CoNRfcF.T  (  rUA  T  |  ON  FROM  F^a  AND  D*  v  I  S  •  1VT2'"* 


“  2  *  -  .  '  l  :  .*  .  8  W 

•*l.  1  .rva  •  i  9« ,  •  l  _*«osm->£  current  e-oat;  <n  from  LONooiT-HiooiNb.  iv7j' 


9<  -  (• :  ’  e  <  s .  1 1  * 

A.  »  F  ..-.A  »  I  VA  .  1  ^  .RE  LORRtNT  Evr^ATJON  FROM  C.E.H.L.  197j*/l 

’  8  *. 

-•w  --.'E  3.8.6  1 

h;«  •  •* M  A  "  9<  .  •fcoN%»SMC«t  vwRREST  *  MOM  »  oMAR  AND  INMAN.  19?0*/l 

•J5  jN  t  J  Swl 

AB  AS.  e-.  *.|  -Ai.il  OF  IHL  S  T  Cam  Ow»’  T^  T  wv  STANDARD  0 1  v  1  A 1  I  ONS .  .  . 
*■-  .Tl'sT  S’  , .!)  TC  rItlw  Tm)  *ua  lENOThS. 

A  -  1 , s*AR 

r4  *  i  .  5#tJR 

RFALt2.9C  ■  -vO  .v^Dt  .Sha2.Sl'-PF  .  I  Si  NO.BLOR 
<#:  •  »■;«>»•  u1  '.0tl2»F«.ii 
‘  M A  »  -OMA A*  RC  OP 

-w^C  /C*v.NA 

>c-/w  .<*CNA 

IF  1  is jPT-1  )  1 jUD.10.J0 

-.All  SoBROoTjnES  F  or  mJNDCASTInc.  AND  F-jRjlAST  ]  NO  STORM 

p . S  J  T ; VNS 

1  ',*n  H  I  NO  <  N  »  .  A  i 

J  J  ’  >  6  C 

}  C  Aw4.  »  J«E  '  S«  .A  «NN 


c 

c 

c 


*0  00  *2  I • 1 iN* 

01  1  >  “U  l  l  )  »lNK. 

4 2  V  (  l  »  •  V  .  1  >  -CNH. 

OLOC  *UlOC / A 
VwOC-VLOC/A 

1  > • A* • 006 
Vl«V l 1  I  • • 006 
OL»ULOC»A 
VL»vlOC*a 

■SITE  SHORE  AND  TlOt  DAT* 


ULl»UL»CAN*.05 
VLl-Vl  *<XN-*#05 
8NECl«BNfCH*CKN 

MR1 TE  I 3.970  »UL1 .VL1 « AH NN )  *A2 <  NN ) tA3 i  NN ) lAdiMM  *SlAT • SnAZ  ♦  SLOPE » 

1  SR*  Cl  . A1 INN?  *A2 < NN ) .A3  I NN> , Aa ! NN ) 

97C  FORMAT!'  SHORE  “  POSITION  COORDINATES  -  X  V  •'.F7.1. 

1  4A4  //9X' SHORE  LATITUDE  “'.F6.U.'  ONSHORE  AZJMjTm 

2  ■*  >Ft • 0 • 1  OEwREtS* /9X, • NE ARShORE  SLOPE  a'»Ft,3.'  AvEHAOL  F£K« 

3 

:r  :  if  T  l ^  -  -9,‘v.te 
SR  «i*ST*C*F 
•.*tN*CvF 

-KITE.  3. 92?.  Si  •  T  *  »ol  >  N.<.»  ,bZ  >  NN  :  .S^Lw.S^hi 
92*  •  riots  -  SF'RIN-J  TlDt  NANOt  •••►'S.2*'  NEAP  TIDt  RANu£ 

ic8«*  »ZAa  yS*  i  'Sitk  At  ..  Co  T  }  j£  •  *  •  F  6  •  3  »  1  SlOPE  A’  N|u«  ’  1  Dt 


**  *  .  ,  .  b.'i4l  7  ;  c  t  —  P  oRM  NUMHfc  K  IS*»cfc«2/l 
>  »«iUl  ►  N 

»l  oEv  I  i* !  .  ***«Al.  Tic.  -  f>v  Nl^EER  lS'»Fs»Z/> 
Avu  -  N.Ll  .  -  «0  I  ’  E  «'*••**•  A  :  *■  N 

.  .  IuhV.„  '■  I  ->t  “  FORM  NUMBER  Is*  » F  b  •  2  /  I 

w  ••*  ' .  t  -  number  is1 * 


i'lUKl, 


!!Jv 
*  -  1 


i  '  .  :  -» 

A  .  *  ?  * 

If  :  !►  ■ 
F JR" A' 


3* -  . »RA' 

: .  ‘uLo 
*  .TP  t 

; »  » i  *• 

-Zi  F  v'^* 
S.S  a  «  I  T  L 
F.Rmj, 


K  •  o 

a :  ■.> 


‘  *i  1  i  :  .  ,  •  1  >  ’  •  !  ;  “i  ‘r  1  ,  1  *1  •  *  !  >  j  ‘  v;  1  »T  J  '  i  ’hAt-.,  1  , 

'  t_'  *  '  i-u*  F  .  '  .  T  r  V  «  '  -NSh  '  .  Tt-6  •  •  AuSh  • .  T  f«.»'fcMtCli'*Tel* 

'  A  v  t  '  i  T  vt  •  '  3'<£AhE  h  ■  «  T  1  .  v  .  *  uSL  '  •  .  •.A.Tjc'PRtvo.' 

'  »  *  S  V  ,  1  a  I  Ni.  '  •  !  it  I  1  »  1  N  j  1  •  I  ?  3  ♦  '  *  1  NO  *  .  :  b  i  .  1  h  '  •  •  V  j  . 

*A\  •  ,  r;.t  .  ’  '  • 

. a -  I T  £ . 3.R6:  : 

*li  . 1 T I. E 1  I 


ts.'KAL'  .  T  >1  •  '  UAL  1  *  T  3  8  .  *  ^  U  '  ,  ’  a  t>  , 
.•M/S'*T1a.'M/S'  .  t  8  3  •  *M‘  •  1  t>  V  *  *  St  w  * 
Pt  L  '  ' 


T  ;  V  ,  ’  NM  '  •  T  .  .  1  '  i  ?  <<.  i  '  RAD  '  I  T  3  1  ♦  '  MAL  *  •  T  36  .  M  NCnt  S  '  •  1  -  - 

•k ! S'  .  160.  'Fli.SIb’t'llyM’H.'  NTS'  .:e2.,FT'•tbV•'sE'-, 
TK\  i'uI.  j'  •’  iv  3i*F  ’  SE  2  '  ) 

.0 '  ..RIU  3.r22 

*  1  i  -  ,  '  F  »  '  i 


L  ,;*lI  SWBhGuTINLS  T.  i.'E  TERMINI  w  uC  A  T  I  ON  AND  COMPUTE  KIND’  w  A  v  t  S  • 

C  AND  TICES  AT  EA'.m  k.uCA!Iv/N 


DO  60  I «1 »N< 

<Ak.L  L!.'.*!  ;  J.  ;  '  »V  i  I  I  tVi.OC  ibMAi  •  *  it*.'  j  ; 

Ca.l  ;  t.  I  p  S  •  A  *  b  ,  i  •  •  j,PMlN.t-MA>,‘-l,f«AD.AA,»A,D*Ai 

C*Lw  .  PM  I  N  .PMAji  ,  Sl  *  !  .  a ^RF  .  X  L  SM  ,  ON  S  H  •  >  A  N  a  .  t  H  A  ^  ,  P  i  .  A  .  *  A  . 

r  A  .  SMAZ  »E  *Z  »C  *A  «  <  Sv*-S  »f  4  a  NO  ' 

IF ,  I.f-.i-  JL  T.  u: 

LAWW  de:a*  i  T  ;  NT  .►'tRCD  .HE  J  I 
H-f.c»-HE  ;.,T 

t  -  :  A»E  <  Ef  aNS  .-rt  ,’uT  «  T  M  *DRA  Tf»  J 

■jL  * 

C  R  A  *  ?»  .  T  IN’ 

E  a-<v-C  .0 

t  N  •  V  •  0 


I  HOUR"  I  STB  I 

6Zw  CAul  ‘ E* Cm' A,B .* • ' .STf CHt 

f l f c h ■ s  rr:~ 

[  *  BNk  0  •«,*.  T  *STF(.H'TtFLHaBNFC.H 

CALw  wAvES 'F y wND.DRATN. »LFCH,tFtTC,Hf IDT, PER -D.il 
I ►  F  I  *  *  I D  >  St  , N6  »  66 

66  :»VL  ?i:i  .  S’ »  TN.FN.  I  S  TRT  »  Th«  ,  T^A*-,  T  ID*  .S^Pm;  ,SlP»-v,  Si.uPE  #  I  .  T  MEAN  > 
‘f  (*ti  NuRF  •  S»*ANv*P<  R  .0  »hE  IOT  ,  ScCPf  •  Hf  •  HR  AN  j  »  Vl  <  .DlP’Hn  SLJ» 

;e  pFAN.-l^O.’  4»9.fc>,*»B 
A  A  -3*.  AN0»  L1  P  AN  --  J  fe  -  • 

Ui  ^ON'iNu- 


'  .R.  ,T(  aAvt  AN.  jSjbH'.HI  tM*  j’ 

EuS».  26**  -.»•••>  •*  2.  *Dl  PT 

if  .  -LS  A*  .4*  .a  * 

4  6  t  N • i N  •  ‘  .  S 

:,  *  a  a 

A  a  S 

A  A  .  Nut 


.1  f  Bu; 


ESUM«ESUM-*E 

OPT  *01 1 J*A»CkN-..S 

VPT«V  (  I  )*A*C<N..  3 

wSRF-i*SuRF*CkTS 

ON»ONSm*OCTS 

Al«AlSh*C*TS 

EF«£FnND«CKTS 
hT 1 *nfc 1 GT*CHF 
MB1»H8*CMF 
TID1*T IDX*CMF 
Vl«VLS*CCMF 

PP-P/CPQ 

*«1TE<  3.510)  1  HOUR  .OPT.  VP  T  . X , Y ,PP , SHANG»**RF »0N *Al . EF  .HU.PtKjj, 

1  HB1.8RANG.V1 

5  30  f  OKMAT  I  3X.  12. 2F  7. 0*2F  7. 2. F  8. 2. «*F  7. 1. F  7.  l.Ftt.Z.F  7,  X. FT. Z*>  /.i,F  7.2  1 
IFt IFTIO.GT.OI  wRITt: 3.53})  T ID1 
531  FOKMATIIm+.IINX.FS..:  ; 

1 HQU R »  1  mOUR  *  I  »■  1*  (  T  !NT  J 

1 F  <  ThOwR.G7«24U^OUR-1hOuR-:4 
50  CONTINUE 
fc.ijCT  »EN*EP 
EGOM»ESoM*CoGol 
LwSCT-tLSCT *C JOUu 
EP»fc  p*C JOUL 
EN-EN^CjOUL 

%  ,  iniLV,3.',lhV  “u"'!-i'NN'.Cil»».).C3l».S..£tiCT,Clcl«l.l.CilN«l.ti1NM 

5>*J  FORMATtlX./'/.lA.TS'MAVE  ENERGY  IN  THE  BREAKER  ZONE  -  1  » E 1 0 . 3  .  3A4  /  / 

* 1 X » T 5 . *  TOT AL  LONG-SHORE  CURRENT  Energy  « * » E 1 0 • 3 » 3A4/ > 

WRJTE13»542I  EP«Ci i NNi *C2lNN» #C3 INN) »EN»Cl INM *C21NNJ »C3»NN) 

5«*2  FORMAT  \  lx  « 15  .*  TOT  Al  POSITIVE  lONG-ShQRE  CURRENT  tNERGv  »'*E10«3* 
13aa//ix»t* total  negative  long-shore  Current  energy  •  * .tio. 3 ,3an ; 
GC  TO  ttOBO 
1000  CALL  EXIT 
ENO 

.  /  OjP 

•DELETE  STrm* 

•store  m S  UA  S’RMa 

n  JOB  FORE  3-Bb-001 

/  /  FOR 

•  1 OCS » C ARt  # i 1 32  PRINTER) 

•ONE  WORD  INTEGERS 
•LIST  SOURCE  PROGRAv 

subroutine  FORtiNx.A.NAwT i 

COMMON  U ( 1301 .v; 130) 

c 

C  FORECASTING  -  COMPUTE  STORM  POSITIONS  FROM  INITIAL  POSITION. 

C  VELOCITY  AND  AZIMUTH 


TINT. 1 .0 
RAD. 57. 2958 

READ l 2  •  90 3 )  SVF L»AZI .g( 1) *V( 

903  FORMAT ;«F 7.0) 

IFlNAuT)  l.i.2 
1  WR I T  £ (  3.9101  SVEL.AZI 

910  FORMA M «X STORM  VELOCITY  «'.F4.0» 
l  9X  » ' sTORm  AZIMUTH  .»fF4.0/l 

GO  TO  5 

\  WRITE i 3.911  )  SvEl.AZI 

911  FORMA! I9X. ' STORM  VELOCITY  ■•.F4.Q. 
1  9X. 'STORM  AZJMgTh  «*.F4,C/) 

5  CONT INUE 
Ui  1  )  *UU  »/A 
Vlll*vll)/A 
01  ST  *SVEL«T INT/A 
AZM.90.-AZ  I 
IF(AZM)  31*32*32 

31  AZM«AZM*360. 

32  DO  33  I *2 • NX 

U!  I  » .U! 1-1  I *DIST.COS l AZM/RAD  I 
Vi  I  ) *V I  I-l)*DIST#slN(AZM/RAD) 

35  CONTINUE 
RETURN 
E  NO 

/  /  DuP 

•DELETE  FORE 

•STORE  «S  v*  FOkt 


KILOMETERS /HOUR’ 


•IOCS! CARD* 1132  PRINTER) 

•ONE  WORD  INTEGERS 
•L  1  S  f  SOuP'.t  PROGRAM 

Subroutine  minoina.ai 
common  Ui 110  »  *v'  130» 
u  I  mf  n  S  I  ■-  N  *13  0'**'  30) 
c 

h:*.;;*ST;nG  -  COMPUTE  STURM  POSITIONS  at  1  HOUR  INTERVALS  *NjM 
c  fc  wuuk  POSITIONS  ON  WEATHER  maps 


DO  15  l  *  1  •  N  x 

RE  AO  I  2 .902  >  A  1 1  :  *T<  ; 

902  FORmat i 2 F  7 . 0 i 
x 1 l . ■ X : I)  /  A 
v  i  ;  \  m  1 ! ; I /A 
I  5  C  .-«r  {NGE 
NX  i  .NX-1 
D‘-  2  -  I  •  j  .  N  x  . 

I  I  •  I  •  *  “  5 


V(I|)>YI!) 

DlF(j-  (  X  (  !♦!  >-X  (1)1/6* 

0 1 FV- ( Y ( 1*1 )-Yi 1 )  )/6. 

DO  20  J-l.S 
M  •  I  l  ♦  J 

U(M)«g(M-l i+DIFu 
V(M)-V(M-I )*01FV 
20  CONI  1 NUt 
I  I »NX«6-S 
UM1)>XINX) 

VI  1  1  )-Y (NX  > 

NX- 11 

RETURN 

END 

//  DUP 

•DELETE  m 1  NO 

•STORE  WS  UA  HIND 

//  JOB  LOCAT  3-66'OOi 

//  FOR 

•10CSICARD.H32  PRINTER) 

•ONE  WORD  INTEGERS 
•LIST  SOURCE  PROGRAM 

SUBROUTINE  LOCAT I  US T *  VS T *ULOC » VLOC *SHA2  .x.r.EAZ.P.J) 

C 

C  U  AND  V  COORDINATES  APE  READ  INTO  thE  PROGRAM  IN  UlOMETERS  AND 

C  CONVERTED  10  UNITS  PROPORTIONAL  TO  Th£  MAjOK  Axis.  ThE  u  AND  V 

C  COORDINATE  SYSTEM  IS  A  RlC TanowLAR  GRID  wITm  U  PGInTINo  LAST, 

c  v  pointing  north,  and  the  origin  located  to  the  suuTh„lst  of  the 
C  SHORE  LOCATION. 

C 

C  SUBROUTINE  LOCAT  IS  USED  TO  ESTABLISH  thE  X  AND  y  COORDINATE 

C  SYSTEM  W I T  M  ThE  OR  I G  l  N  AT  ThL  CEnTL**  oF  Th  t  STORM  AND  POSITIVt 

C  Y  P.j  I  \  T  I  No  .NShOKE,  Nt  GAT  I  v't  v  POINTING  OFFSHORE.  PColTIvE  a 

C  TO  TmL  RIoh*  FAl:*,j  aj-vKt.  AND  NEoAllfE  X  TO  ThE  LEFT  F  AC  i  N 

C  . NSnjH; . 

RAD-S  • 
o*oi T -L_CL 
v-lCl-vSY 
*SoR  «  l  u»»2  *  v  **2  i 
CAcL  A  n  C  T  A  ANui.iV ) 

A«A*,o-LHA.’  /  RA. 

A*-.'  -lCS  «  A  ) 

y  •  :  •  S  I  N  I  A  i 

CI'7*EA2*90,-ShAZ 

IF  .0I»  .CT.36D.  I  DlF-Dn-360. 

A£  *A»DI F /RAD 
P— 2«CCSIA Z  ) 

0  •  Z  •  G 1  N  l  A  i  ) 

RE  * ORN 
END 

/  '  DUP 

•DkLETE  LOCAT 

•STORE  WS  JA  LOCAf 

ft  JOB  T 

//  xtu  CSTRM  2 

•L OCAlCSTRM, lOCAT. WI ND.DECAY.ET I  Mt.FLTCH.wAvt  S»SU*F  . L NR u » . T I Dt ,m 1 ND . 

•localforec 

HOLLAND.  MICHIGAN 
1  JULY  7,  1970  11 

101  2A  ISO.  1.0  *2. 

1000.  0.  MO.  .033 

JOB  ElJPS  j-66-k/Ll 

/ 1  FOR 

•  I OC 5 l CARD  *  1 1 32  PRINTER) 

•ONE  WORD  |NTFGt»S 
•LIST  SOURCE  PROGRAM 

SuBROUT INE  EL IPS( A.B .xl . > 1 ,pmi s ,pma a .»1 .ERAD.XA  .YA.D2A  ■ 

C 

C  SUBROUTINE  EL  IPS  IS  vSEu  T  -  D  E  *  i  A  **  I  N  *  ThI  « I  NO  ANolE  AND  PRLSSw-L 

c  GRADIEN7  at  any  Point  w]7~,N  an  f.LIPMCAC  ;»Torm. 

c 

PAD-57. «.93e 
R-9/A 

AA-A/ A 
HA-b/A 

*  AND  B  *  me.  The  m*j>  An'/  m  « NOR  A*ES  OF  ThE  STORM  ELLlPfct 
C  POINT  Xl.*l  LIES  _-N  A  SlL  ^N.  ELLIPSE  «1th  AXES  A1  ND  01  . 

Al*S^OTcxl«»2*-»**t'/R*w2i 
R  1  •  L  •  A  1 

0-4  1  -H  i  ••  2  -  A1--2- *  1 

C 

c  <.  ;s  1  -«t  INTERLCTION  The  a  Alls  -ITh  A  line  normal  to  The 

C  "ANjEN*  OF  'he  SECOND  Ellipse  AT  >1,11, 


Y2*s wor.eA«»2«<;. -**■••  2 -  aa««2i 


ARE  A  ThE  IN’lAStC'  I  ON  OF  ilM  NoRMAl 

*  'Hi  SEC-‘‘.D  III  ;ri!  <  A  ,c  S'CiMM  ti^IPSt. 


Dl«SJ«n  M1-a0i**2  *  ^1**2  i 
D2«S*JRT  (  I  *2-*C  :  ••i  *  r  2»*2  ) 

; F  i  D  2  A  ,  jl . A  )  Ji  A  «  A 
iKAD*i;i/'J2 
*aC 

?0  I f «AM i n-PMl N 

P  INC-P'.  lF«fcApl-2i**2/2.  I 

Pi-F'MX-PINC 

Pi  r  S  ThE  aAkOMt'KIC  AT  A  0  COMPoTEO  oN  A  noKmAl  Cu«vt 

AloNG  ThE  ha.,,.*  AaIS,  ThIS  PRESSURE  vAlUE  IS  USED  IN  DE T  ERM  I  N I  NG 
ThE  PRtiSURC  o*<  At  I  ■  N  J  NORMAL  K  ThE  l SOBAR  AT  P01NT  A  i  *  <  ,  . 

’.*(  • 1.E-.0.3:  '  i 

:-u:.tu.a.ji  * i « .  .  - . 


CDV-'jTt  T  Hj.  T  AN  jtNT  7  .  Th*t  Ei-IPSE  AT  XI  AND  »i  T  ^  DETER^lNt 

T.it  *  I  No  r  .xh.  >  .  N  . 

IF  i  a  i  ,  j 7  # 0  •  •  A \D  •  >  .  « o  1  •  0 •  C  )  oO  T-;  11 

IF  ;  X  1  .  G’  .0.0  .  AM?.  Y  i  ,L  T  .0*0  I  -V  T0  12 

If  •  ••  .  ..  '.0.  J.ANi*.  r  :  .GT  .0.0  !  O'.'  To  li 

IF  •  .  .l  T,3,  J.nD.Y.  .i-T.O.i?)  oo  T0  11 

*.  A  -  -  .*  A 


20  CwNT INUE 

•  ;  T 

END 

/  Duf 

DELETE  fcLJPS 

S  T  C.  ft  E  aS  UA  ELlPS 


!  DC  SI  CARD.  U32  PRINTED) 

one  **ord  integers 
cIS*  S  >UMCL  PROGRAM 

So'JPs.ljT  I  Nt  WIN.  I  A.Y  .P.FMlN.PMAx.b^AT  . 
lPl.A,<A.YA.b«A2.tA2.D2A.I  SuND.EF>*ND  ) 


OftF ,ALSh*ONSM. SmANG.EHAL  « 


SdBR-.TINE  aIN„  .  J  oit-1  T  0  DETERMINE  ThE  GEGSloPnlC  #  1 ND  SPEED 
AN.  DIRECTION  «  .**  ThE  oARCMt  7  R  |  c  PRESSukE  GRADIENT.  T  hE  SURFACE 
.PfcE"  AN.  .  IKEv'IDN  ARE  Dt  T  t*M INEG  OVER  LAND  AND  SEA  FkCv 
.N«VL  „  *  IDnS  AhI^ied  10  The  of  OS  7  RoPh  I C  *IND»  ThE  ONShoRE  AZJ^-jTh 
I  j  uSuJ  r .)  F  I \J  (i'l  .  NSr  ><•<:  AN,.  A _ONGSK oft£  COMPONENTS  OF  7 Mt  *IND# 
rHt  L^t-TI/t  *.  [  V-  cto  is  v/StD  t  -<  F  l  NO  Thu  WAVE  HEIGHT,  PERIOD. 
AN  j  l*  N  >S-io-E  C  y  -  •  :  N  '  v  :  ^  C  I  T  -  . 


T-s  d  A 

.  t ,  Oi.  A  ’  E  ThE 


_,.nE  4HA.JH.NT  A_0No  A  NORMAL  CoRvE  IS  USED 

TRoF'-;.  ft  IN’.  SPEED. 


C OR  «  2 • # oME G A • D  J Nl A  T / NAD  • 

JO. •ERA D 
OIF  »S>WAX-  ■  i 
»  INC3*.  I  f  •  t  a  p  <  -Z  •  *  2  /  2  •  ' 
i  .  JSC 

: \  =  ;*pinc/o2a 

*.  D-nE  •  ’  =■  '  S  I  -L  ;R  )  *  t  LiP'0N/S00  0.  I 

Awl  A>*CTA'ANj»*A.VA. 

DI  AZ^ftD.-S^AZ 
I  F  :  „  If  .  -;T  ,  *  ->0  ■  )  DSF3DIF-360. 

*  an:,.  -  *  ad«ang-d:  f 

IF  „  ANG.L  r  ,0»r-  1  ftANO«wANG*36C* 

I‘  .  i*N  j.ljT  ,3  c,  D«  :  ft  A  N  .  «  W  A  N  j  -  ?  t>  0  . 

"  AND  SC*A  A-*L  ^•R*tC  *  I  .'N  FAl.  ’jFS  u Sc D  DETERMINE  ACt  »IN 

SP-  f.D  AND  .  IRE  .TI  N  ft'O*'  GE-. STH  -PhU  wiN.  uvtR  i. AND  AND  SLA 

e  .  AN;  -  .  .  I:  .  .  AN.  y‘*A  -  *V 

*  c  ‘  -  P  •  *  .. .  >  r>T;  AN,-  P  E  T  A  -  NO 

..  ■  -N>'.-i  -  -  :  jI  '  I  .E  Y 

»G.  1  A.  N ,4-  w!  -  -  .  :i  T  ;  j E  * 

■  ■■...  *  :f  -  ,->t  -  NF  j*  T  1  .,  t  y 
.  '  .  -  a  .  -  ne  ,a’  :  vt  y 


P’IjN.  V  and  t)E  '  A  vAlUES  *■  -R  The  Sla  A  R  t  wStD 


M*  .  30<,v  *S 

'  A  «  V.  . 

"  so 


A  .  SIN  TOASiF,  fc-ir,  ,S  iS  ^  St  o  N 

s  !  au-  s:-i  1  * ** i  s*4^ - i .  : ne  . 

. .%  1  ,  '  S  I  .  ‘  ?  • 1  ’  x  A  N  s  1  T  :  -  ■  N  V  •  • 


30  OEGREE  ran*,; 
.AND  AND  SEA  L.- 
aM  ,e  .  A  r  ]  ON  A  *  I  .  t  ^ 


AD-A087  858  WILLIAMS  COLL  VILLIAMSTOWN  MASS 
COASTAL  STORM  MODEL. (U) 

APR  76  IT  FOX»  R  A  DAVIS 
UNCLASSIFIED  TR-14 


F/«  4/2 

N00014-69-C-015& 


c 

29  SI NS-S I N  <  SANG/ RAO ) 
iF(SiN$-.5>  35.30  00 

30  B*. 00019 
BETA-29. 

GO  TO  50 

35  IF<51NS*.5)  36* 40.90 

36  B-. 000065 
BETA-50. 

GO  TO  10 

40  JF(SANG-90.)  41*41.42 

41  SANA-SANG*3. 

GO  TO  48 

42  I F ( SANG-1 80 . )  43*43*44 

43  SANA-180. -3. *(180. -SANG) 

GO  TO  46 

44  IF ( SANG-27Q. >  45*43*46 

45  SANA-180. +3. *(SANG-180. I 
GO  TO  46 

4b  SANA-360.-3.*(360.-SANG> 

46  B».O0Ql*(1.275+Q.625*SIN( SANA/ RAO  > I 
BETA- 39.5-10. 5»S INI  SANA /RAO) 

50  6ETR-8ETA/«A0 

C0TA-SIN«BETR)/C0S1BETR)4C0A/(B*C0S(BETR) > 

ALPHA-ATAN1 l./COTA l*RAD 
I F ( NX )  51.52.52 

51  SHANG-WANG490. -ALPHA 
SANG-SANG-ALPHA 

NX-i  * 

GO  TO  27 

52  CONTINUE 

VHC0R-C0R-5IN( ALPHA/RADI/ » 8-COS  1 BETRI  ) 

WSURF -VHCOR*nSGEO 

c 

C  ONSHORE  AND  LONGSHORE  COMPONENTS  OF  ThE  SURFACE  WIND 

C 

ONSH-USURF *COS I SH ANG/RAO ) 

ALSH-nSURF*S1N(SHANG/RAD) 

c 

c  EFWND  IS  The  EFFECTIVE  HIND  SPEED  USED  IN  DETERMINING  HAVES 

C  AN  OFFSHORE  HIND  IS  ASSUMED  TO  BE  *333  TIMES  AS  EFFECTIVE  IN 

C  GENERATING  HAVES  AS  AN  ONSHORE  HIND. 

C 

EFHND-H^jRF* l .6667* • 3333* COS  I SHANG/RAD ) > 

IF(SHANG)  55.60.60 
55  SMANG-SHANG*360. 

60  CONTINUE 

IF ( SHANG-36G. J  70*70.65 
65  SMANG-ShANG— 360. 

70  CONTINUE 
RETURN 
END 

//  DuP 

•delete  wind 

•STORE  HS  UA  WIND 

//  • 

//  JOB  DECAY 

//  FOR 

• l OCS I  CARD *1132  PRINTER  > 

•ONE  WORD  INTEGERS 
•LIST  SOURCE  PROGRAM 

SUBROUTINE  DECAY  I T INT *PEROD . HE IG T» 

C 

C  SUBROUTINE  DECAY  IS  USED  TO  FIND  THE  DECAY  IN  HAVE 

C  THE  HAVES  MOVE  AWAY  FROM  THE  STORM  CENTER. 

C 

c 

C...  FIND  LOGARITHMIC  ATTENUATION  COEFFICIENT 
C 

FREO-l./PEROD 

ATCOF-10.**( I (FREQ-O. 061/(0. 0324>)-l.) 

IF(ATCOF.GT.l.O)  ATCOF-l.O 
IF(ATCOF.LT.O.l)  ATCOF-O.l 
C 

C...  FIND  PROPAGATION  DISTANCE  IN  DEGREES 

C 

01  ST- 1 . 560T»PERQ0*T  I NT* 360. / 400 74. 

C 

C...  FIND  DECAYED  HAVE  HEIGHT 
C 

HE1GT.HEIGT<M2.71BSJ**C-2.4(0.11I14ATCOF)»DIST ) 

RETURN 
END 

//  DUP 

•DELETE  DECAY 

•STORE  WS  UA  DECAY 

//JOB  ETIME  3-86-001 

//  FOR 

*  I OCS ( CARO  *  1 1 32  PRINTER) 

•ONE  WORO  INTEGERS 
•LIST  SOURCE  PROGRAM 

SUBROuT INE  ETIME ( WSNEW.OLDH* HOURS* DRAT N) 

C 

c 

C  SUBROUTINE  ETIME  IS  USED  TO  DETERMINE  THE  EFFECTIVE  FETCH  ANO 
C  DURATION  FOR  ThE  NEW  HIND  SPEED  FROM  ThE  LAST  PREVIOUS  WlNO  SPEED* 

C 

c...  convert  m  to  ft  and  m/sec  to  ft/sec 


3-86-001 


HEIGHT  AS 


90 


c 

M-0lDH/0.»04* 

WS-WSNEW/Q.3040 

c 

Cm«  FlHO  EFFECTIVE  FETCH  USING  NEW  WIND  SPEED  AND  OLD  WAVE  HEIGHT 
C 

U-( ( J2.1«H»/(O.203*WS**2. ) ) 

SUM-U 

DO  10  K-l.51.2 
I-K+2 

SUM-SUM+U**I/l 
10  CONTINUE 

F.<w***2./32.i>*<SUM/0.012*>«*2.38 

C 

C...  CONVERT  FROM  FT  TO  NM  AND  FT /SEC  TO  KTS 
C 

F-F/6Q74. 

WSKTS»WS«3600./4074. 

C 

C...  FIND  EFFECTIVE  DURATION  USING  D6RIVE0  EFFECTIVE  FETCH 
C 

EFDU»-(F/t tWSATS/10*  J •*0.72*10***0«3) )«*0.8 

DRATN-MOURS+EFDUR 

RETURN 

END 

//  DUP 

•DELETE  cTIME 

•STORE  W$  UA  ETIME 

//  JOB  FETCH  3-04** 00 1 

//  FOR 

•IOCS <  CARO *1132  PRINTER ) 

•ONE  WORO  INTEGERS 
•LIST  SOURCE  PROGRAM 

SUBROUTINE  FETCH (A*B*XA«YA*FI 
C 

C  SUBROUTINE  FETCH  IS  USED  TO  FIND  ThE  EFFECTIVE  FETCH  AS  A  STORM 

C  MOVES  ONSHORE  OR  OFFSHORE.  ThE  FETCH  LENGTH  IS  DETERMINED  FROM 

C  THE  FETCH  AREA  WHERE  THE  WIND  IS  BLOWING  TOWARD  THE  SHORE  SITE. 

C  THE  EFFECTIVE  F£TCh  AREA  IS  AN  ELLIPSE  WITHIN  WHICH  THE  WIND  IS 

C  BLOWING  AT  GREATER  THAN  HALF  THE  MAXIMUM  WIND  SPEED.  THE  MAXIMUM 

C  EFFECTIVE  FETCH  IS  THE  MAJOR  AXIS  OF  ThE  ELLIPSE. 

C 

R«MA+B)/2. 

FMAX-0.S8ttl*R 
PI-3.14159 
RAD-57 .2958 
X  *XA*A 
Y*YA«A 

IF  (ABS(X).GE.(0.).AND.  ABS ( X  I . L T . ( 1 • / 3 . *R > )  GO  TO  100 
IF  (ABSfXI ,G£. (1./3.-RI.AND.  ABS I X) *L 7. 1 0.444444«R) j  GO  TO  200 
IF  ( ABSIX) ,GE. t0.444444«R) )  GO  TO  300 
GO  TO  400 
C 

c...  x  greater  than  or  equal  to  zero  and  less  than  1/3  r 

c 

100  Yl-SURT ( « l./3.*R)**2.+<l./2.*FMAX>**2.~X**2«> 

Y2-SORT ( ( l./3.*R>»*2.-X*«2.) 

YJ-O.O 
Y4--Y2 

Y5»-SURT(R*«2.-X*«2. I 
1FIY.GE.YU  GO  TO  110 
IFl  Y.LT.Y1.AND.Y.GE.Y2I  GO  TO  120 
lrt Y.LT.Y2.AND.Y.GE.Y3)  GO  TO  130 
IFl Y.LT.Y3.AN0.Y.GE.Y4)  GO  TO  140 
IF<Y.LT.Y4.AND.Y.GE.Y5)  GO  TO  150 
IF  i  Y.LT.YS J  GO  TO  ISO 
GC  TO  4C0 
110  F-FMAX 

GO  TO  500 

120  RATIO-(Y-Y2»/ <Y1-Y2 > 

CuRVl»l.-l 1 .-COS l P I *RAT 1 0 1 >/2. 

F-1,/2.*FMAX»CURVE«1./2.»FMAX 
GO  TO  *00 

130  RATIO- ABStXI /SORT <  A«*2.*Y2*«2. I 
FMIN-RATIO-l./2.*FMAX 
RANGE -l./2.«FMAX-FM|N 

CURVE-1 *-COS( U.-I  V2-VI /(Y2-V3I »*PI/2.» 

F-CURvE»RANGE-FMIN 

GO  TO  500 

14U  RATI G- A0S (XI/ SON  T l X**2.-Y4*-2.  > 

c 

C...  NOTE...  This  IS  The  SAME  ratio  as  130  SINCE  Y 2  •  Y4 

c 

FHIN.RAT 10-1./2.-FMAX 
*»ANGE-1./2.*FMA»-FHIN 

curve- i i.*ccs<  t y-ya > /  <  y*-v4i  mi n/a* 

F -CURV£»RANGE*FM1N 

OJ  TO  too 

150  *<AT  I  J-  I  Y-Y5  )  /  i  T--Y3I 

CURVE- l  l.-COSI  »  l  .-RAT  IQlMIM/2. 

F-CwRvt«U/2.«FMAX 
*0  TO  *00 
140  F . o • 0 

GO  TO  *00 

c 

C...  X  GREATER  THAN  OR  EOuAi.  TO  1/3  R  AND  Lt*S  TMA 
C 

200  VI- SORT ( I 1./2.4FMAI )4»2««< 1 • / I. -A > **2 .-X*«2 • » 

Y2  — 0GRT(R**I.«X*»2.  f 


R 


wr 


IF  I Y.GT #Vl I  GO  TO  210 
!FtY.LE.Yl.AND.Y.GT.Y2)  GO  TO  220 
IF<V*U.Y2»  GO  TO  230 
GO  TO  400 
210  F-FMAX 

GO  TO  *00 

220  RATIO-lVl-Y)/IYl-Y2> 

CURVE-U.*C0SIRAT10*P1  )  )/2. 

F»FMAX*CURVE 
GO  TO  *00 
230  F-0.0 

GO  TO  *00 
C 

C...  X  GREATER  THAN  OR  EOUAL  TO  0*444444  R 

C 

300  Vl-0.0 

Y2«-SORTIR**2.-(0.4444*R)*»2. ) 

IFIY.GT.YI)  GO  TO  310 
IF  I Y .LE.Y1 .AND. Y.GT «Y2 )  GO  TO  320 
IFtY.LE.Y2)  GO  TO  330 
GO  TO  400 
310  F-FMAX 

GO  TO  500 

320  RAT10*IY1-Y)/|Y1-Y2 ) 

Curve* < i •♦cos ( rat io*pi i >/2. 

F»FMAX*CURVE 
GO  TO  *00 
330  F*0*0 

GO  TO  500 
C 

C...  ERROR  CONDITION 
C 

400  WR|TE(3»410) 

410  FORMAT ( IX • 1  LOGIC  ERROR  ENCOUNTERED  IN  SUBROUTINE  FETCH...*) 

F*0.0 

RETURN 

*00  !F(X.LT.(0*n  F  "0.5*F 

RETURN 
END 

//  DUP 

•OELETE  FETCH 

•STORE  WS  UA  FETCH 

//  JOB  WAVES  3-86-001 

//  FOR 

•IOCS < CARO  *  1 132  PRINTER) 

•ONE  VORO  INTEGERS 
•LIST  SOURCE  PROGRAM 

subroutine  waves(speed*time.fetch.efetc»heigt.perod»iout i 
c 

C  SUBROUTINE  WAVES  IS  USED  TO  DETERMINE  WAVE  PERIOD  AND  HEIGHT  FROM 
c  wind  speed,  duration  and  fetch. 
c 

C...  UNIT  CONVERSION  PACKAGE . .  .KM  TO  NM...-M/SEC  TO  KNOTS 

C 

FCHNM-0.5400*FETCH 
SPDKT ■ 1 . 94394SPEED 
C 

C...  EFFECTIVE  FETCH  PACKAGE 
C 

EFETC* t SPDKT / 10. )••©. 72*10. •*0.S*TIME**1.2* 

IF IFCHNM.LT .EFETC )  EFETC*FCHNM 
C 

C...  UNIT  CONVERSION  PACKAGE. ..NM  TO  METERS. . .KNOTS  TO  M/S EC 

C 

EFETC*EFETC/0. 5400*1000. 

SPDMS-SPDKT/1.9439 

C 

C...  WAVE  PARAMETER  PACKAGE 

C 

CONST-(9.i062*EFETC)/SPDMS**2. 

HEIGT«ISP0MS**2.*0.2B3*TANHI 0.012**CONST**0.42) )/ 9.8062 
PEROD«t2.*3.141*9*SPDMS*i.20*TANHl0.077*CONST**0.25) )/9.8062 
EFETC “EFETC/1000. 

IF i IOUT.EU.1 )  RETURN 
IF(10uT.EU.2)  GO  TO  200 
C 

C...  ERROR  CONDITION 

C 

WRITE  <  3.100 1 

100  FORMAT (IX* > ILLEGAL  OUTPUT  OPTION  IN  SUBROUTINE  WAVES*) 

RETURN 

C 

C...  NON-METRIC  OUTPUT  OPTION 
C 

200  EFETC-EFETC*0.*400 
HEIGT*HEIGT*|.2aO» 

RETURN 

CNO 

//  OuP 

•DELETE  WAVES 

•STORE  WS  UA  WAVES 


//  JOB 

//  FOR 

•  IOCS I  CABO til 92  PRINTER) 
•ONE  WORD  INTEGERS 
•LIST  SOURCE  PROGRAM 


92 


SUBROUTINE  TIOE  (ST.TN.FN.ISTRT .ThR.TDAY.TIDX  »SL  TOP.SLBQT .SLOPE »  i  • 
1  TMEAN) 

c 

C  PERIODS  FOR  THE  MAJOR  DIURNAL  AMO  SEMIDIURNAL  TIDAL  COMPONENT* 

c 

TM2  **12.42 
TS2-12.00 
m«2S*9| 

TOl-25.t2 

HR*24.»TDAY+FLQAT < ISTRT KFLOAT 4|-1 J-THR 
ARG*6.2S32*HR 
C 

C  COMPUTE  THE  AMPLITUDE  OF  THE  TIDAL  COMPONENTS  FROM  THE  SPRING 

C  AND  NEAP  TIDE  RANGES  AND  THE  TIDAL  FORM  NUMBER. 

C 

FM2-($T+TN)/(4.*(1.+FN)  ) 

S2-<ST-TNI/U.*(1.*FN>  I 

FK1-FN*FM2 

Ol«FN*S2 

c 

C  TIDE  LEVEL  IS  THE  SUM  OF  THE  FOUR  MAJOR  TIDAL  COMPONENTS 
C  AT  EACH  HOUR. 

C 

TIDX  •  FM2#COS <  ARG/TM2 )  *  S2*COSlARG/TS2 FKl«COS I ARG/TK1 )  ♦ 

X  01*COS(ARG/T0l ) 

C 

C  BOTTOM  SLOPE  IS  COMPUTED  FROM  THE  TIDE  LEVEL.  AND  THE  SLOPE 
C  AT  HIGH  AND  LOW  TIDES. 

C 

TLOC  ■  (TIDX+ST/2.J  /  ST 

slope  ■  slbot  ♦  tlqc  •  isltop-slbot> 

T 1DX-T IDX+TMEAN 

RETURN 

END 

//  DUP 

•DELETE  TIDE 

•STORE  WS  UA  TIDE 

//  JOB  SURF  3-S6-001 

//  FOR 

•I OCS I  CARD  *1132  PRINTER) 

•ONE  WORD  INTEGERS 
•LIST  SOURCE  PROGRAM 

^  SUBROUT INE  SURF ( SHANG.PEROD.HE IGT .SLOPE. HB.BRANG .VLS.DEPTH .LSCOP ) 

C  COMPUTE  BREAKER  HEIGHT.  ANGLE  AND  LONGSHORE  CURRENT  VELOCITY 

PI0.141S9 
RAO-57. 2956 
G-980.42 
HTCM-100.*MEIGT 

hBCM“ •  38ft*G**0. 2* l PEROD*HT  CM«#2 ) ••0. A 
DEPTH»1.2*HBCM 
WLO-G/(2.*PI >  *PEROD**2 
WL1-WL0 
C 

C  DISPERSION  EQUATION  USED  TO  DETERMINE  SHALLOW  WATER  WAVE  LENGTH 

C 

DO  10  1*1*20 

WLi-WL0«TANH(2.*PI*DEPTH/WLl 1 
10  CONTINUE 
C 

C  SNELL ' S  LAW  USED  TO  DETERMINE  BREAKER  ANGLE  AND  HEIGHT 
C  DUE  TO  WAVE  REFRACTION. 

C 

$INA«SIN(SHANG/RAD)*TANH12.*PI *0EPTM/WL1> 

COSA-SQRTI l.-SINA»*2) 

CALL  ARCTA (8ANG.C0SA.SJNA ) 

rrang-rad*bang 

REFRC-SQRT ( COS (SHANG/RAD) /COS (BANG) I 

HS»RCFRC»H8CM 

C 

C  OPTION  OF  FOUR  DIFFERENT  LONGSHORE  CURRENT  EQUATIONS. 

C 

GO  TO  ( 1.2.J.4) .LSCOP 

C 

C  LONGSHORE  CURHEHl  -  FOX  AND  DAVIS.  1972 

C 

1  *LS>100.'»$LOPE*HB/PCftOD*S!N(4.*BRANG/RAD) 

GO  TO  5 

C 

C  LONGSHORE  CURRENT  -  LONGuET-H I GG I NS.  1970 
C 

2  VLS  •  9.0#SLO®E*SQRT I  G*Hi  ) *S IN ( 2.*BRANG/RAD ) 

GO  TO  ft 

C 

C  LONGSHORE  CURRENT  -  C.E.R.C..  197ft 

C 

ft  VLB  ■20.7*SLOPE*SQAT (  G«H0  I *SIN< 2. *SRANG/RAO > 

GO  TO  ft 
C 

C  LONGSHORE  CURRENT  -  KOMAR  ANO  INMAN.  I 9 TO 
C 

4  Gl-G/100. 

VL I* 100  «*Gl *ME IGT  ••2/i,»Gi*PERQD/<4,*PI )#SIN<2.#»RANG/RAP>* 

1  COS ( 2 .••RANG /RAO ) 
ft  HB-hB/100. 

RETURN 

INC 
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//  DUP 


•DELETE  SURF 

•STORE  WS  UA  SURF 


//  JOB 
//  FOR 

•IOCS (CARD *1132  PRINTER) 

•ONE  WORD  INTEGERS 
•LIST  SOURCE  PROGRAM 

SUBROUTINE  enrgy ih»t •time* i salt .E. IUNIT  ) 


E  NR  GY 


3-86-001 


c 

c 

c 

c 

c 

c 

Cat* 

c 


SUBROUTINE  energy  IS  USED  to  determine  WAV!  energy  during  EACH 
HOUR.  THIS  ENERGY  IS  DERIVED  FROM  THE  ENERGY  WITHIN  A  SINGLE  WAVE 
AND  SUMED  TO  FI  NO  THE  TOTAL  ENERGY  IN  THE  STORMa 


FIND  NUMBER  OF  WAVES  IN  GIVEN  TIME  PERIOD 


WAVES-T  IME*3400, /T 
C 

c...  FIND  PROPER  MASS  DENSITY  OF  WATER  (0-FRESH . 1-SALT I 

IF ( ISALT aEUaO )  DENSE-1.94 
1F( ISALT.EO.l)  DENSE-2.0 
IFt ISALT aGTal >  GO  TO  SOO 
C 

Caaa  CONVERT  WAVE  HEIGHT  FROM  METERS  TO  FEET 

HT-H/O.3Q40 

C 

Caaa  FIND  ENERGY  FOR  A  SINGLE  WAVE 
C 

^  E-S.12*DENSE«92al*tHT«T)**2a/8a 

Caaa  FIND  TOTAL  ENERGY 
C 

E«E*WAVES 

C 

Caaa  IF  OUTPUT  IS  DESIRED  IN  FT-LBS/FT*RETURN 
C 

I F ( I UN I T a EO a 0 )  RETURN 
IF ( JUNITaNEal )  GO  TO  300 
C 

Caaa  CONVERT  FT-LBS/FT  TO  JOULES/Mf TER 

E»E*1. 35382/0. 3048 
RETURN 
C 

Caaa  ERROR  CONDITION 
C 

300  WRITE ( 3*600 ) 

600  FORMAT < IX. 'ILLEGAL  OPTION  IN  SUBROUTINE  ENRGY* t 
RETURN 
END 

//  OUP 

•OELETE  ENRGY 

•STORE  WS  UA  ENRGY 


//  JOB  ARCTA  3-86-001 

//  FOR 

•ONE  WORD  INTEGERS 
•LIST  SOURCE  PROGRAM 

SUBROUTINE  ARCTAI ANGLE #XaY) 

C 

C  ARCTANGENT  ROUTINE  IS  USED  TO  DETERMINE  THE  ANGLE  FROM  X  AND  Y 
C  COORDINATES#  OR  SIN  AND  COSINa 
C 

ANGLE  ■  OaO 
IFlABS(X)-. 001)2. 9»9 
9  IF ( X  )  1.3*3 

3  IFIY)  4,9.3 

4  ANGLE  -  4.283 1825 
GO  TO  S 

1  ANGLE  "  la  1413926 

9  ANGLE  •  ANGLE ♦AT AN ( Y/X l 
RETURN 

2  IF  (ABS(Y)-aOOl)  8.10.10 
10  IFl Y I6#7#7 

6  ANGLE-4, 7123889 
RETURN 

7  ANGLE  •  1.9707963 

8  RETURN 
CNO 

//  OUP 

•DELETE  ARCTA 

•STORE  WS  UA  ARCTA 
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//  JOB 
//  FOR 

•!OCS< CARD* 1132  PR  INTER .TYPEWRITER .KEYBOARD) 

•ONE  WORD  INTEGERS 

COMMON  U( 1 30 ) • V  U  30 ) 

DIMENSION  TITLEI20)  .0AY15) 

DIMENSION  All2).A2(2).A3(2>.A4l2)*Blt2)  .8212)  »CH2)  .C2I2I.C3UI 
DATA  Al/*  KIL*.'  NAU« / *A2/ 'OMET ' . ' TICA' / *A3/ *ERS  *.»L  Ml*/. 

I  AA/*  *.*LCS  ‘/.Bl/*  MET*.*  Fe£*/,a2/»ERS  **'7  •/. 

2C1/*  JOU*.*  FOO'/*C2/'LES  '.'T  PO*/.C3/'  ••'UNOS*/ 

STORM  PROGRAM  LIMITED  TO  FORECAST  MODE 
TYPEWRITER  INPUT  FOR  THE  STORM  PARAMETERS 

CARD  1  -  TITLE 

CARD  2  -  STARTING  DATE.  TIME  AND  INPUT-OUTPUT  OPTIONS 

COLS.  1-2  IS7RT-  STARTING  HOUR 
COLS.  3-22  DAY  -  STARTING  DATE 

COL.  23  INAUT-  INPUT  OPTION 

0  -  METRIC  UNITS 
1  -  NAUTICAL  MILES  AND  FEET 

COL.  24  NAUT  -  OUTPUT  OPTION 

1  ■  NAUTICAL  MILES.  FEET  AND  KNOTS 
Q  -  METRIC  UNITS 

CARD  3  “  STORM  PARAMETERS 

COL.  2  TIDE  PREDICTION  OPTION 

0  -  TIDE  PREDICTION  NOT  INCLUDED  -  OMIT  CARD  4 
1  -  TIDE  PREDICTION  INCLUDEO  -  SEE  CARD  4 

COL.  3  LONGSHORE  CURRENT  EUUATION  OPTION 

1  «-  FOX  AND  DAVIS.  1972 

2  ■  LONGUET-HIGGINS.  1970 

3  ■  C.E.RaC. .  1973 

4  -  KOMAR  AND  INMAN.  1970 

COLS.  4-6  NX  -  NUMBER  OF  STORM  POSITIONS 

SIX  HOUR  INTERVALS  FOR  HINDCASTING 
ONE  HOUR  INTERVALS  FOR  FORECASTING 

COLS.  7-12  BNFCH-  AVERAGE  BASIN  FETCH  IN  KM.  i NAUT •  Ml.) 
COLS.  13-17  TINT  -  TIME  INTERVAL  BETWEEN  STORM  POSITIONS 
NORMAL  SETTING  IS  1.0  HOURS 

COLS.  32-36  SLAT  -  LATITUDE  AT  SHORE  SITE 

CARD  4  -  TIDE  PREDICTION  -  OPTION  FROM  CARD  3 

COLS.  0-5  ST  -  SPRING  TIDE  RANGE  IN  FEET 
COLS.  6-10  TN  -  NEAP  TIDE  RANGE  IN  FEET 

COLS.  11-15  TDAY  -  NUMBER  OF  DAYS  FROM  LAST  SPRING  TIDE 

COLS.  16—2 0  ThR  -  HOUR  OF  LAST  HIGH  SPRING  TIDE 

COLS.  21-25  FN  -  TIDAL  FORM  NUMBER 

0.0  TO  .25  -  SEMIDIURNAL  TIDE 
.25  70  1.5  -  MIXED  SEMIDIURNAL  TIDE 
1.5  TO  3.0  -  MIXED  DIURNAL  TlOE 
GREATER  THAN  3.0  -  DIURNAL  TIDE 

COLS.  26-32  SLPLO-  NEARSHORE  BOTTOM  SLOPt  AT  LOW  TIDE 
COLS.  33-39  SLPHl-  NEARSHORE  BOTTOM  SLOPE  AT  HIGH  TlOE 
COLS.  40-46  TMEAN-  MEAN  UOE  LEVEL  IN  FEET 

CARD  5  -  SHORE  SITE  LOCATION 

COLS.  1-7  ULOC  -  X-COORDINATE  IN  NAUTICAL  MILES 

COLS.  B- 14  VLOC  -  Y-COORDINATE  IN  NAUTICAL  MILES 

COLS.  15-21  SHA2  -  ONSHORE  DIRECTION  -  CLOCKWISE  FROM  NORTH 
COLS.  22-2  0  SLOPE-  AVERAGE  NEARSHORE  BOMUM  SLUPC 
COL.  30  ISLNO-  0  -  CONTINENTAL  COAST  OR  BARRIER  ISLAND 

The  FOLLOWING  INFORMATION  IS  TYPED  IN  FROM  THE  CONSOLE  FOR  EACH 
STORM  SIMULATION  RUN • 

PM|N  -  MINIMUM  BAROMETRIC  PRESSURE  IN  MILLIBARS 
pMAXR—  PRESSURE  at  LARGEST  ENCIRCLING  ISOBAR 

Type  ar  -  major  half-axis  (effective  long  radius) 

SR  -  MINOR  HALF-AXIS  IEFFECTICE  SHORT  RADIUS) 
IN  EA2  -  ORIENTATION  OF  MAJOR  HALF  AXIS  PLUS  OR 

MINUS  90  DEGREES  FROM  NORTm. 

SVEL  -  STORM  VELOCITY  IN  KNOTS 
All  -  STORM  A  2  I  MU  T  h  -  CLOCKWISE  FRQM  NORTH 
XU)  -  INITIAL  X-COORDINATE  FOR  THt  STORM. 

YU)  -  INITIAL  Y-COOROINATE  FOR  Th£  STORM, 

CNK* 1 • S33 19 
CKN*. 53961 
CFM-.304B 
CMFO.2BO0 
CRTS* 1 .9425 
CCMF-.0J2B0B 
CJOUL  ■  • T JT J61 
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c 

RAD-97. 2958 

8080  HEAD  12.915)  TITLE 

919  FORMAT / 20A4 i 

READ (2*917)  I STRT .DAY » INAUT ,NAUl 

917  rORMATU2.9A4.211> 

NN-NAUT-1 

IF < INAUT)  1.1.2 

1  CNR* 1.0 
CFM-1.0 

2  IRiNAUT)  11.11.12 

11  CRN- 1.0 
CMF-1.0 
CRTS* 1.0 
CCMF-1.0 
C JOWL- 1.0 

12  IP  I ISTRT. 60.0.0)  00  TO  1000 

READi 2.901 » 1N0PT .  I  FT  J D .LSCOP .NX .BNFC 1 *T 1  NT .SLAT 
901  FORMAT  191 1.I3.F6.0.F 5. 1 .14X.F 5.1) 

8NFC«-aNFCl«CNK 
IFUFTiO)  5.9.4 

4  READ ( 2 . 920 )  SI  .T l »TOAY . TmR.FN.SLPlO.SLPhI .TM£A1 

920  FORMAT I5F5.2.JF7.2 ) 

ST-$1*CFM 

TN-T1«CFM 

TMCAN.ThEA1.CFM 

9  READ  I  2  «  907  >  ULC1  .VLC1 »ShAZ .SLOPE . ISLND 
907  FORMAT  1 4F7.0 . 1 2 ) 

ULC1»ULC1*CNA 
VLC1-VLC1»CNK 
949  tfRXTEl 1 .947) 

947  FORMAT  (  •  TYPE  (N  ThE  M I  NUMljM  AND  MAXIMUM  BAROMETRIC  PRESSURE  IN  Ml 
ILL  1  BARS  -  FORMAT  AAAA.BBB8.M 

READi 6 .948 )  PMIN.PMAXR 

948  FORMAT (2F5.0) 

«R 1 TE ( 1  .959  ) 

955  FORMAT  1 •  TYPE  IN  ThE  LENGTH  OF  THE  MAJOR  AND  MINOR  HALF  AXES  OF  TH 
IE  STORM  -  FORMAT  XXXX.YVYY.M 
READ  <6*994)  AR1.BR1 
996  FORMAT (2F5.0) 

WRITE  I  1.960) 

960  FORMAT ( *  TYPE  IN  ThE  ORIENTATION  OF  ThE  MAJOR  HALF  AXIS  IN  DEGREES 
1  FROM  NORTH  -  FORMAT  XXX, '» 

READ  1 6 .961  >  EA Z 

961  FORMAT ( F9.0 ) 

C 

C  BAROMETRIC  PRESSURE  AT  LARGEST  ENCIRCLING  ISOBAR  IS  ASSUMED  TO  BE 

C  AT  TWO  STANOARO  DEVIATIONS  FROM  ThE  CENTER  OF  THE  STORM.  TO 

C  FIND  THE  ACTUAL  storm  SIZE.  THE  MAJOR  AND  MINOR  AxES  ARE 

C  MULTIPLIED  By  1,5  AND  ThE  MAXIMUM  PRESSURE  IS  MULTIPLIED  BY 

C  1.145  TO  DETERMINE  TmE  PRESSURE  AT  ThE  MARGIN  OF  ThE  STORM. 

C 

AR-AR 1 *CNA 
BR-0R1-CNA 
rangp-pmaxr-pmin 
PMAX-RANGP-1 . 145.PMIN 

c 

c  Alt  ANO  BR  ARE  RADII  OF  THE  STORM  OUT  TO  TwO  STANDARD  DEV  I  AT  IONS. . . 

C  MULTIPLY  BY  1.9  TO  YIELD  ThE  FULL  LENGThS. 

C 

A-1.5*AP 

B-1.5*BR 

CALL  FQREC INX.A.SVEL.AZ1 t 
DO  42  1-1 .NX 
Ul 1  ) -U ( I  ) *  CNX 
42  VI  I l-VI 1 ) *CNA 
UL0C-ULC1/A 
VL0C-VLC1 /a 
Jl-Ui 1  )  •A..005 
Vl-V  I  1  l-A.,005 
UL -ULOC- A 
VL*VLOC«A 
C 

C  WRITE  SHORE  AND  TIDE  DATA 

C 

UL1-ULPCAN-.05 
VLi*VL*CXN..05 
•NFCl-tNFCH*CXN 
WRITE* 9.916)  TITLE 

918  FORMAT  I 1H1//.1K.2UA4/ ) 

WR| TCI  9.9181 1STRT .DAY 

918  FORMAT ( IX • 'RUN  BEGINS  AT  MOgR  '*12. 1  ON  '.9A4) 

WRITE  I  9.905  i  PMIN.PMAXR, PMAX 

909  FORMAT ( 1M0.' STORM  -  BAROMETRIC  PRESSURE  AT  CENTER  OF  LOW  • '  ,F  7. 1 . 
1'  MILLIBARS* /9* . 'PRESSURE  AT  LARGEST  ENCIRCLING  ISOBAR  -  *  .F  7 , 1 » 

2 '  MILL IBARS* /9x • 'MAXIMUM  PRESSURE  INCLJDEO  IN  STORM  -'.F7.1. 

I'  MILLIBARS* /) 

AR1- ARwCKN+.QS 
BR1-BRWCRN..D9 

WRITE  I  9.928  I  AR1 *A1 INN) *A2 INN ) • A1 1 NN ) * A4 I NN ) *BR 1  *  A 1  I NN ) . A2 I NN ) . 

1  A9IMN) .A4INN) .EAZ 

928  FORMAT  1 9X. 'LENGTH  OF  MAJOR  HALF  AX18  ■'  .F7.1.4A4/ 
l •* • ' LENGTH  OF  MINOR  half  A*|S  •*. FT. 1.444/ 

29X.'0R|ENTATI0N  of  major  axis  - ' »F 7. 1 , *  DEGREES  FROM  NORTh'/I 
IP  I NAyT 1801 *801 .802 
801  wRtTEl 9.9101  SvCL.AZI 

910  FORMAT  I 9X. 'STORM  VELOCITY  -*.F4.0.'  A  I  LOME  TER S/*OUR ' • 

1  9X, ’STORM  AZIMUTH  • '  .F 4 . 0 / i 

GO  TO  809 

602  WRITti 9,911  *  SVlu.AZI 

911  FORMAT i **,• STORM  VELOCITY  •♦.F4.0.'  ANOTS'. 

I  9X  i*  STORM  AUmuTh  -'.Rw.l/f 

809  CON  7  J  Nuf 
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WRITCll.9?Q>ULl.VLl.Al!Nh» » A  2  (  NN  ) »  A  3 { NN ) i  A4 ( NN ) .SLAT • SHAZ .SLOPE * 
l  BNFC1  .AKNN1  .A2INNI  .  A3  ( NM  ,  A4  (NN ) 

•70  FORMAT ( 1  SHORE  -  POSITION  COORDINATES  -  X  •'.F7.1.*  Y  -*.F7.1. 

1  4A4  / /9X 1  SHORE  LATITUDE  -'.F6.0.'  ONSHORE  AZIMUTH 

2  • '  .F6.0 » 1  DEGREES' /9X. 'NEARSHORE  SLOPE  ■••F6.3.'  AVERAGE  FETCH 

3  -• .FA.0.4A4/J 
IF(IFTIO)  59.59.58 

58  Si*ST*CMF 
Tl-TN*CMF 

WRITE  <  3.927  I  S1.T1  .BUNN)  »B2(NN)  .SLPLO.SLPHl 
927  FORMAT  I '  TIDES  -  SPRING  TIDE  RANGE  •••F6.2»*  NEAP  TIDE  RANGE  ■'* 
1F8.2.2A4  /9X. 'SLOPE  AT  LOW  TIDE  »*.F6.3.*  SLOPE  AT  HIGH  TIDE 

1  ■ ' »F6  •  3/ ) 

IFIFN.LE..25)  WR1TE<3.941)  FN 

941  F0RMATI9X.' SEMIDIURNAL  TIDE  -  FORM  NUMBER  IS1. F6.2/) 

IFIFN.GT. • 2 5 • AND .FN • L£ • 1.5 )  WRITE(3.942t  FN 

942  F0RMATI9X. ‘MIXED  SEMIDIURNAL  TIDE  -  FORM  NUMBER  IS'.F6-2/> 
IFIFN.GT-l .5 • AND.FN. LE • 3.0)  WRITEI3.943J  FN 

943  FORMAT (9X» ‘MIXED  OIURNAL  TIOE  -  FORM  NUMBER  IS’»F6.2/J 
IFIFN.GT.3.0)  WRITE! 3 .944  I  FN 

944  FORMAT (9X. ‘DIURNAL  TIDE  -  FORM  NUMBER  lS'.F6.2/> 

59  WRITEU.830J 

830  FORMAT ( '  PLEASE  CHECK  TO  SEE  THAT  The  INFORMATION  ON  The  PRINTER  I 
IS  CORRECT*  /  •  IF  THE  PROGRAM  IS  ALL  SET  TO  GO.  PRESS  THE  1  XEv .  0 
2THERWISE#  PRES S  THE  2  XE Y‘> 

READ  16*831 )  NNN 

831  FORMAT  (ID 
IFINNN-1)  590.590.945 

590  WRITE0.500) 

500  FORMATI1MO.*  HOUR'  .T 1 1  . ' X •  .T 1 8  •  * Y ' *T24 • • XI ' . T 3 1 • » Y 1 • . T 3 7 • • BARO. ' . 
1T45. ‘WIND’ *T52  » ' SURF • • . T59 . • ONSM • » T66 * ' AL$H ' . T 72 . ' EFFLCT . • . T 8 1 » 

2 'WAVE*  .T  88.  'WAVE*  .T98  *  '  BREAKER  '  •  T109.  •  LSC  *  .  /  .  IX  .  T  36  .  '  P  RESS  .  '  .Tww. 
3' ANGLE*  .T52.'WIND' »T59.' WIND' .T66.' WIND* .T73.‘ WIND  1 . T  83  » ' H 1  .T  9  J  « 

4 ' T ' .T97.‘H‘.T 101 .'ANGLE' • T 108 . • VELOC. • ) 

IF < 1FTID.GT.0J  wRI TE l 3.951 ) 

951  FORMAT ( 1H*  »T1 17. ' T IDE < > 

IF(NAUT)  800.800.809 

800  WR1 T£ ( 3.801 1 

801  FORMAT (lHO.T 10. 'AM* *T17.'XM' »T24.'RAD’ .T 31 . *RAD » » T 38 • ’ MB • *T45. 

1  'DEG* .75 3.’ M/S'  *  T  60 .'M/S' .T 4 7 . *  M/S • . T 74 • ‘M/S • #T83»  »M‘ .T 89. 'SEC* 

2  • T97. *M# . T I  02 • 'OEG* .T108. 'CM /SEC' > 

IF  l  IFT  ID.GT.O)  WRITE!  3.802) 

802  FORMAT !lH*.Tli8»'M‘ ) 

GO  TO  815 

809  WRITE! 3.820  I 

820  FORMAT!  1H0.T  10.' NM*  .T17.  *NM‘  .  T24  ,  •  RAD  •  .  T31 .  '  RAD'  .  T  38  .  ‘  MB  •  .  T  «.*->  . 

1  *D£G‘ #753. 'KTS' . T fcO . • XTS • • Tfc7 , *  K7S* ,T 74 • • XTS ' • T82 • • FT • • T89 » * SlL * 

2  .T96. 'FT' .T 102. ’ DEG • .T1O0.' FT/SEC ' » 

IF  t IFTID.GT.O)  WRITEI  3.822J 

822  FORMAT  ( 1H-*  .T110.‘FT‘) 

815  WRITE ! 3.952  > 

952  FORMAT (10X) 

C 

C  CALL  SUBROUTINES  TO  DETERMINE  LOCATION  AND  COMPUTE  WIND.  WAVES. 

C  AND  TIDES  AT  EACH  LOCATION 

C 

00  50  I-l.NX 

CALL  LOCAT  <U ( I I.VIl  I .ULOC.VLOC. SHAZ .X.Y.EAZ.P.O) 

CALL  ELIPSIA.B.P.U.PMIN.PMAX  «P1 .ERAD.XA.TA.D2A) 

CALL  WINO<X.Y.P.PMIN,P*AX.5LAT.WSURP.AL5H.ONSH,SMANG.CRAD#P1 .A,XA, 
1YA.SHAZ.EAZ.02A* ISLND.EFwND) 

IF ! I. EO.lt  GO  TO  510 

CALL  DECAY (TINT .PEROD.Mt I GT ) 

HDECV-MEIGT 

CALL  ET IME IEFWND.HE IGT .TINT .DRATN1 
GO  TO  520 
510  ORATN-TINT 
E SUM-0.0 
EN-0.0 
EP-0.0 
IHOUR-ISTRT 

520  CALL  FETCh!A.B.X.Y.STFChJ 
TLFCh-STFCH 

IF!BNFCM.lT.STFCM)TlFCH»BNFCH 

CALL  WAVES ICFWND.DRATN.TLFCH.EFETC .m£ IGT .PEROD.l I 
ir I  IF T 10 1  56.56.55 

55  CALL  TIDC(ST.TN.FN.ISTRT.THR.TDAY.T|DX.SLPHI . SLPLO. SLOPE • I .TMEANI 

56  CALL  SURF (SHANG.PEROD. HE IGT. SLOPE. MB. BRANG.VLS. DEPTH. LSCOPI 
IF(BRANG-I«0. »  49.49.48 

48  BRANG-BRANG-360. 

49  CONTINUE 

C  COMPUTE  WAVE  ANO  LONGSHORE  CURRENT  ENERGY 

lLB»-O257»(0.6*VLS»*«2.*0CPTHW«2./SLOPE 
IF ( VLS )  45.46.47 
48  EN«EN«ElS 

GO  TO  66 
47  cp«ep*els 

46  CONTINUE 

CALL  CNRGYtHB.PfAOO.TINT.l.C.l I 

ISUM«ESUN‘E 
UPT-U! 1 ) *A*CXN*. 5 
VRT  «V ( 1 14A«CRN*.5 
WSRF-WSURF •CXti 
ON»ONSM*CKT$ 

AL*ALSM*CXTS 
EF-EFWN0WCXT4 
MT 1 «H€ 1 GT  *CMF 
Mtl«HB*CMF 

Vl-VL8*CCMF 
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i*HBl*aRANG°v/H0UR  *UPT 'VPT  *X*Y,P*SHANG,WSRF,ON,At*EF*MTl»PEROO» 

330  FWjJJ‘3X»l2*2p7.0.2F7*2.F«.1.4»F?,lfF7.1.F6.2»F7.X.F7.2*2F7.1l 

IF! IFTID.GT.O)  WR l TE <  3. 53 1 )  TlOX 
531  FORMAT ( 1H+ .114X.F5.2) 

IHOUR-IHOURMFJXITINTI 
I F II HOUR . GT • 2 4 >1 HOU R ■ I HOUR-2 A 
30  CONTINUE 

elsct.en^ep 

ESUM-ESUM*CJOUL 

elsct-elsct*cjoul 

EP-EP*CJOUL 

EN-EN*CJOUL 

WRITE  13 *540)  ESUM.C1 (NN) .C2(NN) .C31NN) .ELSCT.CUNN) »C2 (NN ) »C3 ( NN) 
540  FORMAT ( 1X»//»1X»T5'WAV£  ENERGY  IN  ThE  BREAKER  ZONE  •  '  »E 1 0 . 3 . 3A4 , ✓/ 
1.11X.T5. ' TOTAL  LONG-SHORE  CURRENT  ENERGY  • • E10 • 3 » 3A4 » / ) 

WRITE! 3*542  I  EP.Cl INN) .C2 (NN) ,C3(NN> .EN.Cl(NN) *C2(NN) »C3(NN) 

542  FORMAT ( IX . T3  * 1  TOTAL  POSITIVE  LONG-SHORE  CURRENT  ENERGY  ■••E10.3. 

1 3A4 »//*lX»T5** TOTAL  NEGATIVE  LONGSHORE  CURRENT  ENERGY  -'.E10.3. 
23A4 ) 

WRITE! 3*992) 

992  FORMAT  UM1) 

WRITE! 1 *990) 

990  FORMAT ( *  IF  YOU  WANT  TO  TRY  ANOTHER  SIMULATION  RUN.  TYPE  IN  A  1. 
IOTmERwISE  TYPE  IN  A  ZERO') 

READ l 6 .991 )  NEND 

991  FORMAT (II) 

!F<NEND.eo.l >  GO  TO  943 
GO  TO  8060 

1000  CALL  EXIT 
END 

//  DUP 

•OELETE  CSTRM 

•  STORE  WS  UA  CSTRM 


n  JCB  FOREC  3-66-001 

//  FOR 

•IOCS (CARD. U 32  PRINTER.TYPEwR|TER*KEYBOARD) 

•ONE  wORD  INTEGERS 
•LIST  SOURCE  PROGRAM 
•LIST  SOURCE  PROGRAM 

SUBROUTINE  FORECINX.A.SVEL.A2I ) 

COMMON  U( 130 ) .V( 130) 

C 

C  FORECASTING  -  COMPUTE  STORM  POSITIONS  FROM  INITIAL  POSITION. 

C  VELOCITY  AND  AZIMUTH 

C 

T  I N ).’  ■  1  #  0 
RAD-57. 2950 
WRITE ( 1.973) 

973  FORMAT! •  TYPE  IN  ThE  STORM  VELOCITY  AND  AZIMUTH  -  FORMAT  VVV.AAA. 
1  ) 

READ (6*976)  SVEL.AZl 
976  FORMAT  <  2F4 .0 ) 

WRITE ( 1 .970) 

970  FORMAT! •  TYPE  IN  ThE  INITIAL  X  AND  Y  COORDINATES  OF  THE  STORM  - 
1  FORMAT  XXXX.YYYY. • ) 

READ <6.971)  U I  1 >  . V ! 1 1 

971  FORMAT (2F5.0) 

U(ll«U!l)/A 

V(1  l-Vdl/A 
DIST-SV£L*TINT/A 
AZM-90.-AZI 
IF(AZM)  31.32.32 

31  A2M- A2M+360 • 

32  DO  33  I -2 .NX 

U( I ) -U ( 1-1 )*D1ST«C0S( AZM/RAD) 

V( I ) -V( 1-1 ) ♦OIST*SlN<  AZM/RAD ) 

35  CONTINUE 
RETURN 
END 

/  /  DUP 

•DELETE  FOREC 

•STORE  WS  UA  FOREC 


APPENDIX  B.  HINDCAST  STOW  DATA  AND  OUTPUT 


OAT*  LIST  FOR  HlNOCAST  EXAMPLES 
STEVENSvIlLE.  MICHIGAN 

!SlUl«SS:  ‘I“  ....  “lOW.  .2.  MO.  TOO.  ,0. 

10*3.  A86.  90.  *033 

333  1229 

h1<*  960 

723  TIT 

82*  *«3 

928  6** 

102*  *** 

U7®  678 

1260  7*9 

1*08  913 

1910  1030 

1687  1190 

1*9*  1323 

1773  1*98 

18*3  1*31 

1907  1830 

HOLLAND.  MICHIGAN 

10lUlO  *ii0.7°i»0  1003?  1010.  *2.  330.  200.  30. 

1013  360  90  .033 

961  1030 

106*  9*9 

.098  8*3 

1102  7*1 

11*3  709 

121*  *92 

1271  68* 

1313  683 

1362  686 

1*01  692 

HOLLAND.  MICHIGAN 

101^13  7l»o!7°i*0  1002?  1010.  *2.  *35.  2**.  33. 

1013  3*0  90  .033 

*32  9*9 

388  9*8 

731  913 

88*  8*3 

1027  790 

1102  699 

1131  437 

1151  62* 

1190  603 

12*3  601 

1297  627 

1339  6*7 

138*  6*6 

HOLLAND.  MICHIGAN 

101^12  *  190.191?0  100*?  1012.  *2.  330.  300.  28. 

1013  360  90  .033 

796  703 

832  699 

871  689 

936  686 

101*  666 

1093  630 

1160  6*3 

1196  663 

1300  686 

1371  733 

1**3  816 

15*7  932 

JNCBOrGAN.  WISCONSIN 
7JULT  13,  1972  00 

101  9  200.  1.0  990.  1008.  *2.  1730.  «73.  23. 

1892  1112  270.  .030 

983  179* 

800  1768 

1033  1788 

12*2  1800 

1*69  1639 

1*8*  1976 

179*  2087 

1896  2210 

2028  233* 

SMeBOYGAN,  WISCONSIN 

lin;  X5oo!,T1.0  1002?  1016.  *2.  1030.  700.  20. 

1892  1112  270.  .030 

176  26* 

390  3*3 

398  *10 

793  *9* 

1079  *37 

1387  9*7 

1393  1196 

1833  1*17 

2028  1937 

2080  22*3 

21®1  2313 

Sheboygan.  Wisconsin 

7 AUGUST  1.  1972  00 

101  9  200.  1.0  1002.  1012.  *2.  1*00.  *30.  12. 

1892  1112  270.  .030 

398  1911 

9*3  1203 

1027  1623 

1) 96  1366 

13*3  1340 

1631  16** 

1430  1*00 

2) 61  1983 

2*2*  2203 

mAGOALEN  ISlANOS  -  WEST  StOt 
2NOV1MB16  2*.  197*  00  ,  r 

111  10  200.  1.0  974.  1008.  *8.  86C.  660. 

.86  .37  10.  21.1  3.0  0.026  0,010  .6* 

1663 »  1110.  90.  .026  l 

113*  *46 

1310  348 

1*99  966 

133*  US* 

1587  1221 

1**3  1188 

l  T* I  12’* 

1*3*  133* 

19*3  13*1 

22*2  1920 


MAGDALEN  ISLANDS  -  EAST  SIDE 
2 NOVEMBER  ».  » «•  °°  st,0.  t60.  0.0 

ill  10  200.  1»0  976.  100*. 

!S1  l?37  lo!  21.1  5.0  0.026  0.030  .6* 

1665.  1110*  270.  .026  l 

113*  *66 

1310  588 

1*99  966 

153*  113* 

1587  1221 

16*3  1188 

17*3  127* 

185*  133* 

19*3  15*3 

22*2  1920 

PLUM  ISLAND.  MASSACHUSETTS 
1  APRIL  *•  1975  ini?.  *?,  1200.  600.  15. 

1960  83*  258  .030 

111  139 

610  535 

1091  610 

1258  906 

1720  777 

1850  925 

1961  925 

2072  888 

2090  906 

2202  925 

2350  777 

2*79  666 

2720  592 

3100  550 

3500  500 

CEDAR  ISLAND »  VIRGINIA 

7JULY  l*»  1973  00  .* 

113  13*500.  1.0  1°10*  5®* 

1.72  .90  1*.  18»  ‘60  ‘°35  *039 

2287  1537  285  .037 

195  5*8 

5*2  8** 

932  1H5 

1323  1323 

1738  1312 

1922  1*7* 

2065  1*34 

22*3  1366 

2*19  1399 

2352  1613 

2621  181* 

268*  20*6 

2713  22*9 

SAPElO  ISLAND.  GEORGIA 

1573.  325.  290.  .0155 

619  *08 

93*  385 

1065  373 

1207  373 

1386  *02 

1*85  *20 

1611  *72 

1727  529 

182*  581 

1890  633 

19*2  679 

2000  702 

MUSTANG  ISLAND.  TEXAS 

7JANUA61  12*  1972  00  2S00.  *00.  23. 

Ill  10  1000.  1.0  99C.  1018.  28.  2500.  *uo.  t 

“U  ?«  V.  IB.B  >■>  .•»“  -0i  •** 

2*15.  0.0  280.  *016 

2381  1393 

271*  1535 

298*  .333 

3216  1553 

3501  1*  08 

3  70S  1427 

*0J7  1666 

*023  1742 

*550  1801 

*830  1633 

MUSTANG  ISLAND.  TtXAS 
7JANUAPT  23.  1VT2  UC  ^ 

11)  9  JjOC.  1.0  998.  »  v  i  6  ■  28.  1800.  3j.,  33. 

.83  .*0  25.  14. b  2.5  .C16  .020  .2* 

2*1?.  C.  280.  .Ot« 

l 7 ’ 7  831 

20  3  B  811 

2256  804 

2613  7  v  7 

283*  879 

309*  ee* 

326*  90  3 

3*C2  92’ 

3321  99b 

«GNT£RE7«  CALIFORNIA 
*flBoAB»  13  19b’  Zv 

112  *  3U3‘J  »•  99  8.  1012. 

2.07  .,.9;  6.  18.  l.Oi  » u3b 


\  7.  700.  3*0. 

,0*8  1**9 


8i3.  i  r 
*«5  1233 

1  12  1129 

i lbs  i jc; 

*2* 

;*tc  8  l  ’ 

1*3.  ’10 


«b’6  38* 

«t*C  “  . 

I.Ul.Bl’  1*. 


•  LS’ 
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PI 


HOLLAND*  MICHIGAN 


RUN  BEGINS  AT  HOUR  1  wN  JULY  3*  1970 

STORM  -  BAROMETRIC  PRESSURE  AT  CENTER  OF  LOW  •  1003.0  MILLIBARS 
PRESSURE  AT  LARGEST  ENCIRCLING  ISOBAR  »  1010.0  MILLIBARS 
MAXIMUM  PRESSURE  INCLUDEO  IN  STORM  ■  1011.0  mIllIBARS 

LENGTH  OF  MAJOR  HALF  AXIS  -  550.0  KILOMETERS 

LENGTH  OF  MINOR  half  AXIS  «  200.0  KILOMETERS 

ORIENTATION  OF  MAJOR  AXIS  •  30.0  DEGREES  FROM  NORTH 

LONGSHORE  CuRREM  EUuATION  FROM  FOX  AND  DAVIS.  19  72 

SHORE  -  POSITION  COORDINATES  -  X  -  1O1S.0  Y  •  560.0  KILOMETERS 

SHORE  LATITUDE  •  *2.  ONSHORE  AZIMUTH  -  90.  DEGREES 

NEARSHORE  SLOPE  -  0.033  AVERAGE  FETCH  ■  ISO.  KILOMETERS 


HOUR 

X 

Y 

XI 

Y  1 

BARO. 

KIND 

SURP. 

UNj  1 

ALSh 

EFFtCT. 

WAVE 

WAVE 

BREAKtft 

..SC 

PRESS. 

ANGLE 

M  IND 

WIND 

Ml  NO 

WIND 

T 

H 

ANGLE 

VELOC. 

KM 

KM 

RAO 

RAD 

MB 

DEG 

M/S 

M/S 

M/S 

M/s 

M 

SEC 

M 

DEG 

CM/SU. 

1 

968  • 

1030. 

0.36 

0.05 

1010.98 

210.9 

0.2 

-0.1 

-0*1 

O.0 

0.00 

O.0 

0.00 

—  1  4 . 9 

-0.94 

2 

964. 

1016. 

0.35 

0.03 

1010.95 

211.9 

0.4 

-0.3 

-0.2 

0.1 

0.00 

0.1 

0.00 

-15.7 

-1.03 

3 

1001. 

1003. 

0.53 

0.01 

1020.89 

213.1 

0.7 

-0.6 

-0.3 

0.2 

u.oo 

0.1 

0.00 

-16.5 

-3.24 

4 

1017. 

969. 

0.32 

>0.00 

1010.79 

2  1 4 , 5 

1.2 

-0.9 

-0.6 

J.4 

0.00 

0.2 

0.00 

-17.3 

-5.23 

3 

1033. 

976. 

0.30 

-0.02 

1010.63 

216.3 

1.6 

-1.4 

-1.1 

0.7 

0.01 

0.4 

o.Ol 

-16.2 

-7.85 

4 

1030. 

962. 

0.48 

-0.04 

1013.40 

216.5 

2.6 

-2.1 

-1.6 

1.0 

0  ■  02 

0.5 

0.01 

-19.2 

-11.03 

7 

1066. 

949. 

0.*7 

-0.06 

1010.08 

221.3 

3.6 

-2.7 

-2.4 

1.5 

0.03 

0.7 

0.03 

-20.3 

-14.36 

a 

1071. 

932. 

0.43 

-0.C6 

1009.75 

222.8 

4.4 

-3.2 

-3.0 

1.0 

0.05 

0.9 

0.05 

-21.1 

-17. 77 

9 

1077. 

914. 

0.42 

-0.07 

1009.36 

224,6 

6.2 

-3.7 

-3.6 

2  .< 

0  •  L  7 

1.0 

0.06 

-21.9 

-20.0. 

10 

1092. 

097. 

0*40 

-0.08 

1008.91 

226.9 

3.9 

-4.0 

-4,3 

2.6 

0.09 

1.2 

0 .06 

—  *2.0 

-23. 53 

11 

1087. 

860. 

0.38 

-0.08 

1C08.41 

229.7 

6.4 

-4.1 

-4.9 

2.9 

C.  12 

1.4 

G.ll 

-23*9 

-25.66 

12 

1093. 

862. 

0.36 

-0.09 

1007.87 

233.4 

6.7 

-4.0 

-5,4 

3.1 

0.15 

1.5 

0.12 

-25.2 

—26. 06 

13 

1098. 

843. 

0.34 

-o.io 

1007,31 

238.4 

6.8 

-3.5 

-5,0 

3.3 

0.17 

1 .6 

0.14 

-26.9 

-26.61 

’  4 

1099. 

828  • 

0.32 

-0.10 

1006.64 

242.0 

6.7 

-3.1 

-5.9 

3.4 

0.19 

1.7 

0.15 

-27.9 

-26.12 

1» 

1099. 

810. 

0.30 

•0.10 

1006.37 

246.6 

6.5 

-2.6 

-6.0 

3.5 

0.21 

l.fi 

0.15 

-29.1 

-24. 35 

16 

1100. 

793, 

0.2  6 

-0.10 

1005.92 

232.6 

6.3 

-1.8 

-6,0 

3*6 

0.23 

1.9 

0.14 

-30.  3 

-21*00 

17 

1101. 

776. 

0.26 

-0.10 

1305.49 

260.7 

6.1 

-0.9 

-6.0 

3.7 

0.25 

2.0 

O.U 

-31.3 

-15.26 

IB 

now 

750. 

0.24 

-0.10 

1005.10 

271.4 

6.1 

0.1 

-6.1 

4.1 

0.27 

2.1 

0.05 

-32.0 

-6.27 

19 

1102. 

741. 

0.21 

-0.10 

1004.76 

284,3 

6.1 

1.5 

-5,9 

4.6 

0.31 

2.2 

0.17 

-30.9 

-21.74 

20 

1112. 

736. 

0.21 

-0.11 

1004.70 

306.2 

6.0 

3.5 

-4.0 

5.2 

0.36 

2.4 

0.30 

-25.4 

-41.51 

21 

1122. 

7>C. 

0.20 

-c.13 

1004.72 

330.4 

6.2 

5.4 

-3.0 

5.9 

0.42 

2.6 

0,42 

-15.2 

-47.24 

22 

1132. 

723. 

0.20 

—0*1 4 

2004*03 

350*0 

7.2 

7.1 

-1.2 

7.2 

0.52 

2 . 0 

0.55 

-5.3 

-22*92 

23 

1143. 

720. 

0.19 

-0.13 

1003.02 

3.1 

8.7 

8.7 

0.4 

8.7 

0.67 

3.2 

0.70 

1*6 

0.45 

24 

1133. 

714. 

0.18 

-0.16 

1005*28 

11.7 

10.2 

10.0 

2.0 

10.1 

0.65 

3.6 

0.80 

6.2 

34.13 

l 

1163. 

709. 

0.18 

-0.17 

1003.61 

17.5 

11.5 

11.0 

3.5 

11.4 

1.05 

4.0 

1.08 

9.3 

53.71 

2 

1172. 

706. 

0.17 

-0 . 1 V 

1305.93 

20.8 

12.4 

11.6 

4,4 

12.1 

1.25 

4.4 

1  .27 

10.9 

66.33 

3 

1181. 

703. 

0.17 

-0.20 

1006*27 

23.3 

13.0 

11.9 

5.1 

12.7 

1.43 

4.7 

1.45 

12.2 

76.67 

4 

1189. 

700. 

0.17 

-0.21 

1006.63 

23.3 

13.4 

12*1 

5.7 

13.0 

1.60 

5.0 

1.62 

13.2 

85.07 

5 

1198. 

698. 

O.U 

-0.22 

1007.00 

26.9 

13.5 

12.1 

6,1 

13.0 

1.74 

5.2 

1.75 

14.0 

91.00 

a 

1207. 

693. 

0.16 

-0.23 

1007.36 

28.2 

13.5 

11.8 

6.4 

12.9 

1.86 

5.4 

1.67 

14.6 

97*07 

7 

1216. 

692. 

0.16 

-0.24 

1007.73 

29.4 

13.2 

11.5 

6.4 

12.6 

1.96 

5.5 

1.95 

15.1 

101.04 

8 

1223. 

691. 

0.15 

-0.25 

1008.08 

30.2 

12.7 

10.9 

6.4 

12.1 

2.02 

5.6 

2.01 

15.5 

103.62 

9 

1234. 

690. 

0.15 

-0.26 

1006.41 

30.9 

12.0 

10.3 

6.2 

11.5 

2.06 

5.7 

2.04 

15.8 

105.25 

10 

1243. 

689. 

0.13 

-0.27 

1008.72 

31.5 

11.3 

9.6 

5.9 

;  0  •  7 

1.92 

5.5 

1.90 

16. 1 

102.31 

u 

1233. 

618. 

0.15 

-0.28 

1009.01 

32.0 

10.5 

8.9 

5.*: 

9.9 

1.  71 

5.2 

1.70 

16.3 

97.16 

12 

1262. 

687. 

0.15 

-0.29 

1009.28 

32.3 

9.6 

B.l 

5,1 

9.1 

1.50 

4.9 

1.48 

16.5 

91.36 

13 

1271. 

686. 

0,13 

-0.31 

1009.53 

32.9 

8.7 

7.3 

4.7 

b.2 

1.29 

4.5 

1.28 

16.7 

05.00 

14 

1276. 

686. 

0.13 

-0.31 

1009.70 

33*2 

b.O 

6.7 

4.4 

7.6 

1.14 

4.3 

1.12 

16.8 

7V.  90 

13 

1263. 

68  5. 

0.13 

-0.32 

1009.86 

33.3 

7.3 

6.1 

4.0 

6.9 

0.99 

4 .0 

0.96 

16.9 

74.63 

.8 

1292. 

685. 

0.15 

-0.33 

1010.00 

33.7 

6.7 

5.5 

3.7 

6.3 

0.05 

3.7 

0.04 

17.0 

69.25 

17 

1299. 

684. 

0.13 

-0.34 

1010.13 

34.0 

6.0 

5.0 

3.3 

5.7 

0.72 

3.4 

0.71 

17.1 
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HOLLAND.  MICHIGAN 


Run  begins  at  hour  j  on  july  ?,  iv/o 

STORM  -  BAROMETRIC  PRESSURE  AT  CENTER  OF  LOU  •  1002.0  MILLIBARS 
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LONGSHORE  CURRENT  EQUATION  FROM  FOX  AND  DAVIS.  1972 
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total  negative  lOnu-smcre  Current  energy  »  o.439e  .is  vuulf.s 


r\ 


HOLLAND*  MICHIGAN 


RUN  BEGINS  AT  hOuR  13  ON  JULY*  18.  1970 

STORM  -  BAROMETRIC  PRESSURE  AT  CENTER  OF  LOW  •  1004.0  MILLIBARS 
PRESSURE  AT  LARGEST  ENCIRCLING  ISOBAR  «  1012.0  MILLIBARS 
MAXIMUM  PRESSURE  INCLUDED  IN  STORM  -  1013.1  MILLIBARS 

LENGTH  OF  MAJOR  HALF  AXIS  •  550,0  KILOMETERS 

LE*GTh  wF  MINOR  HALF  AXIS  •  300.0  KILOMETERS 

ORIENTATION  OF  MAJOR  AXIS  «  28.0  DEGREES  FROM  NORTH 

longshore  current  tuuAnoN  *rom  fox  and  davis*  1972 

SHORE  -  POSITION  COORDINATES  -  X  ■  1015.0  V  -  560.0  KILOMETERS 

SmCRE  LATITUDE  •  42.  ONSHORE  AZIMUTH  •  90.  DEGREES 

NEARSHURL  Slope.  ■  0.033  AVERAGE  FETCH  •  150.  KILOMETERS 


HOUR 

X 

r 

XI 

Y1 

BARO. 

WIND 

SURF. 

ONSH 

ALSH 

EFFtCT. 

WAVE 

WAVE 

BREAKER 

PRESS. 

ANGLE 

WIND 

WIND 

WIND 

WIND 

H 

T 

H 

ANGLE 

KM 

KM 

RAD 

RAD 

MB 

DEG 

M/S 

M/S 

M/S 

M/s 

M 

SEC 

M 

DEG 

13 

796. 

705. 

0.17 

0.26 

1011.16 

202.6 

4.7 

-4.4 

-1.8 

1.7 

0.03 

0.7 

0.03 

-11.0 

1A 

805. 

704. 

0.17 

0.25 

1010.95 

202.9 

5.1 

-4.6 

-1.9 

1.8 

0.05 

0.9 

U  *05 

-11.9 

15 

815. 

703. 

0.17 

0.24 

1010.73 

203.3 

5.4 

-4.9 

-2.1 

1.9 

0.06 

1.0 

0.06 

-12.1 

16 

82-. 

702. 

0.17 

0.23 

1010.50 

203.7 

5.7 

-5.2 

-2.3 

2.0 

0.07 

1.1 

0.07 

-12.3 

17 

833. 

7Ul  • 

0.17 

0.22 

1010.26 

204.1 

6.0 

-5.5 

-2.4 

2.1 

0.06 

1.1 

o.ov 

-12.5 

18 

843. 

700. 

0.16 

0.20 

1010.00 

204.5 

6.3 

-5.7 

-2.6 

2.2 

0.09 

1.2 

0.10 

-12.  7 

19 

852. 

699. 

0.16 

0.19 

1009.74 

205.0 

6.5 

-5.9 

-2.7 

2.3 

0.11 

1.3 

0.11 

-12.9 

20 

855. 

697. 

0.16 

0.19 

1009.62 

205.1 

6.6 

-6.0 

-2.0 

2.4 

0.1  J 

1.3 

0.12 

-12.9 

21 

858. 

696. 

0.16 

o.ie 

1009.51 

205.2 

6.7 

-6.1 

-2.8 

2.4 

0.12 

1.4 

O.ie 

-  *  3 . 0 

2? 

861  . 

6  94. 

0.16 

o.ie 

1009.39 

205.3 

6.0 

-6.2 

-2.9 

2*5 

0.13 

1.4 

0.i3 

-13.0 

23 

865. 

692. 

U.  J  6 

0.18 

1009.27 

205.4 

6.9 

-6.2 

-2.9 

2.5 

0.14 

1.5 

0. 14 

-13.0 

24 

868. 

691. 

0.15 

0.17 

1009.15 

205.5 

7.0 

-6.3 

-3.0 

2.5 

0.14 

1.5 

0.  15 

-13.1 

1 

871. 

689. 

0.15 

0.17 

1Q09.02 

205.6 

7.1 

-6.4 

-3.0 

2.6 

0.15 

1.5 

0.  15 

-13.1 

2 

885. 

669. 

0.15 

0.15 

1008.59 

206.7 

7.3 

-6.5 

-3.3 

2.7 

0.15 

1.6 

0*16 

-11.7 

3 

099. 

686. 

0.15 

0.14 

1008.15 

208.0 

7.4 

-6.6 

-3.5 

2.7 

0.16 

1.6 

0.16 

-14.3 

4 

913. 

65  0. 

0.15 

0.12 

1007.71 

209.5 

7.5 

-6.5 

-3.7 

2.0 

0.17 

1*6 

0.17 

-15.0 

5 

928. 

68  '. 

C.15 

0.10 

1007.28 

211.4 

7.4 

-6.3 

-5.9 

2.0 

0.18 

1.7 

0.10 

-.5.9 

6 

942. 

6t  ’  • 

0.15 

0.08 

1006.86 

213.0 

7.3 

-6.0 

-4.0 

2.8 

0. 16 

1.7 

0.10 

-17.0 

7 

956. 

666. 

0.15 

0.07 

1006.46 

216.7 

7.0 

-5.6 

-4.2 

2.0 

O.lB 

1.7 

0.10 

-18.2 

a 

966  • 

663. 

0.14 

0.05 

1006.14 

218.9 

6.7 

-5.2 

“4.2 

2.7 

0.19 

1.7 

0.1b 

-19.0 

i 

975  . 

6?9. 

0.14 

0.04 

1005.64 

221.5 

6.3 

-4.  7 

-4.2 

2.6 

0.19 

1.8 

0.  IB 

-19.9 

10 

985. 

6  76  » 

0.1- 

0.03 

1005.56 

224.9 

5.9 

-4.2 

-*.2 

2.5 

0.18 

1.7 

0.17 

-21.1 

li 

995. 

673. 

0.13 

0.02 

1005.31 

229.1 

5.4 

-3.5 

-4.1 

2.4 

0.  16 

1.7 

0.15 

-22.6 

12 

1004. 

669. 

0.13 

0.01 

1005.08 

234.7 

4.9 

-2.8 

-4.0 

2.3 

0.15 

1*6 

0.13 

-24.5 

13 

1 G  1  A  . 

666. 

0*12 

0.00 

1004.89 

242.2 

4.4 

-2.0 

-3.9 

2.2 

0.14 

1.5 

0.11 

-26.  7 

14 

1028. 

663* 

0.12 

-0.01 

1004.71 

256.9 

3.9 

-0.8 

-3.6 

2.3 

0.14 

1.5 

0.0  7 

-30.3 

15 

1041 . 

661. 

0.12 

-0.03 

1004.59 

2  76.3 

4.1 

0.4 

-4,0 

2.6 

0.15 

1.5 

0.05 

-31.7 

16 

1055. 

658. 

0.11 

-0.04 

1004.54 

296.5 

4.3 

1.9 

-3.8 

3.5 

0.18 

1.  7 

0.  li 

-20.3 

17 

1068. 

655. 

0.11 

-0.06 

1004.55 

319.7 

4.2 

3.2 

-2.7 

3.9 

0.21 

1.0 

0.20 

-20.0 

18 

1081. 

653. 

0.11 

-0.08 

1004.63 

340.5 

4.8 

4.5 

—  1.6 

4.7 

0.27 

2.0 

0.26 

-10.2 

19 

1095. 

650. 

O.io 

-0.09 

J004.78 

355.3 

5.7 

5.7 

-0.4 

5.7 

0.35 

2.3 

o.  36 

-V.4 

20 

1106. 

649. 

0. 10 

-0.11 

1004.95 

3.0 

6.6 

6.6 

0.3 

6. 6 

0*44 

2*0 

1*6 

21 

1117. 

646. 

0.10 

-0.12 

1005.16 

8.7 

7.5 

7.4 

1.1 

7.4 

0.55 

2.9 

0.57 

4.0 

22 

1127. 

6-6. 

0.10 

-0.J3 

1005.39 

12.9 

8.3 

8.1 

1.8 

0.2 

0.67 

3.2 

O.0v 

6.0 

23 

1138. 

645. 

0.10 

-0  •  i- 

1005.66 

16.2 

9.1 

8.7 

2.5 

9.0 

0.79 

3.5 

0.02 

8.5 

24 

1149. 

644. 

0.10 

-0.16 

1005.95 

18. 8 

9.8 

9.3 

3.1 

9.6 

0.92 

3.0 

0.96 

V,  V 

1 

1160. 

643. 

0.10 

-0.17 

1006.27 

20.8 

10.4 

9.7 

3.  7 

10.2 

1.05 

4  .U 

1.07 

10.9 

2 

1166. 

646. 

0.10 

-0.18 

1006.44 

20.9 

10.6 

9.9 

3.7 

10.3 

1.16 

4.2 

1  ■  1  V 

10.9 

3 

1172. 

650. 

0.10 

-0.19 

1006.61 

2C.9 

10.7 

10.0 

3.8 

10.5 

1.27 

4.4 

1.30 

10. V 

4 

1178. 

653. 

0.  1 1 

-0.19 

1006.78 

20.9 

10.9 

10.2 

3.9 

10.6 

1.36 

4.6 

1  .40 

10.9 

5 

1184. 

656. 

O.ll 

-0.20 

1006.95 

20.9 

11.0 

10.3 

3.9 

10.8 

1.45 

4.0 

1.40 

10. v 

6 

1190. 

660. 

0.12 

-0.21 

1007.13 

20.9 

11.1 

10.3 

3.9 

10.6 

1.53 

4.9 

1.56 

10.9 

7 

1196. 

663. 

0.12 

-0.21 

10C7.31 

21.0 

11.2 

10.4 

4.0 

10.9 

1.60 

5.0 

1  .64 

11.0 

a 

1213. 

667. 

0.12 

-0.24 

1007.85 

22.1 

li.3 

10.5 

4,2 

11.0 

1.66 

5.1 

1.70 

11.6 

9 

1231. 

671. 

0.13 

-0.26 

1008.39 

23.0 

11.2 

10.3 

4«- 

10.9 

1.72 

5.2 

1.75 

12.0 

10 

1246. 

674. 

0.13 

-0.28 

1008.93 

23,7 

10.9 

10.0 

4,4 

10.6 

1.76 

6.3 

1.79 

12.3 

11 

1265. 

678. 

0.14 

-0.30 

1009.45 

26.4 

10.4 

9.5 

4.3 

10.1 

1.70 

5.3 

1 . 0  1 

12.0 

12 

1283. 

682. 

0 . 1 4 

-0.32 

1009.94 

24,9 

9.8 

8.9 

4.1 

9.5 

1.64 

5.1 

1,07 

12.  V 

13 

1300. 

686. 

0.15 

-0.34 

1010.40 

25,4 

9.1 

8.2 

3.9 

8.6 

1.47 

4.B 

1  .49 

13.1 

14 

1312. 

694. 

0.16 

-0.35 

1010.67 

25,0 

8.5 

7.7 

3.6 

8.3 

1.35 

4.6 

1.37 

12.9 

15 

1324. 

702. 

0.17 

-0.37 

1010.92 

24.7 

8.0 

7.2 

3.3 

7.7 

1.21 

4.4 

1.23 

u.a 

16 

1335. 

710. 

0.18 

-0.38 

1011.16 

24,5 

7.4 

6.7 

3.1 

7.2 

1.06 

4,1 

1 .00 

12.0 

17 

134?. 

719. 

0.19 

-0.40 

1011.39 

24.2 

6.9 

6.2 

2.8 

6.6 

0.93 

3*0 

0.94 

12.5 

18 

1359. 

72  7. 

0,20 

-C.41 

1011.60 

24,0 

6.3 

5.7 

2.5 

6.1 

0.79 

3.6 

O.ttl 

12.3 

19 

1371. 

735. 

0.21 

-0.43 

1011.79 

23,7 

5.7 

5.2 

2.3 

5.6 

0.67 

3.3 

0,09 

12.2 

20 

1383. 

749. 

0.22 

-0.44 

1011.97 

23,0 

5.2 

4.8 

2.0 

5.0 

0.56 

3.0 

o.5  r 

1  l.v 

21 

1395. 

762. 

0.24 

-0.46 

1012.13 

22.3 

4.6 

4.3 

1.7 

4.6 

0.45 

2 . 7 

0.40 

11.5 

22 

1407. 

776. 

0.26 

-0.47 

1012.27 

21.7 

4.1 

3.8 

1.5 

4.0 

0.37 

2.4 

-  .30 

11.2 

23 

1419. 

789. 

0.27 

•0.48 

1012.40 

21.0 

3.7 

3.4 

1.3 

3.6 

0.30 

2.2 

0.31 

10.0 

24 

1431. 

803. 

0.29 

-0.50 

1012.52 

20.5 

3.2 

3.0 

1.1 

3.1 

0.23 

1.9 

0,24 

10.5 

1 

1443. 

816. 

0.31 

-0.51 

1012.63 

19.9 

2.8 

2.6 

0.9 

2.7 

0. 18 

1.7 

0.  19 

10.2 

2 

1460  . 

839. 

0.33 

-0.53 

1012.76 

18.8 

2.2 

2.1 

0,7 

2.2 

0.12 

1.4 

0.12 

9.0 

3 

1478. 

661. 

0.36 

-0.56 

1012.67 

17.8 

1.7 

1.6 

0.5 

1.7 

0.C7 

1.1 

0.00 

9.0 

4 

1495. 

884. 

0.39 

-0.58 

1012.95 

16.9 

1.3 

1.2 

0.3 

1.2 

0.0  4 

0.0 

0.04 

8.5 

5 

1512. 

907. 

0.42 

-0.60 

1013.02 

16.0 

0.9 

0.9 

0.2 

0.9 

0.02 

0  *6 

0,02 

6.0 

6 

1530. 

929. 

0.44 

-0.62 

1013.07 

15.1 

0.6 

0.6 

0.1 

0.6 

0.01 

0.4 

0,01 

7.6 

7 

1547. 

952. 

0.47 

-0.64 

1013.10 

14,3 

0.4 

0.4 

0.1 

o.<- 

0.00 

0.3 

0,00 

7.1 

WAVE  ENERGY  IN  THE  BREAKER  ZONE  •  0.138E  10  JOULES 
TOTAL  LONG-SHORE  CURRENT  ENERGY  •  0.965E  09  JOULES 
TOTAL  POSITIVE  LONG-ShORE  CURRENT  ENERGY  -  0.962E  09  JOULES 
TOTAi  NEGATIVE  LONG-SHORE  CURRENT  ENERGY  •  0.263E  07  JOULES 


10? 


fii 


lSC 

VElOC, 

CM/StC 

-11.9b 

J 

-it.  Id 

-17.  U 
-IV.  !•* 
-20.50 
"2i .03 
“2  V • 7  i 
-c3.46 
-24.12 
-<4. 71 
-25.25 
-25.76 
-27.04 
-26 • 38 
-29.  7  3 
“31.00 
-32. OB 
-32* 04 
-32.92 
-32.4V 
-31.13 
-2B. 9  1 
-25. V  5 
-21.81 
-14.07 
-9.66 
-24**6 
-36.22 
-29.15 
-b.oo 
6.6- 
*0.62 
32.  79 
63.37 

52.5  V 
60-66 
61. 

66.5  3 
6«.  V2 
71.0* 
12.  Vi 

74.5  7 
70.05 
62.26 
ol.  b2 
66.  55 
6-.il 
60.30 
76.26 
72.60 
66.  73 
61. 79 
56.6- 
51.92 
A*.  37 
*0  •  9  A 

J6.it) 

31.  ’5 

27.5  7 
23.66 
18. 42 
13.63 

?•  V6 
b.  Vi 
A.  57 
2.67 


ShEBOYGAN,  WISCONSIN 


RUN  BEGINS  AT  HOUR  7  ON  JULY  13*  1972 

STORM  -  BAROMETRIC  PRESSURE  AT  CENTER  OF  LOW  -  990.0  MILLIBARS 

PRESSURE  AT  LARGEST  ENCIRCLING  ISOBAR  •  1006.0  MILLIBARS 
MAXIMUM  PRESSURE  INCLUDED  IN  STORM  »  1010.6  MILLIBARS 

LENGTH  OF  MAJOR  HALF  AXIS  •  1750.0  KILOMETERS 

LENGTH  OF  MINOR  HALF  AXIS  «  875.0  KILOMETERS 

ORIENTATION  OP  MAJOR  AXIS  •  25.0  DEGREES  FROM  NORTH 

LONGSHORE  CURRENT  EQUATION  FROM  FOX  AND  DAVIS*  1972 

SHORE  -  POSITION  COORDINATES  -  X  •  1892.0  Y  >  1112.0  KILOMETERS 

SHORE  LATITUDE  •  *2.  ONSHORE  AZIMUTH  ■  270.  DEGREES 

NEARSHORE  SLOPE  •  0.030  AVERAGE  FETCH  -  200.  KILOMETERS 


HOUR 

X 

V 

XI 

Y 1 

BARO. 

WINP 

SURF. 

ONSH 

ALSH 

EFFtCT. 

WAVE 

WAVE 

breaker 

LSC 

PRESS. 

angle 

W  I  NO 

WIND 

WIND 

WIND 

N 

T 

H 

ANGLE 

VELOC. 

KM 

KM 

RAO 

RAO 

MB 

OEG 

M/S 

M/S 

M/S 

M/S 

M 

SEC 

M 

DEG 

CM/SEC 

7 

585. 

1794, 

-0.25 

-0.49 

1010.53 

49,2 

0.1 

0.1 

0.1 

0.1 

0.00 

0*1 

0.00 

22.9 

1.  74 

8 

o2 1 « 

1793. 

-0.25 

-0,48 

1010.51 

49,3 

0.2 

0.1 

0.1 

0.2 

0.00 

0.1 

0.00 

23.0 

2.17 

9 

65  7. 

1792. 

-0.25 

-0,47 

1010.49 

49,4 

0.2 

0.1 

0.2 

0.2 

0*00 

0.1 

0.00 

23.1 

2.68 

10 

692. 

1791. 

-0.25 

-0.45 

1010.46 

49,6 

0.3 

0.2 

0.2 

0.3 

0.00 

0.2 

0.00 

23.3 

3.27 

n 

728. 

1790. 

-0.25 

-0,44 

1010.42 

49.7 

0.4 

0.2 

0.3 

0.3 

0.00 

0*2 

0.00 

23.4 

3.96 

12 

764. 

1789. 

-0.25 

-0,42 

1010.37 

49,0 

0.5 

0.3 

0.4 

0.4 

0.00 

0*2 

0.00 

23.6 

4.75 

13 

800. 

1788. 

-0.25 

—  0 ,4 1 

1010.32 

50.0 

0.6 

0.4 

0,4 

0.5 

0.00 

0.3 

0.00 

23.  7 

5.65 

14 

842. 

1788. 

-0.25 

-0.39 

1010.24 

50.2 

0.8 

0.5 

0*6 

0.7 

0*01 

0*4 

0*00 

23.9 

6.80 

13 

064. 

1768. 

-0.25 

-0.3# 

1010.14 

50.4 

0.9 

0.6 

0.7 

0-8 

0*01 

0.5 

0.01 

24.0 

8.13 

16 

926. 

1780. 

-0.25 

-0.36 

1010.02 

50.6 

1.2 

0.7 

0.9 

1*0 

0.02 

0.5 

0.01 

24.1 

9.64 

17 

96  9. 

1788. 

-0.25 

-0,35 

1009.87 

50.8 

1.4 

0.9 

1*1 

1*2 

0.03 

0*6 

U.02 

24.2 

11.32 

16 

1011. 

1788. 

-0.2  5 

-0,33 

1009.70 

51.1 

1*7 

1*1 

1*3 

1*5 

0.04 

0.8 

0.03 

24.3 

13.16 

19 

1053. 

1780. 

-0.25 

-0.J1 

1009.49 

51.3 

2.0 

1*3 

1*6 

1*0 

0.05 

0.9 

0.04 

24.4 

15.2  2 

20 

1085. 

1790. 

-0.25 

-0.30 

1009.31 

51.5 

2.3 

1*4 

1*8 

2.0 

0.07 

1.0 

0 .06 

24.5 

17.10 

21 

1116. 

1792. 

-0.25 

-0.29 

1009*12 

51.8 

2.6 

1*6 

2*0 

2*2 

0.08 

1*1 

0.07 

24.6 

16.96 

22 

1147. 

1794. 

-0.25 

-0.2  8 

1008.91 

52.0 

2.9 

1.8 

2.3 

2*5 

0.10 

1.3 

0.09 

24.  T 

20.05 

23 

1179. 

1796. 

-0.26 

-0.27 

1008.67 

52.2 

3.2 

1.9 

2.5 

2.8 

0.12 

1*4 

0.10 

2  h  .  a 

22. 7V 

24 

1210. 

1798. 

-0.26 

-0,25 

1008.41 

52.5 

3.5 

2.1 

2.8 

3.1 

0.15 

1.5 

0.13 

24.9 

24.77 

1 

1242. 

1800. 

-0.26 

-0,24 

1008.12 

52.7 

3.9 

2.3 

3.1 

3*4 

0.17 

1.6 

0.15 

25.0 

26. 80 

2 

1280. 

1806. 

-0.26 

-0.23 

1007.76 

53.1 

4.3 

2*6 

3.4 

3*7 

0.20 

1.8 

0.17 

25.1 

28.91 

3 

1318. 

1813. 

-0.26 

-0.21 

1007.38 

53.5 

4.7 

2.8 

3.8 

4.1 

0.24 

1.9 

0.20 

25.3 

31.06 

4 

1353. 

1119. 

-0.26 

-0.20 

1006.96 

53.8 

5.2 

3.0 

4.2 

4.4 

0.20 

2.1 

0.24 

2S.4 

33.27 

5 

1393. 

1826. 

-0.27 

-0.18 

1006.51 

54.3 

5.6 

3.2 

4.5 

4.6 

0.32 

2.2 

0.2/ 

25.5 

35.44 

6 

1431. 

1832. 

-0.27 

-0.17 

1006.04 

54.7 

6.0 

3.4 

4.9 

5.2 

0.37 

2.4 

0.31 

25.7 

37.57 

7 

1469. 

1839. 

-0.27 

-0.16 

1005.53 

55.1 

6.4 

3*6 

5.3 

5.5 

0.41 

2.5 

0.35 

25.8 

39.63 

8 

1503. 

1862. 

-0.28 

-0.14 

1005.22 

55.7 

6.6 

3,7 

5.5 

5.7 

0.4b 

2.7 

0.3B 

26.0 

41.30 

9 

1341. 

1885. 

-0.29 

-0.13 

1004.92 

56.2 

6.8 

3.8 

5.7 

5.8 

0.50 

2.8 

0.41 

26.2 

42.70 

10 

1576. 

1907. 

-0.30 

-0.12 

1004.62 

56.0 

7.0 

3.0 

5,0 

5.9 

0.54 

2.9 

0.44 

26.3 

43.07 

n 

1612. 

1930. 

-0.31 

-0,10 

1004.33 

57.3 

7.1 

•>*8 

6.0 

6.0 

0.56 

3.0 

0.47 

26.5 

44.84 

w 

1648. 

1953. 

-0.32 

-0.09 

1004.05 

57.9 

7.3 

3,8 

6.1 

6.1 

0.61 

3.1 

0.49 

26.7 

45.62 

13 

1684. 

1976. 

-0.32 

-0.07 

1003. 7V 

50.4 

7.3 

3.0 

6.3 

6.2 

0.64 

3.2 

0.51 

26.8 

46.22 

14 

1702. 

1994. 

-0.33 

-0.07 

1003.76 

50.8 

7.3 

3.8 

6.3 

6.1 

0.67 

3.2 

0.53 

26.9 

46.  B 2 

15 

1721. 

2013. 

-0.34 

-0,06 

1003.73 

59.1 

7.3 

3.7 

6*3 

6.1 

0.69 

3.3 

0.55 

27.0 

47.25 

16 

1739. 

2031. 

-0.35 

-0.05 

100J,T2 

59.4 

7.3 

3.7 

6#.» 

6.1 

0.70 

j.3 

0.56 

27.1 

4  1.54 

17 

1757, 

2050. 

-0.35 

-0.05 

1003.71 

59.  7 

7.2 

3.6 

6.2 

6*0 

0.72 

3.4 

0.57 

27.1 

4  7.  70 

16 

1  776. 

2068. 

-0.36 

“0.04 

1003.71 

60.0 

7.2 

3.6 

6.2 

6.0 

0.73 

3.4 

0.58 

27.2 

47.75 

19 

1794. 

2087. 

-0.3  7 

-0.03 

1003.72 

60.4 

7.1 

3*5 

6*2 

5.9 

0.74 

3.4 

0.58 

27.3 

47.  70 

20 

1811. 

2108. 

-0.37 

-0.03 

1003.78 

60.7 

7.1 

3.i* 

6.1 

5*0 

0.  75 

3.5 

0.58 

27.3 

47.57 

21 

1829. 

2128. 

-0.38 

-0.02 

1003.84 

61,0 

7.0 

3.4. 

6.1 

5.8 

0.75 

3.5 

0.58 

27.4 

47.34 

22 

18 46. 

2149. 

-0.39 

-0.01 

1003.92 

61.3 

6.9 

3*3 

6.0 

5.7 

0.75 

3.5 

0.58 

2  7.4 

47.04 

23 

1863. 

2169. 

-0.40 

-0.01 

1003.99 

61.7 

6.0 

3*2 

6.0 

5.6 

0.75 

3.3 

0.58 

27.4 

46.6  7 

24 

1881  . 

2190. 

-0.41 

-o.oo 

1004.08 

62.0 

6.7 

3*1 

5.9 

5.5 

0.75 

3.5 

0.58 

27.4 

46.22 

1 

1098. 

2210. 

-0.4] 

0.00 

1004.17 

62.3 

6.6 

3.0 

5.8 

5.4 

0.  74 

3.5 

0.57 

27.5 

45.66 

2 

1920. 

2231. 

-0.42 

0.01 

1004.22 

62.0 

6.5 

2.9 

5,8 

5.3 

0. 12 

3.4 

0.55 

27.6 

44.57 

3 

1941. 

2251. 

-0.43 

0.01 

1006.28 

63.3 

6.4 

2*9 

5.7 

5.2 

0.  70 

3  .  ** 

0.53 

27.7 

43.41 

4 

1963. 

2272. 

-0.44 

0.02 

1004.34 

63.8 

6.3 

2*8 

5.7 

5.1 

0.68 

3.3 

0.51 

27.8 

42.19 

5 

1985. 

2293. 

-0.44 

0.03 

1004.42 

64.3 

6.2 

2.7 

5.6 

5.0 

0.65 

3.3 

u.49 

27.9 

40.  VI 

6 

2006. 

2313. 

-0.45 

0,04 

1004.50 

64.9 

6.1 

2.5 

5.5 

4.9 

0.63 

3.2 

0.4b 

20.0 

39.55 

7 

2028. 

2334. 

-0.46 

0.05 

1004.60 

65.6 

6.0 

2*4 

5  .  - 

<*.8 

0.60 

3.2 

0.44 

28.1 

38.13 

WAVE  ENERGY  IN  ThE  BREAKER  ZONE  -  0.168E  09  JOULtS 
TOTAL  LONGSHORE  CURRENT  ENERGY  -  0.1 17E  09  JOulES 
TOTAL  POSITIVE  LONG-SHORE  CURRENT  ENERGY  *  0. 1 1 TE  09  JOULES 
TOTAL  NEGATIVE  LONG-SHORE  CURRENT  ENERGY  «  O.OQOl  00  JOULES 
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5HEBQYGAN.  WISCONSIN 


HUN  BEGINS  AT  hOU»  7  ON  JULY  lo»  1972 

STORM  -  BAROMETRIC  PRESSURE  AT  CENTER  OF  LOW  •  1002.0  MILLIBARS 
PRESSURE  AT  LARGEST  ENCIRCLING  ISOBAR  -  1016.0  MILLIBARS 
MAX [mum  PRESSURE  INCLUOEO  IN  STURM  «  1028.0  MILLIBARS 

LENGTH  OF  MAJOR  HALF  AXIS  •  1050.0  KILOMETERS 

t  £  NOTH  OF  MINOR  HALF  AXIS  -  700.0  KILOMETERS 

ORIENTATION  OF  MAJOR  AXIS  ■  20.0  DEGREES  FROM  NORTH 

LONGSHORE  CURRENT  EQUATION  FROM  FOX  AND  DAVJS.  1972 

SHORE  -  POSITION  COORDINATES  -  X  •  1892.0  V  -  1112.0  KILOMETERS 

SHORE  LATITUDE  •  42.  ONSHORE  AZ I MjTh  ■  270,  DEGREES 

NEARSHORE  SLOPt  •  0.030  AVERAGE  FETCH  ■  200.  KILOMETERS 


hour 

X 

Y 

Al 

Y 1 

BARO. 

WIND 

SURF. 

ONSH 

ALSH 

EFFECT. 

WAVE 

WAvE 

BREAKER 

LSC 

PRESS. 

angle 

WIND 

WIND 

WIND 

WIND 

M 

T 

M 

ANGLE 

VtLOC. 

KM 

KM 

RAD 

RAD 

MB 

DEG 

M/S 

M/S 

M/S 

M/S 

M 

SEC 

M 

DEG 

CM/StC 

7 

176. 

286. 

0.52 

-1.08 

1018.02 

29.6 

0.0 

0.0 

0.0 

0.0 

0.00 

0.0 

0.00 

O.o 

0.00 

0 

2 12. 

296. 

0.51 

-1.06 

1018.02 

29.5 

0.0 

0.0 

0.0 

0.0 

0.00 

0.0 

0.00 

0.0 

0.00 

9 

247. 

306. 

0*51 

-1.04 

1018.02 

29.3 

0.0 

o.o 

0.0 

0*0 

0.00 

0.0 

0,00 

o.o 

o.oo 

10 

283. 

313. 

0.50 

-1.02 

1018.02 

29.1 

0.0 

0.0 

0.0 

0*0 

0.00 

0.0 

0.00 

0.0 

0*00 

11 

319. 

323. 

0.*9 

-0.99 

1018.02 

29.0 

0.0 

0.0 

0.0 

0.0 

0.00 

0.0 

0.00 

0.0 

o.oo 

12 

354, 

335. 

0.49 

-0.97 

1018.02 

2B.8 

0.0 

0.0 

0.0 

0.0 

0.00 

0.0 

0.00 

0.0 

o.oo 

13 

390. 

345, 

0.48 

-0.93 

1018.02 

28.6 

0.0 

0.0 

0.0 

0.0 

0.00 

0.0 

0.00 

0.0 

o.  ct 

14 

423. 

336. 

0.48 

-0.93 

1016.02 

28.4 

9.0 

0.0 

0.0 

0.0 

0.00 

0.0 

0.00 

0.0 

O.uv 

13 

459. 

367. 

0.47 

-0.90 

1016.02 

28.3 

0.0 

0.0 

0.0 

0.0 

0.00 

0.0 

0.00 

0.0 

O.oo 

16 

49*  • 

3  77. 

0.46 

-0.08 

1018.02 

28.1 

0.0 

0.0 

0.0 

0.0 

0.00 

0.0 

0.00 

0.0 

o.oo 

17 

329. 

388. 

0.45 

-0.06 

1016.02 

27.9 

0.0 

0.0 

0.0 

o.o 

0.00 

0.0 

0.00 

c.o 

9.0  0 

16 

363. 

399. 

0.45 

-0.04 

1018.02 

27,7 

0.0 

0.0 

0.0 

0.0 

0.00 

0.0 

0.00 

0.0 

0.09 

19 

390. 

410.. 

0.44 

-0.62 

1018.01 

27.5 

0.0 

0.0 

0.0 

0.0 

0.00 

0.0 

0.00 

0.0 

0.00 

20 

630. 

*26. 

0.43 

-0.00 

1C18.01 

27,4 

0.0 

0.0 

0.0 

0.0 

0.00 

0.0 

0.00 

C.u 

o.oo 

21 

663. 

*38. 

0.42 

— u  .78 

1018.00 

27.3 

0.0 

0.0 

0.0 

0.0 

0.00 

0.0 

0.00 

0.0 

0.00 

22 

695. 

-32. 

0.41 

-0.73 

101 7.99 

27.  1 

0.1 

0.1 

0.0 

o.l 

0.00 

0.0 

0.00 

0.0 

0,00 

23 

726. 

*60  • 

0.41 

-0.73 

1017.97 

27.0 

0.2 

0.1 

0.0 

0.1 

0.00 

0.1 

0.00 

13.7 

1.82 

2*» 

760. 

*80. 

0.40 

-0.71 

1017.94 

26.9 

0.2 

0.2 

0.1 

0.2 

0.00 

0.1 

0.00 

13.6 

2.59 

1 

793. 

494. 

J.  39 

-0.69 

1017.91 

26.8 

0.3 

0.3 

0.1 

0.3 

0.00 

0.2 

0.00 

13.6 

3,49 

2 

041  • 

518. 

0.37 

-0.66 

1017.03 

26.7 

0.5 

0.5 

0.2 

0.5 

0.00 

0.3 

0.00 

13.7 

5.09 

3 

088. 

542, 

0.36 

-0.63 

1017.72 

26.5 

0.8 

0.7 

0.3 

0.8 

0.01 

0.4 

0.01 

13.7 

7.13 

4 

936. 

363. 

0.34 

-0.60 

1017.55 

26.4 

1.2 

1.2 

0.5 

1.2 

0.02 

0,6 

0.02 

13.6 

9.52 

5 

904. 

509. 

0.33 

-0.57 

1017.33 

26.3 

1  •  7 

1.5 

0.7 

1.6 

0.04 

0.6 

0.0* 

13.6 

12.  50 

6 

1031. 

613. 

0.31 

-0.54 

1017.03 

2  6.1 

2.2 

2.0 

i.o 

2.2 

0.0b 

1.0 

0.06 

13.5 

15.47 

7 

1079. 

63  7. 

0.30 

-0.51 

1 0  1 6  . 64 

26.0 

2.9 

2.6 

1.3 

2,8 

0.10 

1.2 

0.10 

13.5 

19. 70 

0 

1130. 

609. 

0.26 

-C  .48 

1016.07 

2  7.0 

3.8 

3.4 

1.7 

3.7 

0.15 

4.5 

0.15 

14.0 

24.05 

9 

1102. 

740. 

0.23 

-C.4S 

1015.3? 

26.2 

4.8 

4.2 

2*2 

*.6 

0.22 

1.8 

0.22 

14.5 

30.55 

10 

1233. 

79 2. 

0.20 

-C.41 

102*. 53 

29.5 

3.7 

5.0 

2,0 

5.5 

0.30 

2.2 

0.30 

15.2 

36.6* 

11 

1284. 

044. 

0.17 

-0.38 

1013.58 

31.0 

6.7 

5.7 

3.4 

6*4 

0.40 

2.5 

o.*o 

16.0 

42.93 

12 

1336. 

095. 

0.13 

-3.33 

1012.51 

32.8 

7.6 

6.3 

4.1 

1*2 

0.51 

2*6 

0.50 

16.8 

49.  U 

13 

1307. 

9*7, 

0.10 

-0.32 

1011.30 

34.9 

8.3 

6.0 

4.7 

7.8 

0.61 

3.1 

O.bv. 

17.0 

54.92 

14 

1421. 

909. 

0.07 

-0.29 

1010.61 

37. C 

8.7 

5.9 

5.2 

S*1 

0.71 

3.3 

0.69 

18.7 

59.59 

15 

1456. 

1030. 

0.05 

-0.27 

1009 .85 

39.4 

9.0 

6.9 

5.7 

0*3 

0.00 

3.5 

0.7b 

19.7 

63.31 

2  o 

1*90. 

207?. 

l.JJ 

-0.25 

1009,11 

42,0 

9.2 

.  .  d 

6.1 

0.4 

0.88 

3.7 

0.82 

20.9 

55.9  3 

IT 

1324. 

1113. 

-0. 00 

-0.23 

1008.41 

43.0 

9.3 

6.6 

0.0 

6  •  H 

0.94 

3.0 

0.87 

22.1 

6  7.26 

10 

1359. 

1153. 

-0.02 

-C.21 

1007.77 

40.0 

9.3 

6.2 

6,9 

6.3 

1.0G 

3.9 

0.89 

23.3 

67.3* 

19 

1393. 

1196. 

-0,93 

-0.18 

1007.20 

51.0 

9.2 

5.8 

7.2 

0.1 

1.0* 

-.0 

0.91 

24.  3 

66.36 

20 

1636. 

12  33. 

-0.07 

-0.16 

1006.4* 

53.0 

i.O 

5.3 

7.2 

T.0 

1.06 

4.1 

0  .  >0 

25.3 

6*.  55 

21 

.600. 

12  70. 

-Q.lu 

-0.13 

1005.79 

56.4 

8.6 

4.7 

7.2 

7.3 

1.07 

4.1 

0.89 

26. . 

0«t.  OA 

22 

1723. 

1  307. 

-0.12 

-0.10 

1003.28 

59.0 

6.2 

*•2 

T.v 

6.8 

0.90 

3.9 

0.78 

26.9 

56.33 

23 

1766. 

1343. 

-0.1* 

-0.07 

1004.92 

61.6 

7.7 

3.6 

6.8 

6.3 

0.86 

3.  r 

0.6  T 

27.6 

50*1* 

26 

1010. 

1360. 

-0.17 

-0.05 

100*. 73 

65.1 

7.2 

3.0 

6.5 

5.0 

0.75 

3.5 

0.53 

20.5 

*1.14 

1 

I8t  3. 

1417. 

-0.19 

-0.02 

1004.71 

TO. 9 

6.6 

2.1 

6.2 

5.; 

0.63 

3.2 

o.*o 

2  9.6 

33.34 

2 

1002. 

1304. 

-0.24 

-..CO 

1005.80 

77.3 

6.6 

1.4 

6.5 

4.9 

0.61 

3.2 

9.32 

30.5 

<6.  02 

3 

1911* 

1590. 

-0.30 

0.01 

1007.09 

83.1 

6.5 

0.7 

6.4 

4.5 

0.55 

3.0 

0.22 

30.8 

17.  42 

4 

1940. 

1677. 

-0.33 

0.03 

1000.52 

87.3 

6.2 

0.2 

6.2 

4.2 

0.48 

2.8 

o.ll 

20.9 

9.95 

5 

1970, 

17  6*. 

-0.41 

0.04 

1010.00 

90.9 

5.1 

-0.0 

5.0 

3.0 

0.40 

2.6 

0.06 

30.  1 

5.65 

6 

1999, 

1630, 

-0.46 

0.06 

1011.47 

93.6 

5.3 

-0.3 

5.2 

3.4 

0.32 

2.4 

0 . 0  v 

30.  * 

10.03 

7 

2026. 

1917. 

-0.32 

0.00 

1012.67 

95.0 

4.6 

-0.4 

4.6 

2.9 

0.2* 

2.1 

0.09 

30.1 

11.13 

9 

2057, 

1900, 

-0.33 

0.09 

1 01 3.65 

95.0 

4.3 

-0.* 

4.2 

2.7 

0.20 

1.9 

0.07 

29.  V 

10,3  1 

9 

2043. 

2039. 

-0.38 

0.09 

101*. 37 

95.9 

3.8 

-0.4 

3.8 

2.4 

0.17 

1.0 

0.06 

29,0 

4,43 

10 

203*. 

2090. 

-0.62 

0.10 

1015.03 

95.9 

3.4 

-0.3 

3.4 

2.1 

0.13 

1.6 

v  .  05 

29.6 

0  .  *  * 

11 

206  3  . 

21*1. 

-0.43 

0.  10 

1013.61 

95.9 

3.0 

-0.3 

3.0 

1.9 

0.10 

1  .* 

o.  03 

29.3 

7#*l 

12 

:  o  t  i . 

2192. 

-0.66 

0.11 

1016.12 

95. V 

2.5 

-0.2 

2.5 

1.6 

0.07 

1.2 

0  «  o2 

2 V.  ' 

4.3* 

13 

2000. 

2243. 

— u.  n 

o.n 

1016.36 

96,0 

2.1 

-0.2 

2.1 

1.3 

0.05 

1  .0 

0.02 

26.9 

5.30 

14 

2099. 

2288. 

-J.  7* 

0.13 

1016.89 

9  7,0 

1.7 

-0*2 

1.7 

1.1 

0,03 

0.0 

0.01 

20.  7 

*  .  7y 

13 

2117. 
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RUN  8EGJNS  AT  hOU«  4  ON  DECEMBER  14.  1973 
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longshore  Current  equation  from  fox  and  uavis*  1972 
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MIXED  SEMIDIURNAL  TIDE  -  FORM  NUMBER  IS  0.90 
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